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I 


Summary 


In  order  to  clarify  the  effects  of  vibrational  excitation  on  shock-wave  transi¬ 
tions  of  weak,  spherical  N-waves,  which  were  generated  by  using  sparks  and  exploding 
wires  as  sources,  the  compressible  Navier-Stokes  equations  were  solved  numerically, 
including  a  one-mode  vibrational-relaxation  equation.  A  small  pressurized  air-sphere 
explosion  was  used  to  simulate  the  N-waves  generated  from  the  actual  sources.  By 
employing  the  random-choice  method  (RCM)  with  an  operator-splitting  technique,  the 
effects  of  artificial  viscosity  appearing  in  finite-difference  schemes  were  eliminated 
and  accurate  profiles  of  the  shock  transitions  were  obtained.  However,  a  slight  ran¬ 
domness  in  the  variation  of  the  shock  thickness  remains.  It  is  shown  that  a  computer 
simulation  is  possible  by  using  a  proper  choice  of  initial  parameters  to  obtain  the 
variations  of  the  N-wave  overpressure  and  half-duration  with  distance  from  the  source. 
The  calculated  rise  times  are  also  shown  to  simulate  both  spark  and  exploding-wire 
data.  It  was  found  that,  in  addition  to  the  vibrational-relaxation  time  of  oxygen, 
both  the  duration  and  the  attenuation  rate  of  a  spherical  N-wave  are  important  factors 
controlling  its  rise  time. 

The  effects  of  the  duration  and  the  attenuation  rate  of  a  spherical  N-wave  on  its 
rise  time,  which  are  designated  as  t kv  tl-vave  effect  and  ik*  ncnsiat lonary  effect, 
respectively,  are  discussed  in  more  detail  pertaining  to  Lighthill's  analytical  solu¬ 
tions  and  the  RCM  solutions  for  nonstationary  plane  waves  and  spherical  N-waves'k.  It 
is  also  shown  that  the  duration  and  the  attenuation  rate  of  a  spherical  N-wave 
affected  by  viscosity  and  vibrational  nonequilibrium,  so  that  they  can  deviate  from 
the  results  of  classical,  linear  acoustic  theory  for  very  weak  spherical  waves. 
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Ihe  pressure  waves  generated  by  supersonic 
transport  aircraft  (SST)  and  from  explosions  in  air 
are  often  observed  as  weak  X -waves  far  from  the 
source.  Such  pressure  waves  are  heard  as  sonic 
booms.  The  loudness  of  these  waves  depends  on 
their  maximum  overpressures  and  rise  times  (Ref.  1). 
The  N- waves  with  short  (microseconds )  rise  times 
are  perceived  as  louder  and  more  startling  than 
the  ones  with  long  (milliseconds)  rise  times.  As 
a  consequence ,  \-wave  rise  times  were  investigated 
extensively  for  SSI'  sonic  booms  and  for  explosions 
in  air  (Refs.  2-4).  However,  the  observed  SSI 
rise  times  were  often  found  to  be  larger  than  those 
wh  i  e  li  w  e  re  e  s  t  i  rna  t  ed  f  rom  classical  t  h  o  o  ry  f  o  r 
viscous  shock  structures  of  steady,  plane  waves, 
derived  by  Taylor  (Ref.  3i.  A  recent  review  of 
this  matter  may  he  found  in  Ref.  o. 

This  d i sc  repancy  wa  s  a  1 1  r i but  ed  ma inly  to  t  h  e 
effects  of  atmospheric  turbulence  ikefs.  '-101 , 
and  real -gas  effects  arising  from  the  vibrational 
excitation  of  the  oxygen  and  nitrogen  air  molecules 
(Refs.  11,  1  -2 )  .  However,  the  decisive  factor  for 
this  increased  rise  time  was  still  in  question. 

The  re  we re  d i f  f i cu 1 1 ies  in  p ro  v i d i n g  c  or  re 1  a  t i on  s 
between  the  observed  and  analytically  estimated 
rise  times,  owing  to  a  lack  of  information  regard¬ 
ing  the  ambient  temperature,  humidity  and  air 
turbulence.  Such  quantities  are  not  always  readily 
available.  It  was  therefore  necessary  to  carry  out 
some  simulation  experiments  under  controlled  con 
Jit  ions  where  known  atmospheric  conditions  could 
be  obtained. 

Hoi  st -Jensen  iRef.  t>  i  was  able  to  generate 
well -formed  weak  spherical  N-wavos  by  using  sparks 
or  exploding  wires  as  a  source  in  a  still -air  dome, 
usually  used  for  air-cushion  experiments  (Ref.  13). 
In  this  manner  he  wanted  to  clarify  the  vibrational 
effects  on  the  rise  time  of  SSI  N  waves.  He  found 
that  the  observed  rise  times  were  much  shorter 
than  the  rise  times  estimated  from  the  analysis  of 
plane,  fully  J i spersed  waves  iRef.  Id).  The  results 
could  not  be  explained  h\  any  existing  analysis. 

The  object  of  this  report  is  to  j>r.  .  ;  le  a  theoretic 
al  basis  for  explaining  Holst  Jensen1'*  data,  which 
will  be  outlined  in  vect ion 

The  processes  involved  in  the  generation  of 
N-wavos  by  exploding  sparks  and  wires  are  very  com 
plex  and  are  not  readily  predicted.  Consequently, 
it  is  necessary  to  assume  a  reasonable  source  model 
in  order  to  simulate  the  explosions.  In  this  paper 
it  is  assumed  that  the  expanding  plasma  can  be 
simulated  by  a  pressurised  sphere  of  small  radius 
at  room  temperature.  The  computer  simulation 
requires  adjusting  the  radius  of  the  pressurized 
sphere  and  the  imaginary  diaphragm  pressure  ratio 
to  fit  the  experiments  for  maximum  overpressure 
and  hal f-durat ion  of  the  N-wave  with  distance  from 
the  source.  It  is  then  possible  to  determine  the 
•  //*  "  '  '•  initial  energy  of  the  source.  The  latter 
is  of  academic  interest  as  it  is  not  possible  to 
determine  the  actual  energy  release  from  the  voltage 
and  capacitance  of  the  discharge  without  a  great 
deal  of  additional  time-dependent  measurements. 

The  nonstationary,  spher ica I  -symmet ric  Navicr- 
Stokes  equations  were  solved  numerically,  including 
the  equation  of  one-mode  vibrational  relaxation  for 
explosions  of  pressurised  phen^  in  atmospheric 


air.  An  operator  spl j t t ing  technique  was  used  in 
which,  at  the  first  stage  of  calculation,  the  solu¬ 
tions  for  inviscid,  fro -on  flow  were  obtained  by 
applying  the  Random  Choice  Method  (RCM)  and  then 
the  effects  of  viscosity  and  vibrational  nonequili¬ 
brium  were  evaluated  by  using  an  explicit  finite- 
difference  method. 

The  RCM  is  a  numerical  method  which  was  devel¬ 
oped  hv  til  i  mm  (Ref.  14},  Chonn  'Ref.  13)  and  Sod 
(Ref.  lb j  for  flow  problems  including  shock  waves. 

In  this  method,  a  Riemann  problem  is  solved  for 
each  spatial  mesh  at  each  time  step  and  then  one 
of  its  solutions  is  chosen  at  random  as  a  solution 
for  the  next  time  step  by  using  a  random  sampling 
techni  ,ue .  It  is  the  great  merit  of  this  method 
that  shock  waves  and  contact  surfaces  can  be  ex 
pressed  as  discontinuous  surfaces  without  smearing 
arising  f rom  art i f ic ia 1  viscosities  inherent  in 
all  finite-difference  methods.  Hus  is  the  main 
reason  for  adopting  the  RO!  for  the  present  analysis. 
The  algorithm  is  based  on  a  program  developed  by 
Saito  and  Glass  IRef.  1").  The  application  of  the 
operator-splitting  technique  for  analysing  the 
Nav ier-Stokes  equations  was  first  introduced  by 
MacCormack  (Ref.  18  .  In  his  analysis,  the  inviscid 
solutions  were  obtained  using  a  characteristic 
method.  Recently,  Satofi,ln  and  Shimizu  (kef.  \9 ) 
have  tried  to  solve  the  N  ier-Stokes  cqua* }ons 
for  a  shcck-ti.  >e  problem  ny  applying  the  RCM  with 
an  operator-splitting  technique.  In  the  present 
analysis,  the  RCM  with  an  operator-splitting  tech¬ 
nique  was  extended  to  include  vibrational  relaxation 
effects  for  spher ica 1 ly- symmet ric  waves. 

It  will  be  shown  subsequently  nat  the  rise 
times  of  weak,  spherical  N- waves  .crated  by  sparks 
and  exploding  wires  are  seriously  affected  by  two 
factors  which  never  appear  in  steady  plane  waves. 
These  are  designated  as  an  V-:/:-’-  effect  and  a 
’•  ■  ;V’  "*  respectively.  The  N-wave 

effect  means  that  the  rise  times  of  weak  N-waves 
are  affected  by  the  expansion  of  the  flow  immedi- 
a t e 1 >  behind  the  shock  front.  The  nonstat ior ary 
effect  means  that  the  rise  times  of  weak  shock 
waves  respond  to  chai  cos  in  shock  strength  so  slowly 
that  their  transient  behaviours  must  be  considercJ. 

I  he  fundamental  analytical  ideas  about  these 
effects  were  provided  by  l.  ighthill  (Ref.  20)  for 
both  viscous  N-waves  and  impulsively-generated 
viscous  plane  waves.  In  Section  3,  his  results 
are  re -examined  for  use  in  the  present  study. 

In  order  to  consider  the  effects  of  vibrational 
excitation  of  oxygen  and  nitrogen  air  molecules,  the 
papers  of  Polyakova  et  al  (Ref.  21)  and  Johannsen 
and  Hodgson  (Ref.  12)  for  plane,  dispersed  waves 
are  also  re-examined  in  Section  3,  and  an  approxi¬ 
mate  relation  is  rived  for  the  rise  time  of  a 
fully  or  partly-dispersed  wave,  furthermore,  the 
modified  Taylor  and  I. ighthill  solutions  for  fully- 
dispersed  waves  are  discussed. 

In  Sections  4.1  and  4.2,  the  basic  equations 
and  the  numerical  method  of  solution  are  described. 

In  Section  4.3,  to  validate  the  method  of  solution 
for  nonstationary  shock  transitions,  RCM  solutions 
for  nonstationary  viscous  and  dispersed  plane  waves 
arc  compared  with  analytical  solutions  described 
in  Section  3.  As  for  solutions  for  spherical  waves 
(Section  4.4),  some  numerical  results  for  weak 
spherical  N-waves  in  air  are  presented  for  the 
following  five  cases:  (i)  formation  of  N-waves  in 


the  near-field  of  a  pressurized  sphere,  f i i )  com¬ 
parison  between  perfect-inviscid,  perfect-viscous, 
real -inviscid  and  real-viscous  solutions,  iiii) 
effects  of  vibrational  relaxation  time  or  ambient 
temperature  and  humidity,  (iv)  effects  of  N-wave 
duration  or  radius  of  pressurized  sphere,  and  (v) 
effects  of  nitrogen  vibrational  relaxation.  The 
observed  rise  times  of  spark  and  exploding-wire 
generated  N  -waves  are  also  compared  with  those 
obtained  from  the  analytical  simulations. 

In  this  report,  the  usual  definition  of  rise 
time  is  followed,  and  is  taken  as  the  time-interval 
for  the  overpressure  to  vary  from  10 to  90"»  of  its 
peak  value.  This  definition  is  quite  arbitrary  and 
is  especially  useful  for  actual  SST  signatures,  as 
discussed  in  Ref.  6.  Figure  1.1  illustrates  the 
definition  of  an  N-wave  rise  time  tr  and  its  half¬ 
duration  tj.  Figure  1.2  also  illustrates  the 
definition  of  a  plane-wave  rise  time  tr.  The 
corresponding  shock  thickness  2.x  and  half-duration 
length  xj  may  approximately  be  given  by 

2x  =  a ,  t  ,  x,  =  a.t. 

1  r  d  Id 

where  aj  is  the  undisturbed  speed  of  sound,  since  we 
consider  only  very  weak  waves. 

2.  SPARK  AND  FXPl.OPING-WIRF  DATA 

In  this  section,  the  spark  and  exploding-wire 
experiments  which  were  carried  out  by  Holst -Jensen 
(Ref.  b)  and  the  resulting  data  are  summarized. 

The  purpose  of  these  experiments  was  to  generate 
weak,  fully-developed  N-waves  with  overpressure 
below  100  Pa  in  air,  which  would  have  interference 
free  shock  fronts.  This  was  accomplished  by  using 
sparks  and  exploding  wires.  The  dome  containing 
the  UTIAS  air  cushion  vehicle  (ACV)  circular  track 
facility  (Ref.  13)  was  used  as  a  still-air  reservoir 
for  part  of  the  experiments.  Its  major  internal 
diameter  is  about  42. “m.  This  provided  waves  free 
from  interference  with  walls  and  other  objects. 

For  detecting  weak  shocks  in  the  overpressure 
range  5-100  Pa,  a  condenser  microphone  was  used 
[Bruel  &  Kjaer  4135  free  field  6.3  mm  (1/4  in)  dia|. 
Amplification  of  the  microphone  signal  was  provided 
by  a  preamplifier  B8K  2619.  The  response  of  the 
microphone  system  was  tested  in  the  III  IAS  Travelling- 
Wave  Sonic-Boozi  Simulator  (Ref.  22).  When  measuring 
without  its  protective  grid  at  zero  angle  of  inci¬ 
dence,  the  microphone  has  an  approximate  minimum 
rise  time  tr  =  2.9  usee.  The  oscilloscopes  used 
were  Tektronix  types  555  and  535  with  a  type  I) 
plug-in  that  has  a  bandwidth  better  than  300  KHz. 

The  microphone  was  calibrated  with  a  BfiK  pistophone 
type  4220,  which  gives  a  sound  pressure  level  at 
250  Hz  of  124  dB. 

In  the  first  series  of  experiments,  sparks  were 
used  as  a  source  of  N-waves.  The  sparks  were  gener¬ 
ated  by  the  energy  released  from  a  charged  ’’.5  uF 
capacitor.  The  maximum  charging  voltage  was  8  KV 
and  the  discharge  device  was  a  thyratron.  A  micro¬ 
phone  was  placed  ahead  of  the  measuring  microphone 
in  parallel  to  get  the  trigger  signal  for  the 
oscilloscope.  The  source  and  microphone  were  set 
up  at  1.8m  above  the  floor  to  avoid  interference 
from  reflected  signals. 


Fairly  extensive  measurements  were  done  by 
using  sparks  at  temperatures  of  273-277  K  and 
relative  humidities  of  50-73”,,.  Five  source-receiver 
distances  (4.1m,  4.9m,  9.8m,  15. bm  and  21.6m)  were 
employed  with  four  different  charging  voltages  of 
4.4  KV,  5.0  KV,  5.4  KV  and  6.0  KV.  This  series  of 
measurements  is  termed  Series-I.  Another  series  of 
measurements  (Series-1 1)  was  also  done  at  a  temper¬ 
ature  of  289  K  and  relative  humidity  of  503.  for  the 
distance  range  of  11. 8-19. 0m  and  a  charging  voltage 
of  4 .4  KV. 

Fxploding  wires  were  used  to  produce  N-waves  by 
replacing  the  resistor  in  the  sparK  circuit  by  a 
thin  nickel  wire  0.125  mm  dia  and  optimum  length  of 
5  cm.  The  sudden  discharge  of  energy  vaporized  the 
wire.  The  expansion  of  the  metal  vapour  generated 
an  N-wave  in  the  far  field.  The  measurements  were 
done  at  two  conditions  for  Series-111  (T|  =  2"  k, 

Rll  =  753,,  r  =  6.7m,  12.8m,  24.3m,  S  =  4.6  KV,  6.0 
KV),  and  Series -IV  (Tj  =  280  K,  Rll  =  8'. 53,,  r  = 

24.3m,  29.3m,  S  -  4.6  KV,  6.0  KV ) ,  where  T]  is  the 
room  temperature,  Rll  the  relative  humidity,  r  the 
distance  from  the  source  and  S  the  charging  voltage. 

The  vibrational  relaxation  times  for  oxygen  and 
nitrogen  were  evaluated  by  using  the  empirical  re¬ 
lation  obtained  from  the  absorption  of  sound  waves 
by  Bass  and  Shields  (Ref.  23),  as  tabulated  in 
Table  2.1.  The  vibrational  relaxation  time  at  room 
temperature  strongly  depends  on  the  absolute  humidity 
of  the  atmosphere,  as  water  molecules  significantly 
reduce  its  value. 

Representative  oscillograms  from  sparks  and 
exploding  wires  are  shown  in  Fig.  2.1.  It  can  be 
seen  that  both  a  spark  and  an  exploding-wire  source 
make  it  possible  to  produce  well-established  N-waves 
far  from  the  source.  In  the  exploding-wire  experi¬ 
ments,  the  N-waves  were  much  cleaner  than  those 
generated  by  a  spark,  especially  with  regard  to  the 
rear  shock.  It  was  found  that  the  wire  length  1. 
plays  a  significant  role  in  shaping  the  rear  shock 
pressure  profile.  After  testing  several  wire 
lengths,  a  wire  length  1.  =  5.0  cm  proved  to  generate 
the  most  symmetrical  N-waves,  and  was  used  in  all 
subsequent  runs.  The  microphones  were  set  up  normal 
to  the  wire  to  minimize  any  line-source  effect. 

In  Figs.  2.2  -  2.4,  the  maximum  (peak)  over¬ 
pressure  (>p)max.  the  half-duration  tj  and  the  rise 
time  tr  are  plotted  against  the  distance  from  the 
source  r.  Figure  2.5  shows  plots  of  tr  vs  ( ,'.p )  ma  x . 
The  data  for  different  series  are  represented  by 
different  symbols,  which  arc  common  through  Figs. 

2.2  -  2.5.  For  the  Series-I  experiment,  the  data 
are  plotted  only  for  S  =  4.4  KV  and  6.0  KV  to  avoid 
confusion . 

In  Fig.  2.2,  the  lines  indicate  the  curves  of 
(■‘Plmax  '  r  n  ’  "'h'Fh  are  drawn  from  the  arbitrary 
points  to  fit  the  experimental  data,  where  n  is 
termed  the  decay  index  of  maximum  overpressure. 

The  solid  and  broken  lines  correspond  to  the  curves 
for  n  =  1  and  1.4,  respectively.  For  100  -  ( Ap ) max 
20  Pa  both  spark  and  exploding-wire  data  show  that 
maximum  overpressures  decay  nearly  inversely  pro¬ 
portional  with  distance  from  the  source,  as  estimated 
from  linear-acoustic  theory.  On  the  other  hand,  the 
spark  data  show  that  the  decay  index  increases  below 
20  Pa.  This  deviation  from  linear-acoustic  theory 
can  be  attributed  to  real-gas  effects  arising  from 


I 
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vibrational  excitation  of  oxygen  i see  Section  4.4). 
It  is  noted  that  the  same  input  energy  does  not 
result  in  the  same  decay  of  l-iplnuix  *or  different 
energy  sources.  The  exploding-wi re  source  makes 
for  a  stronger  explosion  in  air  than  the  spark 
source  for  the  same  discharge  voltage.  It  should 
also  he  noted  that  the  overpressure  decays  are 
different  for  the  different  series  of  spark  exper¬ 
iments  despite  the  same  discharge  voltage. 

In  figs.  2.3  -  2.5,  the  broken  lines  indicate 
the  tendency  of  the  experimental  data.  The  half- 
du ra t i on  t  j  i nc  rea s e s  with  r .  Th e  du ra  1 1 on  s  f o r 
the  exploding-w ire  experiment  (85-155  ..sec)  are 
longer  than  those  for  the  spark  experiments  (50 -"*5 
..sec).  Ihe  rise  times  tr  also  increase  with  r, 
while  the  maximum  overpressure  decreases  with  r. 
it  should  be  noted  from  Jig.  2.5  that  the  rise 
times  tr  are  different  for  the  different  series  of 
experiments  and  supply  voltages  at  the  same  maximum 
overpressure . 


5 .  50M1  ANALYSES  IQR  WHAk  SHOCK  TRANSITIONS 

In  this  section,  some  analytical  solutions  for 
weak  shock  transitions  are  reviewed  and  discussed 
in  connection  with  the  spark  and  exploding-wire 
data,  which  were  shown  in  Section  2.  In  Sections 
5.1  -  5.5,  some  analytical  solutions  for  viscous- 
shock  transitions  are  shown  in  cases  of  steady 
planar  waves,  quasi -stat ionary  N -waves  and  nonsta¬ 
tionary  planar  waves,  respectively.  The  analytical 
solution  for  steady  planar  waves  was  derived  by 
I  ay  lor  (Ref.  3),  and  will  be  designated  as  the 
lay  lor  solution  or  the  Taylor  shock  transition. 

The  analytical  solutions  for  quasi -stat ionary 
N-waves  and  r.onstat  ionary  planar  waves  were  defined 
by  I.ighthiM  (Ref.  JO),  and  will  he  designated  as 
the  lighthill  solutions,  or  the  l.ighthill  \-wave  and 
the  l.ighthill  shock  transition,  respectively.  In 
Section  5.4,  solutions  for  dispersed  waves  with 
vibrational  excitation  are  shown  for  a  steady 
plane  wave,  and  an  approximate  expression  is 
derived  for  the  rise  time  of  a  fully  or  partly- 
dispersed  wave.  The  Taylor  and  l.ighthill  solutions 
are  extended  to  dispersed  waves  with  vibrational 
relaxation  by  using  a  bulk-viscosity  concept,  and 
the  extended  solutions  will  be  designated  as  the 
modified  Taylor  solution  and  the  modified  l.ighthill 
solution,  respectively.  Some  insight  is  also  given 
into  the  structures  and  rise  times  of  weak  spherical 
\ -waves . 


5.1  l :_las s  i c al  Taylor  Plane  Shock  -Wave  Trans  it  ions 

In  the  following  three  sections,  Sections 
5.1  5.5,  the  ’ •  «.•  or  ;V  :;•* -shock  transitions 

are  considered,  where  the  vibrational  node  of 
molecular  internal  energy  is  assumed  to  be  /V  £•  *: . 
Viscous,  steady  shock  waves  are  formed  as  a  result 
of  a  balance  between  the  wave  *  form- steepen ing 
tendency  due  to  the  finite-amplitude  compression 
(convection)  effects  and  the  wave  form-easing 
tendency  due  to  the  vi scons -d if fus ion  effects. 

This  balancing  determines  the  thickness  of  a  steady 
shock  wave  and  depends  on  the  shock  strength. 


v  , 


i 

-  (1 


1  -  +  1  )v  .ix-il^t ) 
exp - y? - 


(3.1  ) 


for  a  shock  w.  travelling  with  steady  profile  at 
a  constant  .>pecd  Us ,  where  v  =  flow  velocity 
re ! a t i ve  to  the  g round ;  v ’  -  f 1  ok  ve 1 oc i t v  a t 
x  •  •  ,  :  =  ratio  of  specific  heats,  x  --  distance, 

t  =  time,  -  dif'fusivitv  of  sound,  defined  by 


where  =  kinematic  viscosity,  ..  =  viscosity, 

..r  =  bulk  viscosity  due  to  rotational  relaxation, 
I’r  =  Prandt  1  number.  All  the  thermodynamic  and 
transport  coefficients,  ,  ,  ,  ,.r  and  Pr,  may 

be  assumed  to  be  constant  throughout  the  flow, 
since  the  shock  waves  are  weak.  The  original 
Taylor  solution  did  not  include  the  bulk  viscosity 
due  to  rotational  relaxation  as  it  appears  in  Tq. 
(5.2).  However,  in  the  present  paper,  the  term 
rr«’cr  .-*■  ‘s'-’  v  is  used  when  it  includes  only  the 
effects  of  rotational  relaxation  in  order  to  dis¬ 
tinguish  from  the  -  :  .**  .*  *  <*  *  •;  which 

includes  both  the  effects  of  rotational  and  vibra¬ 
tional  relaxation. 

l-'rom  the  weak-wave  assumption,  we  have 


v/aj  ?  i  ‘.pm  fpj ) 


(  a  .  5  I 


where  'p  is  the  overpressure  (.’p  =  p  -  pj )  ;  aj, 
the  undisturbed  speed  of  sound;  pj ,  the  undisturbed 
pressure.  Then  Jq.  1.3.1)  can  be  rewritten  as 
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where  t'p)j  is  the  overpressure  at  x  •  Define 

a  d l mens  ion  1  ess  variable, 


ajlx-l^t)  (‘.pi. 


Then 
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figure  3.1  exhibits  the  Taylor  velocity  or 
pressure  profile  in  a  plot  of  v/vj  or  (Ap)/(.'.p)j 
against  Z.  lhe  variable  Z  is  a  similarity  variable, 
since  the  velocity  or  pressure  profile  can  be 
obtained  as  a  unique  curve  against  Z  for  shock 
waves  with  different  strength  ( p ) j/p | .  and  it  will 
he  termed  the  distance  parameter. 


I'he  classical  Taylor  solution  iRef.  3)  for  weak, 
plane  shock  wave  transitions  is  expressed  In 
1  lghthill  (Ref.  .’ill  as 


Three  different  definitions  of  shock  thickness 
for  Z  are  also  shown  in  Tig.  3.1.  The  thickness 
l  \Z),(  is  defined  by 


w 


v:  (pij 

1  -’it  =  'Jv/j:’,=0  =  'di  '.p)/J:'_=0 

lhis  thickness  corresponds  to  the  velocity  or  den¬ 
sity-based  thickness,  and  it  has  been  used  in  some 
literature  for  shocks  of  moderate  strength.  The 
thicknesses  and  are  defined  by  the 

distances  for  the  overpressure  to  vary  from  10®  to 
90 o  and  from  5*  to  95‘1>,  respectively,  of  its  equi¬ 
librium  value  behind  the  shock.  The  last  definition 
was  used  by  I.ighthill  (Ref.  JO)  for  the  shock  thick¬ 
ness  derived  from  the  velocity  profile.  From  Hq. 
l3.oj  or  13."),  we  can  evaluate  the  values  of 
(■‘■3)o  and  l •.-)[)  as 


C.:)0  =  4  .66  \  i:z)^  =  5 . 1 J-* ,  UM)JJ  =  6.870 

These  will  be  termed  the  thickness  parameters.  The 
second  definition  of  the  shock  thickness  (10-90°s 
overpressure)  is  used  throughout  this  report  because 
it  can  give  a  reasonable  criterion  for  evaluating 
the  thickness  of  a  shock  wave  with  an  antisymmetric 
structure,  which  is  found  in  N-waves  and  in  partly 
or  fully  dispersed  plane  waves. 

The  actual  Taylor  thickness  (\x)A  and  the 
Taylor  rise  time  t£i)  (10-901)  overpressure)  can  be 
related  to  the  Taylor  thickness  parameter  (.\I)q  as 


'■■'V 


‘•-’o 


♦1  I '•!>), /p 


(3.8) 


from  !q.  1.5.51,  where  tr(,  is  the  Taylor  rise  time 
corresponding  to  the  Taylor  thickness  (,5x)(j.  We 
assume  tr()  =  ('.xlo/aj,  since  the  wave  speed  is 
nearly  equal  to  a]  for  very  weak  waves. 


be  realized  because  it  requires  an  infinitely  long 
time  for  the  wave  to  reach  a  steady  state  through 
viscous  diffusive  action.  In  the  case  when  the 
shock  strength  increases,  the  nonlinear  effects  are 
strengthened,  while  the  diffusive  effects  remain 
unchanged.  However,  the  shock  thickness  cannot  be 
less  than  the  molecular  mean-free-paths,  since  the 
shock  compression  process  is  after  all  a  result  of 
molecular  collisions.  In  other  words,  for  strong 
shocks,  the  shock  thickness  has  a  lower  limit  which 
is  controlled  by  molecular-collision  processes. 

Figure  3.3  shows  a  comparison  between  the  ex¬ 
perimental  and  theoretical  (Taylor)  rise  time  tr 
vs  the  maximum  overpressure  (Ap)max.  The  Taylor 
curves  shown  in  Fig.  3.2  are  reproduced  for  Tj  = 

273  K  and  290  K.  As  seen  from  Fig.  3.3,  the  rise 
times  for  the  spark  data  (Series  I  and  II)  are 
shorter  than  the  Taylor  rise  times  for  the  same 
maximum  overpressure,  while  the  rise  times  for  the 
exploding-wire  data  (Series  III  and  IV)  are  longer. 
Both  data  do  not  coincide  with  the  Taylor  curves. 

It  is  clearly  seen  that  the  Taylor  rise  times  for 
steady  viscous  shocks  can  give  no  reasonable  explan¬ 
ation  for  the  observed  rise  times  for  weak  spherical 
N-waves.  Therefore,  another  analysis  is  required 
for  this  purpose. 

3.2  Viscous  Plane  N-Waves 


In  this  section,  consideration  is  given  to  the 
case  of  a  kautnecd  S-ixtOe ,  which  is  produced  by 
moving  a  piston  forward  and  then  retracting  it  to 
its  original  position  in  a  tube.  The  generated 
plane  N-wave  gradually  decays  due  to  viscous  effects 
as  it  proceeds.  1. ighth ill  (Ref.  20)  solved  this 
problem  and  obtained  a  similar  solution  for  weak 
plane  N-waves,  where  the  velocity  profile  is  given 
as 


In  Fig.  3.2,  the  Taylor  thickness  (/>x)y  or  the 
Fay  lor  rise  time  try  are  plotted  in  a  nondimensional 
form  against  l.'.plj/pi  for  a  range  of  ( .'.p ) _/pi  = 

10’5  -  It'"-'  or  I'plj  =  1  Pa  -  100  Pa  in  the  atmos¬ 
phere.  At  NIT  for  air  =  1.333  x  10-5  m-/s, 

=  2/3,  ■  =  1 .4 ,  Pr  =  0.7  and,  from  F.q .  (5.2), 

■  -  5.43  x  1U'3  m-’/s.  Using  a )  =  331."  m/s,  the 
characteristic  length  and  time  are 

■/aj  =  1.03  x  10'  m,  '/a^  =  3.1  x  10"l°  sec 

Therefore,  for  ( ’.pjj/pj  =  I0~^  or  (J.p)j  =  10  Pa  at 
NTP,  Then  (\x)o  =  5.3  mm  and  tro  =  16  -sec,  from 
Fig.  3.J.  The  Taylor  thickness  or  rise  time  is 
inversely  proportional  to  the  shock  strength 
l  P)j/Pl-  As  the  shock  speed  is  weakened,  the 
Taylor  thickness  increases  and  tends  to  infinity 
as  (  '.p)j  -  0. 


u  - - ^ -  i •  V  ) 

1  ♦  exp  ( X“/J  *  t )  /  •.  exp  ( Re )  -  ll 

where  X  is  a  coordinate  measured  in  3  frame  of 
reference  which  moves  in  the  same  direction  as  the 
waves,  with  an  undisturbed  speed  of  sound  aj  and  is 
defined  as  X  =  x  -  a j t ;  u  is  the  excess  wavelet 
velocity  whose  variations  are  responsible  for  the 
\  *?  j*;v'  effects  and  is  defined  as  u  =  a  ♦  v  -  a] 

(a  is  the  local  speed  of  sound,  v,  the  particle 
velocity);  Re  is  a  Reynolds  number  of  each  half  of 
the  N-waves,  which  is  defined  in  terms  of  the  mass 
flow  in  that  half.  For  example,  for  the  front  half 

Re  =  i  J  udX  (5.10) 

X 

n 


As  mentioned  at  the  beginning  of  this  section, 
the  balance  between  the  finite-amplitude  (nonlinear) 
compression  effects  and  the  viscous-dif fusion 
effects  determines  the  thickness  of  a  steady  shock 
wave.  As  the  wave  is  weakened,  the  nonlinear 
effects  are  gradually  diminished,  while  the  viscous- 
diffusion  effects  remain  unchanged  regardless  of 
the  shock  strength.  Therefore,  for  very  weak  shocks, 
the  diffusion  effects  exceed  overwhelmingly  the 
nonlinear  compression  effects  and  broaden  the  shock 
thickness  to  very  large  values.  In  the  limit  of 
(Ap)>  "  0,  the  nonlinear  effects  disappear  and  only 
the  diffusion  effects  remain,  so  that  the  thickness 
tends  to  infinity.  However,  in  an  actual  case,  the 
steady  structure  of  such  a  very  weak  wave  would  not 


where  Xn  is  the  node  u  =  0  and  *  is  the  diffusivity 
defined  by  Kq.  (3.2).  Note  that  Re  is  not  invariant, 
but  varies  with  time  as  the  mass  flow  varies  with 
the  decay  of  the  wave.  The  riLi'isr-j  .V-uaiv  means 
that  its  total  mass  flow  always  vanishes  as 

|  udX  =  0 

From  the  nonlinear  wave  relation, 

U  *  v  (3.11) 
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Using  hqs.  (.3.3)  and  (5.11)  and  defining  the 
similarity  variables 


r  =  ^ia,  i  m, 

-■  1  '  •  p, 


(  3 . 1  ) 


then  from  Fq.  (3.9), 


P  .  :  I  1  .  I 

exp (Re)  -  1 


1-1 


(3.131 


figure  3.4  shows  the  pressure  profiles  for  several 
different  Reynolds  number  Re  in  a  plot  of  P  against 


For  a  given  Reynolds  number  Re,  we  can  obtain 
I’max  Uhe  maximum  value  of  P) ,  3;  (the  shock  thick¬ 
ness  defined  by  1 0 - 9 0 “a  overpressure)  and  "j  (the 
half  length  of  the  N-wave  measured  from  the  origin 
to  the  point  of  10°a  overpressure  in  the  wave  front) 
Then  the  following  parameters  can  be  obtained: 


AZ  =  ii_r  (“‘"OP 

>♦1  max 


VAX)  (Ap)n 


(3.14) 


a.X  (ip) 

1  d  r  max 


“d  "  -,♦!  'd  max 


“f'd  '  ‘P’max 


(3.13) 


where  AX  is  the  shock  thickness  corresponding  to 
i\r, ,  AX  =  Aff't;  Xg,  the  half  length  of  the  N-wave 
corresponding  to  Ad.  Xj  =  AdV^t;  (Ap)max>  the  maxi¬ 
mum  value  of  (Ap) .  The  parameters  AZ  and  ij  corre¬ 
spond  to  the  shock  thickness  and  the  flow  duration 
of  the  N-wave  with  reference  to  the  dimensionless 
variable  Z,  which  is  defined  similarly  to  Eq.  (3.5) 
as 


aIX  W-ax 
pl 


(3.1b) 


AZ  is  the  thickness  parameter  defined  in  the  previ¬ 
ous  section  and  Zj  will  be  termed  the  duration 
parameter.  Details  of  the  derivation  of  P,1K1X,  AZ 
and  Zd  are  given  in  Appendix  A. 

Figure  3.5  exhibits  the  pressure  profiles  for 
the  same  cases  as  shown  in  Fig.  3.4  in  a  plot  of 
(Ap)/(Ap}max  against  Z  -  Zg,  where  Zo  is  the  Z 
at  (Ap)/(Ap)ma/  =  0.5.  The  solid  line  indicates 
the  Taylor  solution  for  steady  plane  waves,  which 
is  given  by  Eq.  (3.6)  or  (3.7).  The  l.ighthill 
N-wave  solution  approaches  the  Taylor  solution  as 
Zd  -*•  00  or  Re  -»  This  can  also  be  shown  from  F.q. 

^3 . 13)  as  follows.  Assume  that  P  reaches  its  maximum 
pmax  at  ’  =  Am  f°r  large  Re.  Then,  approximately. 


same  form  as  Fq.  (3.b),  the  Taylor  solution,  since 
Pmax\'  can  be  replaced  by  Z  -  Zg,  where  Zg  is  the 
Z  at  f  =  r,m.  It  should  be  noted  that  the  shock 
thickness  decreases  as  the  Reynolds  number  Re  or 
the  duration  parameter  Zd  decreases  for  the  same 
maximum  overpressure. 

In  Fig.  3.6,  the  ratio  of  the  thickness  para¬ 
meter  (AZ)/ (.'  Z)0  is  plotted  against  the  duration 
parameter  Zd,  where  (  ,Z)g  is  the  (AZ)  for  Zd  *  • 

( lay  lor  solution)  and  is  given  by  (AZJg  -  S.127. 
this  figure  clearly  shows  the  dependence  of  the 
shock  thickness  on  the  duration  of  the  N-wave.  As 
the  duration  or  the  maximum  overpressure  increases, 
the  shock  thickness  approaches  the  Taylor  value.  As 
the  duration  or  the  maximum  overpressure  increases, 
the  shock  thickness  approaches  the  Taylor  value.  As 
the  duration  or  the  maximum  overpressure  decreases, 
the  deviation  from  the  Taylor  value  increases. 

In  Fig.  3.7,  the  normalized  shock  thickness 
(AX ) / ( ' /a i )  or  the  normalized  rise  time  tr/[/a.-l 
is  plotted  against  the  normalized  maximum  over¬ 
pressure  ( Ap ) max ,  p  j  for  the  normalized  duration 
Xd/t  7a,  i  or  tg/(7a,-l  -  constant.  It  can  also 
be  seen  from  Fig.  3.'  that  the  shock  thickness 
or  rise  time  decreases  for  a  fixed  maximum  over¬ 
pressure  (Aplmax  as  the  duration  of  N-wave  de¬ 
creases.  This  is  the  .V- •  ‘‘f-  ‘ *  described  in 
the  Introduction. 

In  Fig.  3.8,  the  experimental  data  of  Ref.  6 
are  compared  with  the  l.ighthill  solutions  for 
N-waves.  The  rise  time  tr  is  plotted  against  the 
maximum  overpressure  (Ap)max.  The  solid  lines 
exhibit  the  N-wave  solutions  for  tj  «  50  ..sec  and 
'0  „sec  which  correspond  to  the  half-durations  in 
the  spark  experiments.  The  Iaylor  rise  time  for 
T i  =  2"3  k  is  also  plotted  against  (Ap)max.  The 
figure  shows  that  the  rise  times  obtained  in  the 
spark  experiments  are  adequately  explained  by  the 
Lighthill  model  of  viscous  (frozen)  N-wave  shocks 
though  the  measured  rise  times  slightly  deviate 
from  the  theoretical  curves  in  the  range  of  the 
lower  overpressure. 

In  Fig.  3.9,  the  experimental  data  are  plotted 
on  a  figure  showing  the  ratio  of  the  thickness 
parameters  (AZ)/(AZ)g  vs  the  duration  parameter  Zg, 
shown  in  Fig.  3,b.  The  data  cover  the  range  of 
Zg  =  10-100,  in  which  the  spark  data  lie  between 
Zd  =  10  and  bO  and  the  exploding-wire  data  lie 
between  Zd  =  50  and  100.  Using  the  duration  para¬ 
meter  Zd,  the  data  may  be  categorized  into  three 
domains.  Above  Zg  '  50,  the  measured  (AZ) -values 
deviate  from  the  Lighthill  curve  and  steeply  in¬ 
crease  with  increasing  Zd .  In  the  range  Zd  =  15-50, 
the  measured  (AZ)-values  nearly  coincide  with  the 
Lighthill  curve,  a  scatter  of  the  data  exists. 

Below  Zg  15,  the  measured  (AZ)-values  again 
deviate  from  the  curve  and  steeply  decrease  with 
decreasing  Zg .  The  broken  lines  are  drawn  to 
stress  the  tendency  of  the  data. 


P  A  ,  Re  =  A  2/2 

max  m  m 


Put  A  =  Am  ♦  A'  (A*  •<  Am) •  then 


P  -  A  [1  ♦  exp (A  A'))'1  ^  P  |1  ♦  exp(P  A’))'1 
m1  r  m  max 1  ‘  max  ' 

(3.P) 

in  the  limit  of  Re  *  ■».  Equation  (3.17)  has  the 


Figure  3.10  shows  a  comparison  between  the 
observed  and  Lighthill  N-wave  pressure  profiles. 
Typical  profiles  in  the  Series  I -IV  are  plotted  by 
the  broken  lines  in  comparison  with  the  correspond¬ 
ing  analytical  ones,  which  are  evaluated  from  Fq . 
(3.13)  to  have  the  same  maximum  overpressure 
(Ap),,,^  ant*  the  same  hal  f -durat  ion  tg  as  the  experi¬ 
mental  ones,  and  plotted  by  the  solid  lines  to  fit 
each  other  at  the  nodes  of  the  N-waves.  As  seen 
from  the  figure,  the  pressure  profiles  observed  in 
the  spark  experiments  (Series  1  and  II;  Figs. 

3.10(a)  and  (b)]  nearly  coincide  with  the  analytical 
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ones,  while  the  pressure  profiles  observed  in  the 
exploding-wire  experiments  [Series  111  and  IV;  Fig. 
3.10(c)  and  (d)]  deviate  from  those  predicted  ana¬ 
lytically.  The  main  difference  between  both 
experiments  is  that  of  the  half-duration  of  the 
N-wave.  Figure  3.10,  as  well  as  Figs.  3.8  and  3.9, 
suggests  that  the  lighthiil  viscous  N-wave  model 
does  not  always  explain  the  rise  tunes  of  N-waves 
over  the  entire  range  of  tj  or  Zj. 

3 . 3  Nonstationary  Viscous  Plane  Waves 


In  this  section,  consideration  is  given  to  a 
nonstationary  plane  wave,  which  is  generated  by 
the  impulsive  motion  of  a  piston  in  a  tube.  The 
initially  discontinuous  wave-front  is  smoothed  out 
due  to  viscous  diffusion  and  it  tends  to  form  a 
final  steady  profile.  It  will  be  shown  in  the 
succeeding  sections  that  this  process  of  shock 
thickening  effcsr)  plays  an  important 

role  in  determining  the  rise  times  of  weak  spherical 
N-wa ves . 


LdLL  = 
i  ‘p) , 


exp 


(3.22) 


which  is  the  Taylor  solution  for  steady  plane  waves, 
F.q .  i  3  .  h  |  . 


figure  3.11  shows  the  pressure  profiles  for 
several  different  time  parameters  in  a  plot  of 
('p)/ip)j  against  the  distance  parameter  Z.  The 
pressure  profile  approaches  the  Taylor  profile  as 
-  It  can  he  seen  that  the  shock  thickness 
i  Z)  increases  as  increases  [whether  based  on 
maximum  slope  or  10-90&  of  <  ’ p) / (  \p) 2I • 

In  Fig.  3.12,  the  ratio  of  the  thickness  para¬ 
meters  |'Zj/|’Z)o  is  plotted  against  »T.  If  we 
define  a  character  1 st ic -t ime  parameter  of  shock 
thickening  s  as  at  (  Z)/(‘.Z)q  =  0.99,  then 


>  =3.3 

s 


or 


s 
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Lighthiil  (Ref.  20)  has  given  a  solution  for 
the  nonstationary  plane  wave  by  solving  Burger's 
Equation.  He  obtained  the  following  result: 

u-, 

U  ( X ,  t )  =  - = - : - 


-v~/2  *t 


u,  X  -  4u,t; 


dv 


1  exp 


-X 


-v  /2‘t  . 
e  dy 


X-u,t 


(3.18) 


in  which  the  initial  wave  form  is  given  by 

u(X,  0)  =  u,  for  X  •  0,  and  :ero  for  X  0  ('•l'-1) 


from  which  the  corresponding  time  ts  and  distance 
xs  are  obtained  from  F.q.  (3.21)  as 

t  x  -I  \p  |,  -2 

s  _  s  _  .  '2  ,,  . 


which  are  designated  as  the  shock-thickening  time 
and  distance,  respectively.  These  are  inversely 
proportional  to  the  square  of  the  shock  strength 
(  p ) 2/ PI •  This  means  that  it  takes  a  progressively 
longer  time  and  distance  to  reach  a  final  steady- 
state  for  weaker  shock  waves  or  for  lower  lip)’/p]. 
Physically,  this  tendency  of  longer  shock-thicken¬ 
ing  time  or  distance  for  weaker  shocks  is  attributed 
to  the  decline  of  shock  steepening  due  to  nonlinear 
(convective)  effects. 


where  u2  is  the  excess  wavelet  velocity  for  X 
Using  F.qs.  (3.3)  and  (3.11), 


erf 

!  ,♦!  1  C 


(,'.p)  ,  ,*1  „ 

WTJ  ■  1  *  eXP  :  7~  “ 


:  .-1  /:  -1 

,-tt  ‘  4.  «/  2 


erf 


"  4/  1/2  ,  j 


(3.20) 


where  Z  and  :  are  the  distance  parameter  and  the 
time  parameter,  respectively,  defined  by 


% 

Pi  ’ 


j  (  '>p)2 


(3.21) 

The  complementary  error  function  is  defined  by 

2 

erfc(X)  =  |  e  J  dy 
X 

Note  that  the  shock  strength  ("plr/pi  depends  on 
the  piston  velocity  v2  [-  2u2/( > ♦!)] and  is  invari¬ 
ant  throughout  the  process.  When  t  •»•»,  f.q.  (3.20) 
becomes 


In  Fig.  3.13,  the  normal  iced  shock-thickening 
t  ime  ts/(  '/a) ‘ )  or  distance  x s/ (  6/ a  1 )  is  plotted 
against  the  shock  strength  (.'p)2/pi-  The  time 
scale  on  the  right  hand  side  indicates  the  shock 
thickening  time  at  NIP  in  air.  For  ( ' p ) 2/p 1  =  10"^ 
or  (\p)2  =  10  Pa,  ts  =  1  sec  or  xs  =  330m.  These 
values  suggest  that  the  nonstationary  effect  on  the 
rise  time  or  the  shock  thickness  becomes  very 
important  for  weak  shock  waves,  for  it  takes  a  long 
time  or  a  large  distance  to  reach  a  steady  state. 

This  result  is  of  value  in  interpreting  Fig.  3.4 
or  3.5,  which  provides  solutions  for  quasi-station- 
ary  N-waves  at  the  final  values  after  a  very  long 
time  without  specifying  how  long  it  may  actually 
take.  The  above  solution  quantifies  the  time  or 
distance  in  specific  cases.  The  spark  and  exploding- 
wire  generated  N-waves,  described  in  Section  2,  are 
also  expected  to  be  affected  by  this  nonstationary 
effect,  since  the  maximum  overpressures  are  below 
20  Pa  only  over  a  distance  of  10m. 

3.4  Shock  Transitions  with  Vibrational  F.xcitation 


The  structure  and  thickness  of  shock  waves 
with  vibrational  excitation  in  air  will  be  consid¬ 
ered  now.  The  analytical  results  of  Polyakova, 
Solyan  and  Khokhlov  (Ref.  21)  and  Johannsen  and 
Hodgson  (Ref.  12)  for  plane  dispersed  waves  are 
re-examined  and  compared  with  Mol st-Jensen ' s  data 
(Ref.  b ) .  Furthermore,  extensions  of  Lighthiil 
solutions  for  N-waves  and  nonstationary  waves  to 
shock  transitions  with  vibrational  excitation  are 


i 


I 


made  possible  by  using  a  bulk-viscosity  concept. 

For  weak  shock  waves  with  vibrational  excita¬ 
tion,  steady  shock  waves  are  formed  as  a  result  of 
a  balance  between  the  wave-form-steepening  tendency 
due  to  finite-amplitude-compression  effects  and 
the  wave-easing  tendency  due  to  both  effects  of 
viscous  diffusion  and  vibrational  relaxation.  For 
very  weak  waves,  the  compression  effects  diminish 
and  the  wave-form-easing  effects  become  predominant. 
As  discussed  in  Section  3.1  for  vt'aooua  or  frjzcn 
shock  transitions,  in  the  limit  of  (.'pb'  0,  the 
nonlinear  compression  effects  disappear  and  the 
wave-form-easing  effects  remain,  so  that  the  wave 
thickness  tends  to  infinity.  For  weak  shocks 
whose  strengths  are  slightly  above  the  limit  of 
zero  overpressure,  the  vibrational  relaxation  is 
more  effective  than  the  viscous  diffusion  for  the 
wave-easing  tendency.  In  this  case,  the  compres¬ 
sion  process  is  so  slow  that  the  energy  dissipation 
due  to  vibrational  nonequilibrium  becomes  predomin¬ 
ant  compared  with  that  due  to  translational  and 
rotational  nonequilibrium  which  requires  a  more 
rapid  change  of  the  flow  properties.  As  the  wave 
strength  increases,  the  shock  thickness  decreases 
owing  to  the  increase  in  nonlinear-compression 
effects.  When  the  nonlinear-compression  effects 
overcome  the  wave-easing  effects  due  to  vibrational 
relaxation,  the  frozen  shock  transition  appears  in 
the  compression  process  of  the  wave. 

Figure  3.14  illustrates  these  two  types  of 
shock  transition  with  vibrational  excitation 
through  pressure  and  temperature  profiles.  The 
vibrational  temperature  Ty  is  also  plotted  to  show 
the  process  of  vibrational  energy  excitation.  The 
former  wave  dominated  by  the  vibrational  excitation 
is  called  a  fully  dispersed  wave,  and  the  latter 
wave  including  the  fpczcn  (relatively  sharp,  vis¬ 
cous)  shock  transition  is  called  a  partly  dispersed 
wave.  For  strong  shocks,  the  nonlinear  compression 
mainly  balances  with  the  viscous  diffusion,  though 
it  is  accompanied  by  the  slower  process  of  vibra¬ 
tional  excitation.  As  shown  in  Fig.  3.14,  for 
stronger  shocks ,  the  temperature  goes  up  to  the 
maximum  (Rankinc-Hugon iot )  value  through  the  a.:n 
shock  compression  and  then  it  falls  to  the  final 
equilibrium  state  through  the  tv/ajwf  n  z  re  as 
vibration  attains  its  share  of  energy. 


Polyakova  et  al  (Ref.  II)  have  obtained  an 
analytical  solution  for  the  structure  of  steady, 
plane  dispersed  waves  for  nonviscous  and  noncon¬ 
duct  ive  gases  as 


,  .k-1  2 

y  *  y0  <v0  *  vj  vQ 

- ; -  =  «.n  - 


J 


(V0  -  V) 


k-1 


(3.24) 


where  y  =  t  -  C/ae;  r.  =  l.agrangian  coordinate, 
ae  «  equilibrium  speed  of  sound;  yp  =  constant  of 
integration;  .j  =  vibrational  relaxation  time  for 
j-molecule;  v  =  velocity  in  a  moving  coordinate 
system,  vp  =  absolute  value  of  the  velocity  at  the 
spatial  coordinate  '  *  •  <»;  k  =  mae/(2vpt.); 
m  =  (af^  -  ae’)/ae^;  af  =  frozen  speed  of  sound; 

=-  =  (7  +  1 )  /2 . 


In  order  to  rewrite  Fq .  (3.24)  using  the 
normalized  overpressure  ( Ap )/ ("p ) i  and  the  distance 
parameter  Z,  which  were  introduced  in  the  previous 
sections,  introduce  two  quantities;  the  bulk  vis¬ 
cosity  and  a  critical  overpressure. 


The  bulk  viscosity  ( u v k j  for  the  .J-molecule 
can  be  expressed  as 


(u  )  .  = 

V  J 


;o(af 


-  a  “)  = 


(3.25) 


for  processes  sufficiently  slow,  where  .  q  is  the 
equilibrium  density  of  the  medium.  Then  the  dif- 
fusivity  ( •  v ) j  for  i-molecule  with  a  bulk  viscosity 
(Uv)j  can  be  expressed  as 


‘Vj  ■ (- 


j/iJ0 


.ma 
J  e 


(3.2b) 


This  diffusivity  will  be  used  as  a  reference  physic¬ 
al  property.  It  should  be  noted  that  the  use  of 
this  property  does  not  mean  that  the  vibrational 
relaxation  processes  can  always  be  replaced  by  the 
bulk  viscosity,  which  is  valid  only  for  processes 
sufficiently  slow. 

The  critical  overpressure  is  defined  as  the 
equilibrium  overpressure  behind  a  plane  dispersed 
wave  whose  wave  velocity  is  equal  to  the  frozen 
speed  of  sound.  When  the  equilibrium  overpressure 
exceeds  the  critical  overpressure,  the  steady  plane 
wave  is  a  partly  dispersed  wave  with  a  frozen 
(viscous)  shock  front,  which  is  followed  by  the 
vibrational  relaxation  region.  When  the  equilibrium 
overpressure  is  below  the  critical  overpressure,  the 
steady  plane  wave  is  a  fully  dispersed  wave  with  a 
smooth  transition,  which  is  controlled  by  the  vib¬ 
rational  excitation  of  the  molecules. 


The  equilibrium  overpressure  across  a  normal 
shock  wave  with  vibrational  excitation  can  be  given 


(Ap),  2,(Mf2-l)  ♦  2C,-l)(Mf2,-l)c. 

~p^~~  =  (7-1)  *"2(,-T)c  " 


(3.27) 


where  Mf  is  the  frozen  Mach  number,  c;  the  vibra¬ 
tional  specific  heat  for  j-molecule  normalized  by 
the  gas  constant  Cj  =  cj/R  in  which  cj  is  assumed 
to  be  constant  across  the  shock  wave.  Details  of 
the  derivation  of  Fq.  (3  .27)  are  given  in  Appendix 
B.  If  the  harmonic  oscillator  approximation  is 
applied  to  the  vibrational  energy  level,  the  vib¬ 
rational  specific  heats  for  Oj  and  in  air  may  be 
written  as 

”0  j 

c„  =  0.209  r  I  exp  -  —  J  (3.28a) 

.»  -  2  ,  ■-)  >, 

c  =  0.781  =—  i  exp  I  -  (3.28b) 

N  1-  T1  -  '  T1  1 


where  I j  is  the  initial  gas  temperature  (room  tem¬ 
perature),  -j  the  vibrational  characteristic  tem¬ 
perature  "q  '=  2239.1  K,  ej.  =  3352  1C.  For  Mf  =  1  , 
we  have  the  critical  overpressure  for  the  j-molecule 
as 


(Ap).,  ■ 
*  cr , j 


2(i-l)2c 
(-*11  -  2(,-l)c. 


>♦1  j 


(3.29) 


for  Cj  1,  which  is  usually  valid  for  atmospheric 
air,  as  very  little  vibrational  excitation  can  exist 
at  nearly  room  temperature.  The  critical  overpres- 


I 


r 


sure  (Apjcr  ,  depends  on  the  gas  temperature  Ti, 
since  the  vibrational  specific  heat  Cj  depends  on 
I'l  ■ 


In  Fig.  3.15,  the  critical  overpressures 
('P)cr.O  and  (-'p)cr,0*N  are  plotted  against  Tj. 
The  lines  denoted  by  O;  and  O2+N2  are  calculated 
from 


t;'p)cr,0  .  2(y-l)2 
Pj  ’  <+l  C° 

( tp)er,0»N  2(,-l)2 

p,  , +1  tc0  CN 


(3.30a) 


(3.30b) 


respectively.  That  is,  in  the  former  case,  only 
the  vibrational  excitation  for  O-molecules  in  air 
is  taken  into  account.  For  (ip) 2  S  (ip)cr  j,  the 
steady  plane  wave  is  fully  dispersed,  and  tor 
(.\p)2  >  (ip)cr  j  it  is  partly  dispersed. 

The  diffusivity  (5y)j  can  be  expressed  by  the 
critical  overpressure  as 


Vj  = 


1*1 


I . 

J 


WCr, 


for  c .  <c 
} 


(3.31) 


The  parameter  k,  which  appears  in  Eq.  (3.24),  can 
be  rewritten  as 


1  (-'"p)2 
F  t'Plcr.J 


for  c  .  << 
J 


(3.32) 


tion  of  Eq.  (3.34)  are  also  given  in  Appendix  B. 


Johannesen  and  Hodgson  (Ref.  12)  have  also 
obtained  an  exact  solution  for  steady  plane  dis¬ 
persed  waves  for  nonviscous  and  non-conduct ive 
gases,  as  follows: 


Mf  [0  +  1 )  ♦  2Q-l)c.] 
—  - 


-(,+l)Mf2  — 


U  7 

1 - - 


f  I  ,  U 

: -  x.n  1  -  — 


u 
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|  >  7  ^2 

|1  +  -  (>  +  l)M/  ~  I  — 

1.  1  1  U1 


2 

U1 


1  -  Q2/Ul 


U2  1 

-  —  !  (3.37) 


where  u  is  the  flow  velocity,  uj ,  Uj  are  the  flow 
velocities  at  x  ■*  *  .  Using  the  relations 


1 


u 


u  U2  1  (2p)2  j-  (ip)  I 

Uj  Uj  -,Mf2  pl  -  (ip)2  J 


(3.38a) 


(3.38b) 


and  neglecting  the  higher  order  terms  of  0(Cj),  the 
same  equation  as  Eq.  (3.34)  is  obtained,  which  was 
derived  from  the  Polyakova  et  al  (Ref.  21)  formula, 
by  using  the  distance  parameter  defined  by 


That  is,  the  parameter  k  is  the  ratio  of  the  critic¬ 
al  and  equilibrium  overpressures.  For  k  1,  the 
wave  is  a  partly  dispersed  wave,  and  for  k  1  the 
wave  is  a  fully  dispersed  wave.  The  derivations  of 
Eqs .  (3.31)  and  (3.32)  are  given  in  Appendix  B. 

Using  the  relation 


V  (ip)2  2y  x  (ip)2 

(A  )  .  p,  ’  ,  +1  a,  t .  (Ap) 

v  j  1  j  1  cr,j 


(3.39) 


Further  details  can  be  found  in  Appendix  B.  Equa¬ 
tion  (3.34)  will  be  used  as  a  solution  for  steady 
plane  dispersed  waves. 


1  +  2-  -  2  jdE) 

v  ('P)2 
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(3.33) 


then 
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(3.34) 


from  F.q.  (3.24),  where  the  distance  parameter  Z  is 
defined  as 


In  the  limit  of  a  weak  wave  ( ip ) 2  "*  0,  bq. 
(3.34)  tends  to 


>*1 

2  > 


m  .  -  m. 

L  (Ap) - 


- 


!  m. 

L  (AP); 


(3.40) 


This  has  the  same  form  as  the  Taylor  solution,  Eq. 
(3.7),  in  which  the  diffusivity  f  is  replaced  by 
(iy):.  In  the  limit  of  weak  shocks,  the  shock 
compression  process  is  infinitely  slow,  so  that 
the  bulk  viscosity  concept  may  be  applied  to  the 
vibrational  relaxation  process.  The  solution,  in 
which  the  diffusivity  6  is  replaced  by  (\0j  or 
+  +  ( <5 v) j ,  will  be  called  the  modified  Taylor 
solution. 


aj2y  ( Ap) 2 


(3.35) 


in  a  similar  way  to  Eqs.  (3.16)  and  (3.21)  in  the 
previous  sections,  it  can  be  rewritten  as 


Z 


2y _ 


(Ap). 


(Ap) 


cr,j 


(3.36) 


Zq  is  an  arbitrary  constant.  Details  of  the  deriva- 


Figure  3.16  shows  the  pressure  profiles  for 
several  different  values  of  (Ap) 2/ (Ap)cr, j  in  a 
plot  of  ( Ap ) / ( Ap) 2  against  Z-Zq.  The  curve  for 
(Ap)2/(Ap)cr  ,  -»  0  corresponds  to  the  modified 
Taylor  solution.  For  partly-dispersed  waves 
( (Ap ) 2  •  (Ap)cr  j],  there  appears  a  discontinuous 
shock  front.  Die  overpressure  (Ap) f  immediately 
behind  the  frozen  shock  is  given  by 


(Ap)f 

WT 


(Ap) 


(Ap)  - 


cr  <} 


(3.41) 
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In  Fig.  5.1b,  the  chain  curve  indicates  the  pressure 
profile  for  lAp)>/( Ap)cr  j  =  2,  in  which  the  discon¬ 
tinuous  shock  strength  at' I  =  Zo  is  (.\p)f  =  0 . S ( Ap )  > - 

The  thickness  parameter  (AZ)  is  defined  by  the 
10-901,  equilibrium  overpressure,  and  can  be  related 
to  the  rise  time  tr  as 


*  r 


13.42) 


For  fully-dispersed  waves  where  [('41)2  •  l  p)cr  jl> 
then  from  liq.  (3.34) 


•  n  9 


5. 12" 


l-'o 


13.43.) 


regardless  of  the  value  of  l'p ) j/ (.'.p)cr  j-  That  is, 
the  thickness  or  rise  time  of  a  fully-dispersed 
wave,  which  is  based  on  the  10-90“,.  equilibrium 
overpressure,  has  the  same  value  of  the  thickness 
parameter  as  the  Taylor  thickness  or  rise  time,  if 
the  diffusivity  ( ' v ) j  is  used  instead  of  '. 

In  Fig.  3.1',  the  ratio  of  the  thickness  para¬ 
meter  (  Z ) / ( ' Z ) 0  is  plotted  against  the  equilibrium 
overpressure  normaliced  by  the  critical  overpressure 
for  fully  and  partly-dispersed  waves.  It  can  be 
seen  in  the  figure  that  the  effect  of  dispersion  on 
l.'.Z)  remains  up  to  iplj  =  10(.'.p)cr  j .  This  means 
that  the  rise  times  for  steady  plane  waves  arc 
affected  by  the  vibrational  relaxation  up  to  l.'p)j  = 
500-1,000  Pa  in  air,  since  ( 'p)cr  j  =  30-100  Pa  in 
the  usual  range  of  ground  temperatures  (see  Fig. 
3.15). 

The  Lighthill  solutions  for  N -waves  (Section 
3.2)  and  nonstationary  waves  (Section  3.3/  mac  he 
applied  to  fully-dispersed  waves  for  small  ('P)v 
f p)cr  ;  by  replacing  the  diffusivity  •  with  the 
vibrational  diffusivity  (v)j  in  order  to  provide  a 
rough  estimate  of  the  .'.-'.’.-'.'and  r.  j-  ' 
effects  on  the  thickness  or  rise  time  of  dispersed 
waves  with  vibrational  excitation. 


13.23).  For  t  .p)  2/Pl  *  1°~4  or  l  P  J  2  =  1(1  *’a  • 
ts  =  15  sec  or  xs  =  5  km.  This  shows  that  it  is 
very  difficult  to  obtain  a  plane  dispersed  wave 
in  a  steady  state  on  a  laboratory  scale. 

As  for  the  X-wave  effect,  the  dispersed  wave 
is  affected  by  the  expansion  behind  the  shock 
front  more  seriously  than  the  viscous  wave,  since 
the  former  has  a  larger  thickness  than  the  latter 
for  the  same  duration  and  maximum  overpressure. 
Therefore,  both  the  .V-u Ktve  and  nonj ta t ionaru 
effects  will  seriously  modify  N-waves  with  vibra¬ 
tional  nonequilibrium. 

In  Fig.  3.18,  the  exp] oding-wi re  data  are 
compared  with  several  theoretical  curves  in  a  plot 
of  the  rise  time  tr  against  the  maximum  overpressure 
(’P)max*  The  chain  lines  indicate  the  Taylor  and 
the  modified  Taylor  rise  times.  The  broken  lines 
indicate  the  modified  I. i ph thill  rise  times  for 
N-waves  of  tj  =  100  and  1  JO  ..sec.  The  vibrational 
diffusivity  |  * e ! i)  for  oxygen  is  used  for  the  modi¬ 
fied  la'  1  >r  and  the  modified  lighthill  solutions. 

All  cm  '  are  evaluated  for  the  gas  temperature 
II  =  J80  k  and  the  relative  humidity  RH  =  8  .  5  ^5. 
iSeries  IV).  The  corresponding  vibrat ional -relaxa¬ 
tion  time  and  the  critical  overpressure  for  oxygen 
are  about  3.  "3  ..sec  and  bl  Pa,  respectively.  The 
measured  rise  times  are  much  shorter  than  the 
modified  lighthill  rise  times  for  fully-dispersed 
N-waves.  This  discrepancy  can  be  attributed  to 
the  \  •  ’  •*»;•/!*  .  effect. 


figure  5.19  shows  a  comparison  between  the 
observed  and  modified  lighthill  N'-wave  pressure 
profiles  in  a  similar  way  to  fig.  3.10  for  viscous 
N-waves.  Typical  profiles  from  Series  I  - 1 V  are 
plotted  using  broken  lines  in  comparison  with  the 
corresponding  analytical  ones  shown  as  solid  lines, 
which  are  evaluated  from  Fq.  (5.  15),  with  •  replaced 
by  i ’ v l p •  Ihe  profiles  have  the  same  maximum  over¬ 
pressure  i  'p)max  and  the  same  ha  1 f -duruc ion  tj  as 
the  experimental  ones,  and  fit  at  the  nodes  of  the 
\  waves.  By  contrast  to  fig.  5. lit,  the  discrepancy 
between  the  observed  anJ  analytical  profiles  is 
c l ea  r . 


Assume  that, 

a.  -  531."  m/s,  =  10  ^  sec, 

1  J 

(ApJ  .  =  50  Pa,  p,  =  101.3  KPa 
r  cr,.i  1  1 

then,  from  f.q.  (3.31), 

(  •  )  .  =  4'  x  llT3  m'/ s 
'  J 

(Compare  with  *  above  for  translation  and  rotation 
of  3.43  x  10~*>  m^/s,  that  is.  the  dispersed  shock 
structure  is  entirely  controlled  by  the  vibrational 
relaxat ion . ) 

-?  2  -in 

(‘v)./aj  =  14  x  10  m,  =  43  x  10  luSec 

(Compare  with  fl/aj  =  1.03  x  I0-7m  and  =  5.1  x 

1 0 " 1 O  sec  noted  above.)  These  values  arc  about  ten 
times  as  long  as  the  ones  evaluated  for  viscous 
shocks  in  Section  3.1.  This  means  that  the  thick¬ 
ness  or  the  rise  time  of  a  plane  dispersed  wave  is 
about  ten  times  as  long  as  that  of  a  viscous  shock 
wave  for  the  same  shock  strength  ( ' p ) j/p i •  The 
shock -thickening  time  or  distance  of  an  impulsive 
step  wave  is  also  tenfold  greater  for  a  dispersed 
wave  than  for  a  viscous  wave,  as  seen  from  Fq. 


lo  conclude  this  section,  consideration  is 
given  to  .»  characteristic  feature  of  weak  N-waves 
with  vibrational  nonequilibrium.  figure  5 . JO  illus¬ 
trates  a  classification  of  weak  N-waves  bv  their 
degree  of  vibrational  nonequilibrium.  The  profiles 
of  gas  and  vibrational  temperatures  are  plotted 
under  the  following  assumptions:  (if  the  maximum 
(peak l  overpressures  are  below  the  critical  over¬ 
pressure  for  steady,  plane  waves;  (ii)  the  maximum 
overpressure  is  the  same  for  all  cases  in  Fig.  3. JO; 
liii)  only  one  mode  of  vibrational  excitation  is 
considered.  As  seen,  the  N-waves  can  be  classified 
into  five  categories:  i  a  )  quasi -equi  1  ibnum  wave, 

(b)  moderate  1 y-nonequ 1 1 ib rium  wave,  ( c )  highly -non 
equilibrium  wave,  id)  nearly -frozen  wave,  (e)  quasi- 
frozen  wave. 

The  degree  of  excitation  of  vibrational  energy 
is  denoted  by  the  vibrational  temperature  Tv,  which 
is  plotted  by  broken  lines  in  fig.  5. JO.  The  time 
lag  between  the  gas  and  vibrational  temperatures 
corresponds  to  the  vibrational  relaxation  time 
In  a  quasi -equ l l ibr mm  wave,  the  vibrational  temper 
at ure  nearly  follows  the  gas  temperature.  This  is 
the  case  where  the  concept  of  bulk  viscosity  is 
valid  and  thr  modified  lighthill  solution  for  N-waves 
may  be  applied.  The  structure  of  the  shock  front  is 
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I 


controlled  by  the  vibrational  relaxation,  that  is, 
the  wave  is  a  ful lv -dispersed  wave.  In  a  moder- 
at ely -nonequi 1 ibrium  wave,  an  appreciable  deviation 
of  the  vibrational  temperature  from  the  gas  tempera¬ 
ture  can  be  seen.  In  this  case,  the  concept  of  bulk 
viscosity  cannot  be  applied  to  the  vibrational  re¬ 
laxation,  though  the  front  structure  is  still 
controlled  by  the  vibrational  relaxation.  This  wave 
can  also  be  considered  as  a  fully-dispersed  wave. 

In  a  highly -nonequi librium  wave,  the  front  structure 
is  controlled  by  both  processes  of  vibrational  ex¬ 
citation  and  viscous  dissipation.  The  wave  becomes 
a  partly-dispersed  wave  in  the  sense  that  the  front 
structure  is  partly  controlled  by  viscous  effect. 

The  structure  of  a  nearly-frozen- flow  frozen  wave 
is  mainly  controlled  by  viscous  effect,  though 
v Lb rational  excitation  still  remains  in  the  rest  of 
the  flow  field.  In  a  quasi -frozen  wave,  the  vib¬ 
rational  excitation  is  marginal  so  that  the  whole 
flow  field  can  be  considered  as  frozen. 

The  discrepancy  between  the  observed  and  ana- 
1  >  1 1  c  a  1  rise  times  and  pressure  profiles  described 
in  the  preceding  sections  may  be  explained  by 
considering  the  above  classification  for  \-waves. 

The  N-waves  generated  by  sparks  could  be  highlv- 
nonequilibrium  waves  or  nearly-frozen  waves,  since 
the  front  structures  seem  to  be  mainly  controlled 
by  viscous  effect.  The  N-waves  generated  by  explod¬ 
ing  wires  could  be  moderately-nonequ i 1 ibrium  waves. 
The  coupling  of  the  N-wavo  and  nonstationary  effects 
would  make  the  situations  even  more  complex. 


4.  RANUOM-aiOICT  ANALYSTS  TOR  WIAK  SHOCK  TRANS I - 
HONS 

4 . 1  basic  Tquat ions 


The  analysis  is  based  on  the  following  as sump - 

t ions : 

(a i  The  flow  is  a  non stat ionary  one-dimensional 
( planar  or  spherically  symmetric)  viscous, 
c omp res s l b 1 e  air  f 1 ow . 

(b)  The  viscosity  ..  and  thermal  conductivity  are 
assumed  to  be  constant,  as  the  shock  waves  are 
weak . 

(c)  The  gas  is  assumed  to  be  thermally  perfect; 
the  equation  of  state  for  a  thermal ly-perfect 
gas  is  used. 

id)  Both  cases  of  calorical lv-perfect  and  imperfect 
gases  are  analysed.  Tor  calorical lv- imperfect 
cases  [ referred  to  as  real  gases),  the  vibra¬ 
tional  relaxation  of  air  molecules  are  taken 
into  account.  However,  for  most  cf  the  analyses, 
only  the  vibrational  relaxation  of  oxygen  is 
taken  into  account,  since  the  vibrat icnal -relaxa¬ 
tion  time  of  nitrogen  is  much  longer  than  the 
duration  of  most  N-waves  analysed  in  this  study. 
The  effects  of  nitrogen  vibrational  relaxation 
are  discussed  only  in  the  last  part  of  this 
section.  The  harmonic-oscillator  approximation 
is  applied  to  the  vibrational  energy  level. 


Then  the  basic  flow  equations  can  be  written  as: 


jU  +  *F 
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where  S  -  0  for  plane  flows  and  .i  =  2  for  spherical 
flows,  .  -  density,  v  -  velocity,  p  -  pressure, 

T  -  temperature,  I  -total  energy,  e  -  internal 
energy,  R  -  gas  constant,  -  vibrational  energy 
for  the  j -molecule  ( .i  -  U  for  oxygen  and  \  -  N  for 
nitrogen  i,  ;  • ,  )  c.  equ  i  I  ib  rj  urn  vibrational  energy 
for  the  j ■molecule,  tj  vibrational  relaxation 
time  for  the  '-molecule. 


Based  on  the  harmonic  oscillator  approximation, 
the  equilibrium  vibrational  energy  for  the  i  molecule 
(  ,)  can  he  expressed  as 


where  ,  is  the  characteristic  vibrational  tempera 
ture  for  tin-  i  molecule:  “0  -  22.v)  k,  \  =  3  3  SO  k, 
t.  j  molar  concent  rat  i on  for  the  i  molecule:  eo  - 
l>\201>,  -  0.-H1.  Du-  vibrational  temperature 

i  1\ ) j  for  the  i  molecule  can  also  be  defined  as 

R 

,  -  .Vpp/^rTT  i  1  ’  • 1 


e)  I'he  rotational  relaxation  is  taken  into  account 
through  the  bulk-viscosity  concept.  The  bulk 
viscosity  due  to  the  rotational  relaxation  is 
assumed  to  he  ,.r  =  (2/3);.. 


The  vibrational  relaxation  times  for  oxygen  to 
and  nitrogen  r\  are  evaluated  using  The  empirical 
relations  obtained  from  the  absorption  o4‘  sound 
waves  by  Bass  and  Shields  iRet.  23i,  as  follows; 
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*0  ,  . . 4,.  0.05  > 

—  24  +  4.4  x  10  h  -  - 

p  |  0.391  + 


\  =  4-  —  4-  1 9  ♦  3S0h  exp  i-b.142  3  -  1 

N  2"  P  ,k  1  o  1.  l  K  T 


where  pp  =  101.3  KPa,  Tq  -  293.15  K,  h  -  absolute 
humidity  of  air  (%) .  As  seen  in  Eqs.  (4.5)  and 
14.6),  the  vibrational  relaxation  times  for  oxygen 
and  nitrogen  strongly  depend  on  the  absolute  humid¬ 
ity  of  air.  In  Fig.  4.1,  Tq  and  Tv  are  plotted  as 
functions  of  the  absolute  humidity  for  p  =  101.3 
KPa  (Ref.  12).  The  relaxation  time  for  nitrogen  is 
two  or  three  orders  longer  than  the  relaxation  time 
for  oxygen.  The  relative  humidity  is  defined  as 

RH  =  h(p/psat)  (4.7) 

where  psat  is  the  partial  pressure  of  water  vapour 
at  saturation,  and  given  by  the  Goff-Gratch  equation 
(Ref.  24)  as 

,o*10<P»at/V 

=  10. 79586 [  1  - (T0/T)  ]  -  S.02808  log1()(T/T0J 

*  1.50474xl0-4l-l(r8-i9692^/ro)-ll' 

♦  0.42873X10'3  V- 78955  11 -Wl  -  , 


The  RCM  is  applied  to  step  (1)  and  the  explicit 
method  of  finite  difference  is  applied  to  steps  (2) 
and  (3).  In  step  (4),  the  integrated  relation  was 
used.  If  one  step  is  passed  over  among  steps 
(2)-(4),  then  the  following  solutions  result:  plane 
flow,  an  invisc id-nonequi 1 ibrium  flow  or  a  viscous- 
frozen  flow,  respectively.  These  are  termed  as  a 
plane  solution,  a  real -inviscid  solution  and  a 
perfect-viscous  solution,  respectively.  The  full 
solution  including  both  effects  of  vibrational  ex¬ 
citation  and  viscosity  is  called  a  real-viscous 
solution . 

An  outline  of  RCM  is  described  below.  Figure 
4.2  shows  an  illustrative  diagram  for  grid  con¬ 
struction  and  sequence  of  the  sampling  procedure. 

The  notations  r  and  t  are  increments  of  space  and 
time,  respectively.  For  arbitrary  integers  n  and 
i,  the  properties  lJh  +  '  at  time  ( n-w  1 ) '  t  are  cal¬ 
culated  from  the  properties  U”  at  time  n't.  The 
intermediate  values  arc  evaluated  at  time 

(n  +  l/2).'t.  In  the  region  of  i  r  r  (,i»l)’r 
and  n't  t  (n  »  1/2)  t  (surrounded  by  the  broken 
lines  in  Fig.  4.2),  the  Riemann  problem  (shock-tube 
problem)  is  solved  for  the  discontinuous  initial 
values 


4 . 2  Numerical  Method 

An  operator-splitting  technique  was  applied  to 
Fq.  (4.1).  The  calculation  is  done  for  each  spatial 
mesh  in  each  time  step  using  the  following  procedure: 

il)  The  hyperbolic  equations  are  solved  for  an  in¬ 
viscid  frozen  flow, 


where  the  subscript  I  indicates  the  solut ion  of 
step  ( 1) . 

12)  The  spherical  corrections  arc  made  by  using  the 
values  of  the  physical  properties  evaluated  in 
step  (1). 


13)  The  viscous  diffusion  equations  are  solved  by 
using  the  values  of  the  physical  properties 
evaluated  in  step  (2). 


♦  1  —  !  C,  -  J  ( H  )  , 
r  T  .  2  J  v  2 


in  this  region.  Then,  for  example,  the  solution 
consists  of  a  shock  wave  S,  an  expansion  wave  R 
and  a  contact  surface  as  shown  in  Fig.  4.2.  At 
tine  in  ♦  1/2) 't,  the  region  i'r  ■  r  (i  +  l).'r  can 
he  divided  into  four  subregions  ())-|4)  (or  five 
subregions,  if  the  interior  of  the  expansion  fan  is 
taken  into  account ) ,  and  the  physical  properties  in 
each  subregion  can  he  determined  from  the  solution 
cf  the  shock-tube  problem  for  the  initial  condition 

i  1.13).  l'he  values  are  equated  to  those  of 

ii"  1  "  at  point  P  (r  -  r  ;  i.'r  r  (i*ll'r|, 
which  is  chosen  at  ranuom.  That  is,  we  assume 
lln  !  ",  =  l/n  The  choice  of  I’  is  made  bv  a 

1*12  p 

random- samp) ing  technique  in  such  a  way  that  the 

sampled  points  are  uniformly  distributed  within  a 

finite-sampling  frequency.  In  a  similar  way,  the 

values  of  lln  !7~  are  obtained  from  the  initial 
i  -  1 /  2 

values  of  lln  and  ll"  At  the  second  half-time 
J  i  - 1  j 

step,  the  values  of  ll.  are  evaluated  from  the 

,  ..  ..11*1/2  ,  , ,n *  1  / 2  .... 

values  ol  ll  and  ll.  ,  ,  ,  as  initial  ones. 

1-1/2  i ♦ 1 / 2 

Godunov's  iterative  technique  is  applied  to  solve 
tile  Riemann  problem.  As  the  vibrational  energies 
are  assumed  to  he  frozen,  they  are  invariant 
across  the  waves,  and  keep  their  initial  values, 
whose  boundarv  is  the  contact  surface. 


(4 |  The  vibrational  relaxation  equations  are  solved 
by  using  the  values  of  the  physical  properties 
evaluated  in  step  (3). 


Tt~  -  “V.3 


The  final  solutions  are  obtained  in  step  (4). 


In  the  second  and  third  steps  of  the  operator  - 
splitting  technique,  explicit  finite-difference 
schemes  are  employed.  Flic  finite-difference  forms 
of  l.qs .  (4.10)  and  (4.11)  reduce  to 


'«v. ir1' 


, 


..  ,n»l  ,,  ,n  +  l  ,t  ,n»l  ,n  +  l  .n  +  1, 

1  Vi  i  i  *  2  1  ((V  i»l  “  •(L21i  ^2  i-1 

(  r) 

♦Tf-V  lit,)"*}  -  (C2)"*{]  -  jlUy.,]^1  2,t  (4.15) 

1 


Ihe  multiple  time  step  is  used  to  evaluate  (U3 )  ^ 
to  improve  accuracy.  At  the  intermediate  substep 


4.3  Solutions  for  Plane  Waves 


As  a  check  on  the  method,  the  one-dimensional 
shock -tube  problem  was  solved  for  a  perfect-inviscid 
flow,  perfect -viscous  flows  and  rea 1 -invisc id  flows. 
The  thickness  and  structure  of  the  shock  waves  are 
compared  with  those  obtained  analytically  in  Chapter 
3. 


-"'Viu.r''/"'  *  [(c,)n’('/k) 

•’  1  (  r)  - 


3'i  +  l 


j  -'-t/k  n. ( ../k ) 

2r .  ;.r  ‘ 1  3 J  i-1 


-v'":;  k,i  -  .-i«v3ir(‘/k)  t  ■  °> i— 


k) 
(4.16) 


where  the  time  increment  .t  is  subdivided  by  k. 

Most  of  the  calculations  were  carried  out  for  k=  10. 

In  the  4th  step,  the  vibrational  relaxation 
equations  for  air  molecules 


‘Vo 


(.1  =  *\  \) 


(4. 17) 


are  solved  in  each  spatial  mesh  under  the  assumption 
of  constant  temperature  and  pressure,  thereby  yield¬ 
ing  the  analytical  relation 


4.3.1  Perfect-inviscid  Solution 


Figure  4.3  shows  a  computer  plot  of  a  perfect- 
inviscid  solution  of  overpressure  (.\p)  against  dis¬ 
tance  x  for  several  time  intervals  for  a  diaphragm 
pressure  ratio  P4 \  =  2  and  initial  temperature  ratio 
T4I  =  1.  The  overpressure  (.\p)  is  normalized  by  the 
initial  pressure  pj ,  and  the  distance  x  is  normalized 
by  the  length  of  the  high-pressure  chamber  xq  (x*  = 
x/xq) .  The  diaphragm  is  placed  at  x*  =  1 .  The  time 
t  is  normalized  by  xq/hj  (t*  =  ajt/xy)  .  After  start¬ 
ing  the  calculation  or  the  removal  of  the  diaphragm, 
a  shock  wave  as  a  discontinuous  front  propagates 
towards  the  right  hand  side,  and  a  rarefaction  wave 
propagates  towards  the  left  hand  side.  When  t*  =  1 
the  head  of  the  rarefaction  wave  arrives  at  the  end 
wall  of  the  high-pressure  chamber.  The  shock  Mach 
number  Ms  is  about  1.16,  and  the  normalized  equili¬ 
brium  overpressure  or  the  shock  strength  ( ^ p ) 2 / P 1 
is  about  0.403.  It  should  be  noted  that,  unlike 
finite-difference  schemes,  the  shock  wave  as  a  dis¬ 
continuous  front  occupies  one  mesh  jump  without 
smearing,  where  the  normalized  one -mesh  size  \x*  - 
1/40. 


h  \ 


vr 


* .  1 

\  e  a 


n>  1 
i 


4.3.2  Perfect -Vi scous  Solutions 


,  n  ♦  1  t 

exp  -  — 
1  a 1  1  1 


(j  =  0,  N) 


14.18) 


i he  finite-difference  schemes  with  multiple  time 
steps,  similar  to  Fq.  (4.1b),  were  also  applied  to 
lq.  >1.1"),  and  found  to  give  the  same  results  as 
lv{.  i4.18).  In  order  to  reduce  the  computation  time, 
I q .  (1.18)  was  used  for  most  of  the  calculations. 

\ s  described  in  Section  4.1,  in  the  present  study, 
only  the  vibrational  relaxation  equation  for  oxygen 
was  solved  (except  Section  4.4.6).  Furthermore,  the 
hulk  viscosity  concept  was  applied,  instead  of  iiq. 

1.1"),  to  the  vibrational  relaxation  for  oxygen  in 
Sections  4.4.3  and  4.4.6,  in  which  the  N -waves  with 
long  durations  were  analysed. 


I  he  condition  of  symmetry  is  imposed  on  the  wall 
boundary  and  at  the  centre  of  the  sphere.  That  is, 
at  the  boundary  i^’r 

IF  ,  I  v )  =  11  t  l  -  v )  (4.19) 

0  3  0  1 


Ihe  condition  of  continuity  of  gradient  is  imposed 
on  the  free  boundary.  That  is,  at  the  free  boundary 
ll  r 


(  4 . 20 ) 


The  entire  programs  are  given  in  Appendix  C. 


Figure  4.4  shows  a  computer  plot  of  a  perfect- 
viscous  solution  for  the  same  case  as  Fig.  4.3.  The 
rarefaction  wave  reflects  at  the  end  wall  (x*  =  0) 
and  proceeds  towards  the  right  hand  side.  As  expected, 
smooth  shock  transitions  due  to  actual  viscosity  are 
obtained.  In  order  to  show  these  smooth  transitions 
clearly,  a  hypothetical  chamber  length  x  =  0.001  cm 
was  assumed  at  an  initial  pressure  and  temperature  of 
Pl  =  101.3  KPa  and  \'\  =  275.15  K.  Consequently, 
t  =  0.106  ..sec  for  t*  -  3.506. 

Here,  it  was  not  necessary  to  obtain  the  whole 
flow  field.  The  fine  structure  of  the  shock  front 
was  important.  Therefore,  in  order  to  save  computa¬ 
tion  time,  the  calculation  was  done  only  in  a  con¬ 
fined  region  near  the  front  for  the  wave  far  from 
the  diaphragm,  neglecting  the  behaviour  of  the 
rarefaction  wave.  Figure  4.5  illustrates  the 
region  of  calculation  and  a  plot  of  the  shock-front 
path  in  the  x*-t*  plane.  In  the  calculation,  30-80 
mesh  points  around  the  front  were  used,  and  the 
physical  properties  at  each  mesh  point  were  trans¬ 
ferred  hack  to  two  points  in  the  computational 
space  as  the  wave  proceeds  over  two  points  in 
physical  space.  The  condition  of  continuity,  Iiq. 
(4.20),  is  imposed  on  the  free  boundary  of  the 
region  of  calculation.  In  Fig.  4.5,  the  white 
circles  indicate  the  perfect -viscous  RCM  solutions, 
in  which  the  position  of  the  shock  front  is  defined 
as  the  po:  it  ion  of  50"*  of  (  p)j.  The  solid  and 
broken  lines  indicate  analytical  shock  and  sound¬ 
wave  paths,  respectively.  Ihe  numerical  solution 
for  the  shock  path  is  in  excellent  agreement  with 
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ana  lysis . 

Figures  4.6(a)-lc)  show  perfect -viscous  numeri¬ 
cal  solutions  for  the  shock- tube  problem  described 
above  by  comparison  with  Taylor's  and  Lighthill's 
analytical  solutions  for  the  shock  thickness,  which 
is  defined  by  10-90%  of  (/ip) 2-  The  ratio  of  the 
thickness  parameters  (AZ)  /  (/.Z)q,  which  corresponds 
to  the  thickness  or  rise  time  normalized  by  the 
Taylor  thickness  or  rise  time,  is  plotted  against 
the  time  parameter  r  defined  by  F.q.  (3.21).  The 
figures  indicate  that  the  step  wave  with  zero  thick¬ 
ness  is  reduced  to  a  plane  wave  with  a  smooth  transi¬ 
tion  owing  to  viscous  action,  as  the  wave  proceeds. 
The  broken  and  solid  lines  indicate  Taylor's  and 
lighthill's  solutions,  respectively.  The  various 
numerical  solutions  are  indicated  by  symbols.  All 
calculations  were  carried  out  for  the  same  case  as 
Fig.  4.4  [P41  =  2,  T41  =  1,  Tj  =  273  K,  p]  =  101.3 
KPa,  Ms  =  1.1b,  (/p)2  =  40.7  KPa] . 

Figure  4.0(a)  shows  the  effect  of  multiple  time 
step  for  viscous  correction.  The  mesh  size  is 
...\*  =  1/40  (/.x  =  2.5xlO“S  cm).  The  black  and  white 
circles  indicate  the  cases  for  k  =  1  and  10  in  F.q. 
(4.1o),  respectively.  The  k  =  10  result  for  the 
transient  behaviour  of  the  shock  thickness  is  closer 
to  lighthill's  solution.  It  is  seen  that  the  mul¬ 
tiple  time  step  for  viscous  correction  improves 
the  result  for  the  transient  behaviour  of  the  shock 
thickness.  The  random  walk  due  to  the  random 
sampling  in  RCM  and  the  overshoot  of  the  thickness 
value  above  Taylor's  value  can  be  seen.  The  mul¬ 
tiple  time  step  of  k  =  10  was  used  for  all  calcula¬ 
tions  described  below. 

Figure  4.6(b)  shows  the  effect  of  the  choice 
of  random  numbers.  The  mesh  size  is  \x*  =  1/80 
(  x  =  1.25x10”5  cm).  The  black  and  white  circles 
indicate  the  cases  using  the  random  numbers  by 
maximum-length  linearly  recurring  sequence  and 
linear  congruent ial  sequence,  respectively.  It 
can  be  seen  that  the  latter  method  is  in  better 
agreement.  Therefore,  linear  congruent ial  sequence 
was  used  for  all  other  calculations  in  the  present 
study,  as  well  as  by  Saito  and  Glass  (Ref.  1").  It 
is  also  seen  in  Fig.  4.6(b)  that  the  result  is  im¬ 
proved  by  reducing  the  mesh  size  by  half,  in  com¬ 
parison  with  the  result  in  Fig.  4.6(a). 

Figure  4.6(c)  shows  the  comparison  between  the 
RGM  and  MacCormack's  finite-difference  method  (MFM), 
which  is  shown  in  Appendix  D.  The  MFM  solution  is 
in  poor  agreement  with  Lighthill’s  solution.  Its 
thickness  or  rise  time  values  are  much  larger  than 
the  analytical  ones  owing  to  the  effect  of  art  if i- 
cial  viscosity.  The  RCM  solution  with  operator¬ 
splitting  techniques  is  superior  to  the  MFM  solution 
for  the  same  mesh  size,  although  random  scattering 
of  the  thickness  or  rise  time  values  do  occur. 

Better  agreement  with  Lighthill's  solution  was 
attained  by  using  a  finer  RCM  mesh  as  shown.  Com¬ 
puter  costs  would  limit  the  ultimate  mesh  size  to 
he  used . 

In  Fig.  4.^,  the  normalized  overpressure  (Ap)/ 
('p)i  is  plotted  against  the  distance  parameter  Z 
at  times  r  =  0.99,  45.0  and  58.3  for  cases  of 
x  =  1.25x10"^  cm  |white  circle  in  Fig.  4.6(c)].  The 
origin  of  Z  is  taken  at  the  place  of  (Ap)/(Ap)2  = 

0.5.  The  solid  lines  indicate  Lighthill’s  solution 
for  the  transient  state  at  :  =  0.99  and  Taylor's 
solution  for  the  final  steady  state  at  :  *  The 


RCM  pressure  profiles  show  very  good  agreement  with 
the  analyses.  This  result  suggests  that  the  RCM 
with  the  operator-splitting  technique  may  be 
applied  to  analyse  the  transient  behaviour  of  a 
viscous  shock  structure,  though  some  random  walks 
and  overshoot  above  the  Taylor  value  were  observed 
for  the  thickness  or  rise  time  data. 

4.3.3  Rea  1  - Inv i sc  id  Solution 

The  initial  conditions  ( P4  j  =  1.0018,  T4j  = 
1.0,  pi  =  101.3  KPa,  T 1  =  303.15  K  and  RH  *  90%) 
were  chosen  so  as  to  give  a  fully-dispersed  wave 
in  the  final  steady  state  for  a  real -invisc id  flow, 
and  to  give  a  fast  approach  to  the  steady  state 
in  order  to  reduce  the  computational  cost.  Only 
the  vibrational  excitation  for  oxygen  molecules 
was  taken  into  account  for  atmospheric  air.  The 
corresponding  relaxation  time  for  oxygen  is  T 0  = 

1.04  Msec  and  the  characteristic  time  using  the 
bulk  viscosity  (*v)n  **or  oxygen  is  (‘  v )  ()/ a  1  ^ 

8. 4x10“ 10  sec.  The  equilibrium  shock  Mach  number 
Me  =  1.0004  and  the  equilibrium  overpressure  is 
(\p)>  =  91.1  Pa,  which  is  less  than  the  correspond¬ 
ing  critical  overpressure  for  oxygen  (  p)cr  q  =  95-5 
Pa,  so  that  the  wave  may  become  a  fully-dispersed 
wave  in  the  final  steady  state. 

Figure  4.8  shows  the  transient  hehaviour  of 
the  pressure  and  temperature  profiles  of  the  dis¬ 
persed  wave  obtained  for  the  condition  described 
above  (x0  =  0.5  cm  and  ‘x  =  0.0125  cm).  The  solid 
lines  indicate  the  pressure  and  temperature  profiles 
which  are  the  same  in  normalized  plots  of  l'p)/(p)2 
and  (/T)/(‘T)2  as  the  wave  is  very  weak.  The  broken 
lines  indicate  the  normalized  vibrational  tempera¬ 
ture  profiles  (/Tv)/(‘T)j.  Ten  profiles  are  shown 
for  the  time  parameter  '  =■  0.0003,  0.41,  0.81,  1.6, 
3.3,  4.9,  6.5,  8.1,  9 .  ~  and  11.4  or  the  normalized 
distance  x*  -  1.2,  30,  60,  120,  238,  360,  476,  593, 
716  and  830,  where  *  is  defined  using  the  bulk  vis- 
cosity  (  •  V 1 0  for  oxygen  as  -  =  [a  1  -t  /  (  ;  v  )0)  [.'  p)  i/pi  ]  2 
The  calculation  was  also  carried  out  only  for  a 
confined  region  near  the  front  for  the  wave  far 
from  the  diaphragm,  similar  to  the  perfect-viscous 
flow  as  shown  in  Fig.  4.5.  The  initial  step  wave 
is  smoothed  out  owing  to  the  dissipative  effect 
of  the  vibrational  relaxation.  It  should  be  noted 
that  this  process  which  smears  the  wave  is  largely 
different  from  that  of  the  viscous  wave.  This 
tendency  of  smoothing  has  been  shown  analytically 
for  linear  waves  (Ref.  25)  and  for  nonlinear  waves 
(Ref.  2b).  In  a  transient  state,  the  wave  is  a 
partly-dispersed  wave  with  a  frocen  shock  front, 
even  if  the  equilibrium  shock  pressure  is  below 
the  critical  overpressure.  This  suggests  that  the 
nonstationary  effect  is  more  important  for  dispersed 
waves  than  for  viscous  shocks. 

Figure  4.9  shows  plots  of  ( '.Z  >  /  ( ‘  Z )  q  vs  '  for 
real -inviscid  shocks.  The  solid  and  broken  lines 
indicate  the  modified  I.ighthill  solution  and  the 
modified  Taylor  solution,  respectively.  The  symbols 
indicate  the  RCM  solutions  for  ,'x  =  0.025  cm  and 
0.0125  cm,  respectively.  The  latter  case  correspond 
to  the  one  in  Fig.  4.8.  The  RCM  solutions  of  shock 
thickness  show  random  walks  and  overshoot  above  the 
Taylor  value,  similar  to  the  viscous  solutions  shown 
in  Fig.  4.b.  The  thickness  tends  to  approach  the 
modified  Taylor  value  using  the  bulk  viscosity  for 
oxygen  vibrational  relaxation.  It  should  be  noted 
that  the  shock-th  ickcn ing  time  of  the  RCM  solution 
is  nearly  the  same  as  that  of  the  modified  L i gh thill 
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solution,  although  the  (.Mj ,  l.-'.Z vs  :  plot  of  the 
Ki'M  solution  deviates  from  Lighth ill's  solution 
owing  to  the  difference  in  the  t ransient -wave 
profiles  between  the  two  solutions  shown  in  Figs. 

5.11  and  4.8.  That  is,  the  shock-thickening  time 
based  on  the  modified  Lighthill  solution  provides 
a  reasonable  estimate. 

In  I'tg.  4.10,  the  normalized  overpressure 
i ,‘pi  i'pij  is  plotted  against  the  distance  parameter 
Z  at  -  25.0  and  27.0  for  the  case  of  .‘.x  =  0.0125 
cm  (white  circle  in  Fig.  4.0),  where  Z  is  also 
defined  using  the  bulk  viscosity  for  oxygen  vib¬ 
rational  relaxation.  The  solid  line  indicates  the 
analytical  solution  for  l  •  p)  V  ( '.p )  L-r  q  =  0.954 
evaluated  from  l.q.  (,5.54)  for  steady  dispersed  waves. 
1'he  ROM  pressure  profile  for  '  =  27.(>  shows  very 
good  agreement  with  anal>sis,  hut  the  ’  =25.0 
solution  shows  a  slight  deviation  from  the  analyti¬ 
cal  one  at  the  upstream  side  of  the  front.  This 
deviation  would  be  attributed  to  the  randomness 
associated  with  the  KOI  solution.  However,  in 
general,  the  ROM  solution  for  real -inviscid  flow 
;  t  'vules  very  reasonable  results. 

1.1  Solutions  for  -Spherical  haves 

\ s  described  in  i.hapter  5 ,  the  shock  structures 
of  spherical  waves  may  be  affected  by  N-wave  and 
nonstationary  effects  and  would  be  different  from 
those  of  plane  waves  in  some  situations.  The  pur¬ 
pose  of  this  section  is  to  show  some  characteristic 
features  of  transient  behaviour  of  shock  structures 
of  spherical  waves  through  the  ROM  analysis  associ¬ 
ated  with  the  spark  and  explod ing-w i re  experiments. 

Iwcnty-three  cases  of  numerical  results  are 
presented  in  this  section  for  spherical  waves,  and 

termed  as  cases  41,  \2 .  HI,  B2,...,  (1,  C2 . 

1 1 1  ,  I»J .  respectively.  The  A-series  lAl,  A2....I 

corresponds  to  perfect  - inv isc id  solutions;  B-series, 
perfect -viscous  solutions;  C- series,  real -inv isc id 
solutions;  and  H- series,  real -viscous  solutions. 

Ihe  parameters,  which  should  be  given  as  initial 
conditions,  are  the  radius  of  the  pressurized 
sphere  rq,  the  pressure  and  temperature  ratios 
I'll  and  1 4 1  across  the  initial  inner  pressurized 
air  and  the  ambient  atmosphere,  the  atmospheric 
pressure  pi  and  temperature  Ij,  ami  the  relative 
humidity  KH.  These  are  tabulated  for  each  case 
in  l able  4.1.  We  assumed  p]  =  101.5  KPa  for  all 
cases.  The  relaxation  time  ,  and  the  spatial 
meshes  r*  and  r  are  also  tabulated  in  Table  4.1. 

The  atmospheric  conditions  ilj  and  RH)  are  chosen 
from  data  in  the  spark  and  exp lod ing-w i re  expert 
ments  described  in  ihapter  2  (Series  1  IV). 

In  the  C  and  P-series  analyses  (real  gases  i  ,  the 
vibrational  excitation  is  taken  into  account  onl\ 
for  oxygen  except  case  PS,  in  which  both  vibrational 
excitations  for  oxygen  and  nitrogen  are  included. 

In  cases  P5  through  PS,  the  vibrational  relaxation 
for  oxygen  is  evaluated  by  using  the  bulk  viscositv 
concept  instead  of  solving  the  relaxation  equation 
for  oxygen. 

1.1.1  Ne.i  r  - 1  i  *.•  1  d  Solutions  for  lYrfVct  Inviscid 
I  lows 


!:i  this  stvf  i  or. .  perfect  -  inv  I  sc  id  solutions 
lot  spberual  wave''  are  shown  in  the  near  field  of 
the  pre  < '-u  t  i  ;«'d  sphere.  The  process  of  \  wave 


format  ion  from  an  explosion  of  a  pressurized  air 
sphere  and  the  effects  of  the  pressure  and  temper¬ 
ature  ratios  are  discussed. 

Case  A1  is  a  perfect- inviscid  solution  for 
1*41  -  2  in  the  near-field  of  a  pressurized  sphere. 
Figures  4. 11 (a) -(c)  exhibit  computer  plots  of 
overpressure  distribution  at  various  times  after 
the  exp  I  os* on. 

Figure  4 . 1 , . a  i  shows  the  initial  process  of 
explosion  of  a  pressurized  air  sphere.  The  front 
shock,  which  is  formed  immediately  after  bursting 
the  sphere,  decays  as  it  propagates  outwards, 
leaving  an  expanding  low  behind  it.  The  rare¬ 
faction  wave,  which  propagates  inwards  into  the 
sphere,  reflected  at  the  centre  of  the  sphere  and 
produces  a  highly  rarefied  region  behind  it.  A 
second  imploding  shod  wave  of  ever  increasing 
strength  i>  formed  at  the  boundary  between  the 
inner  and  outer  expansion  regions.  Some  "noise" 
in  the  pressure  profiles  in  the  expansion  region 
can  be  attributed  to  the  random  walk  inherent  in 
the  RCM.  Ihe  comparison  between  near-field  solutions 
of  the  explosion  of  a  pressurized  air  sphere  using 
Lax,  MacCormaek  and  Random-Choice  methods  for  a 
perfect -inviscid  flow  is  given  in  Appendix  T . 

The  succeeding  process  of  N-wave  formation  is 
shown  in  Fig.  4.11(b).  The  imploding  second  shock 
reflects  at  the  centre  of  the  sphere  and  produces 
a  highly  compressed  region  around  it.  The  reflected 
second  shock  is  initially  very  strong,  but  rapidly 
decays  at  it  proceeds  outwards.  It  follows  the 
front  shock  and  forms  the  rear  shock  of  an  N-wave. 
Figure  4.11(c)  exhibits  the  propagation  of  an 
established  N-wave,  which  maintains  a  similar 
profile  as  it  propagates  outwards.  Its  maximum 
i peak)  overpressure  decays  gradually  and  the  dur¬ 
ation  increases  slowly. 

Figures  4.12ial-idi  show  a  comparison  of 
established  N- waves  for  cases  A1-A4.  Figure  4.12(a) 
exhibits  a  pressure  profile  for  the  same  case  as 
Fig.  4.11(c),  though  the  mesh  size  is  increased  to 
‘ r*  -  1,10  to  be  compared  with  cases  A2-A4.  In 
case  A2,  the  temperature  ratio  T.j  j  is  twice  that 
for  case  A1  .  In  case  A  5 ,  the  pressure  ratio  P.j]  is 
increased  from  2  to  9.  In  case  A4,  both  pressure 
and  temperature  ratios  are  higher.  Figure  4.12(b) 
shows  that  case  A2  results  in  a  more  symmetric 
N-wave  than  case  A1  owing  to  the  hotter  sphere, 
which  enables  the  second  shock  to  form  sooner. 

This  suggests  that  the  half-duration  of  the  nega¬ 
tive  overpressure  of  an  N-wave  can  be  controlled 
through  a  choice  of  I p  .  Figures  12(c)  and  12(d) 
show  that  for  higher  I'.j  i  and  T4  \  the  S -waves 
generated  by  a  spark  or  an  exploding  wire  cannot 
be  simulated  using  a  pressur i zed-sphere  explosion 
mode  1  . 

As  seen  in  Tigs.  4. 12(a) -(d),  the  overpressure 
profiles  of  the  positive  phase  show  only  a  slight 
change  in  shape  regardless  of  P.jj  and  \.\\  (although 
p  and  t  lit*  durations  are  different).  However,  the 
negative  phases  strongly  depend  on  these  ratios. 

Ihe  length  of  the  positive  side  is  of  the  order  of 
r*  1  or  r  ■  ro  in  each  case.  That  is,  the  half¬ 
duration  of  an  \  wavi  is  determined  mainly  by  a 
choice  of  the  sphere  radius  within  the  range  of 
i’ll  and  [4 1  considered  here.  In  the  following, 
use  is  made  of  1 ; \  1,  in  order  to  simplify  the 

analysis,  since  attention  is  focussed  on  the  front - 
shock  structures  of  the  \  waves  in  this  work. 


1 . I . J  Comparison  Between  Perfect-lnvisciJ,  Perfect  - 

V 1  sc  ou  s ,  ReaJ-  In  \i  sc  id  and  Real -Vi  sc  ou  s  , 

I -a  r  - F  1  e  1 d  So l u 1 1 on s 

l‘hc  calculation  for  cases  A5,  B1  ,  ('1  and  1)1 
were  carried  out  tor  the  same  parameters  in  order  to 
make  the  comparison  clear  between  perfect  -  inv  i  sc  id , 
perfect -v  iscous ,  real  - inviscid  and  real -viscous 
Si  hit  ions  in  the  far  field.  1'he  vibrational  excita¬ 
tion  for  oxygon  was  taken  into  account  for  real 
cases  il  l  and  1)1).  The  ambient  conditions  corre¬ 
spond  to  the  series- l  experiment,  and  the  relaxation 
t  i me  q  =  15.0  ..sec 

Ihe  results  are  shown  in  Figs.  4.15-4.1’’. 

Figure  4.15  shows  the  path  of  the  shock  front  by 
plotting  the  centre  of  the  front  [ 0 .3 ( ‘p)max ] .  The 
normalized  radius  r*  and  the  normalized  time  t*  are 
defined  by  r*  =  r/r,j  and  x*  -  ajt  'r(),  respectively. 
Ihe  sol  id  line  indicates  the  path  of  a  sonic  line. 

It  i>  seen  that  away  from  the  explosion  the  front 
path  nearly  coincides  with  the  sonic-line  path, 
lius  result  indicates  tfu  validity  of  the  method  of 
solution  with  regard  to  the  propagation  of  the  wave. 
Ihe  calculations  were  also  carried  out  only  in  a 
vonfineJ  computational  region  near  the  front  in  the 
tar  field  as  well  as  the  calculations  for  plane 
wa i es  shown  i n  F i g .  I . 5 . 

Ihe  maximum  o\ erpressuie  i  pnux I  for  spherical 
waxes  deca\  with  increasing  distance  r  from  the 
centre.  According  to  classical  acoustic  theory 
.  p:!;-iix  r  1  for  weak  spherical  waves.  However, 
as  shown  in  the  following,  the  decay  of  the  maximum 
overpressure  can  deviate  from  classical  theory  if 
the  effects  of  viscosity  and  vibrational  nonequ  i  1  i  - 
brium  are  taken  into  account  .  In  order  to  roadilv 
evaluate  the  decay  rate  of  the  maximum  overpressure, 
the  Jecuv  index  n  i>  /ntroJuceJ,  where  n  is  defined 
locally  as  p>nux  r  n.  In  general,  the  value  of 
n  varies  with  r,  while  n  -  1  applies  to  spherical 
ac oust ic  wavi > . 

figure  111  shows  the  decay  of  the  maximum 
overpressure  for  four  cases  as  a  funct ion  of  the 
distance  r.  In  the  perfect  -  invisc id  solution  (case 
\."> ) ,  the  maximum  overpressure  decays  at  a  rate 
inverselv  proportional  to  r  in  -  1)  for  (’p)max 
loo  Pa,  though  n  I  for  i  pim  100  Fa.  In 
other  cases,  HI  ,  1*1  and  PI,  the  decay  indices  n 
increase  for  (p),nax  l*h)  Pa  due  to  the  dissipative 
effects  of  viscosity  and  vibrational  nonequ i 1 ibrium 
m  comparison  with  case  A5.  While  almost  the  same 
overpressures  are  obtained  for  i  pimax  100  Pa  for 
all  cases  including  case  A3,  at  r  =  20m,  n  =  1.25 
for  case  HI  and  n  =  1  .4l>5  for  cases  Cl  and  ol  .  The 
deviation  from  the  classical  acoustic  theory  for 
(  ’"max  H'O  Pn  is  attributed  to  r hr  nonlinear 
effects  in  a  wave  of  finn«  ...  •  1  i  *  n  *  1 » ■ . 

f  igure  4.15  exhibits  the  ha  1 f  jurat  ion  tj  as  a 
function  of  distance  r.  The  rapid  increase  of  tj 
near  the  centre  is  attributed  to  nonlinear  effects. 
In  case  \5,  tj  is  constant  for  (p)max  100  Pa, 
while  tn  cases  HI,  Cl  and  f > ! ,  tj  increases  with  r 
di i e  t  o  d  i  s  s  i  pa  t  i  ve  e  f  fee  t  s  o f  v  i  s c  os  i  t  y  a  n d  v  i  b  ra  - 
tional  nonequ  i  I  ibruim. 

Figure  1 . In  shows  the  rise  times  tr  as  a  func¬ 
tion  of  distance  r,  and  Fig.  4 . I "  shows  the  pressure 
profiles  at  several  locations  for  cases  A5,  HI,  Cl 
and  01.  The  perfect - invisc id  solution  results  in 
a  discontinuous  front  so  that  tr  =  0  in  this  case, 


unlike  the  smoothing  causes  by  artificial  viscosity 
in  finite-difference  methods.  As  seen  in  case  Cl, 
the  effect  of  vibrational  nonequilibrium  contributes 
to  tr  only  for  weak  waves.  The  rise  times  for  the 
real -viscous  case  PI  arc  almost  the  same  as  the 
rise  times  of  the  perfect -viscous  case  HI,  until 
the  effect  of  vibrational  nonequilibrium  becomes 
noticeable.  The  viscous  effect  plays  a  dominant 
role  in  determining  the  rise  time  in  these  cases. 
However,  the  vibrational  nonequ i ) ibri urn  plays  an 
important  role  in  reducing  the  maximum  overpressure. 


The  profile  of  the  pci  feet -viscous  transition 
at  r  =  21.6m  |Fig.  4.17(1))]  is  not  similar  to 
either  the  profile  for  a  steady  plane  wave  Section 
5.1),  the  quasi -stat ionary  N-wave  for  moderate 
Reynolds  number  (Section  7.2),  or  the  nonstationary, 
plane  wave  (Section  5.2).  This  shows  a  character¬ 
istic  feature  of  the  nonstationary  effect  for 
spherical  N- waves.  Figure  4.1  "(c)  indicates  that 
the  wave  is  a  part ly  -dispersed  wave  with  a  discon¬ 
tinuous  front,  even  though  the  stead}'  plane  wave 
becomes  a  ful ly ~d i spersed  wave  wit)  a  smooth  tran¬ 
sition  for  the  corresponding  overpressure  at  r  = 
Jl.om  (Section  5.4  i.  Again,  this  is  a  nonstat ionary 
effect  for  dispersed  waves,  which  is  discussed  m 
Section  4.5.5.  Ihe  nonstat ionan  dissipative 
effects  due  to  viscosit}  and  vibrational  nonequi li¬ 
brium  are  coupled  in  the  real  -  \  i -icons  solution 
|  Fig.  l.l"ul)|. 


Ihe  results  for  cases  A3,  HI,  Cl  and  PI  show 
that  the  decay  behaviour  of  the  maximum  overpressure, 
the  ha  1 f -du rat  ion ,  the  rise  time  lor  the  shock  thick¬ 
ness!  and  the  pressure  luot’ilo  of  a  weak  spherical 
N-w;ne  can  be  affected  by  both  xiscosity  and  vibra¬ 
tional  nonequ  1 1  ibrium.  Hus  shows  that  both  effects 
must  he  taken  into  account  when  analysing  the  shock 
structure  of  a  weak  spherical  N-wave. 

1.4.5  Simulations  foj  Spark  and  1  xpl nd ing -W i re  - 

Hone rated  N -waves 

In  this  section,  the  numerical  simulations 
are  shown  for  the  spherical  N- waxes,  which  were 
genera t ed  from  spa rk  and  explod i ng -wire  sources , 
described  in  Chapter  2.  A  requirement  was  set 
for  the  *.  of  weak  spherical  N-wnvcs  that 

the  calculated  maximum  overpressure  I  and 

the  ha  1 f  - jurat  ion  tcj  should  coincide  with  the 
experimental  values  at  a  specified  location  r.  This 
requirement  can  be  fulfilled  by  giving  appropriate 
values  to  the  initial  pressure  ratio  Pj \  and  the 
radius  of  the  pressurized  sphere.  However,  in 
practice,  the  adjustment  of  the  values  of  P4 j  and 
rp  is  a  laborious  task  in  order  to  match  required 
values  of  (  pimax  an^  ld  at  11  spec  i  lied  location. 
Several  trial  calculations  were  needed  to  get  the 
final  result.  Cases  H2,  C2,  PI,  P2  and  P5  are 
the  results  of  .*  V. '  v  for  the  spark  and 
exp  1 od 1  ng  - w i re  da t  a . 

In  Figs.  4.18-4.21,  the  results  of  the  numerical 
calculations  are  compared  with  the  experimental  data 
by  plotting  ('p)max  vs  r,  ttj  vs  r,  tr  vs  r,  and 
tr  vs  l  pimax-  1h  these  figures,  the  experimental 
points  are  plotted  by  white  symbols  and  the  numerical 
ones  by  blacl  symbols.  The  solid  and  broken  lines 
denote  the  interpolated  lines  for  the  numerical  and 
experimental  data,  respectively.  In  Fig.  4.21,  the 
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broken  1  mes  denote  the  I.ighthill  rise  times  tor  N 
waves  with  tj  50  and  "0  ..sec,  and  the  chain  line 
denotes  the  fay  lor  rise  time  tor  steady,  plane  waves. 
I  he  abrupt  changes  in  rise  lime  are  attributed  to 
t‘ne  randomness  inherent  *n  the  Rl'M. 

I  i  cures  4 .  1 S  and  4.1*>  show  that  one  can  simu] 
at-  tile  change  of  and  td  against  r  hy  a 

proper  choice  ot'  vo  and  I’.j  \  .  fuvves  BJ  and  t'J  in 
t  ics.  1. Jo  and  l.Jl  indicate  that  the  perfect - 
viscous  ,BJi  and  rea  1  -  im  i  sc  id  U'Jl  solutions 
c.innot  simulate  the  nonstationary  behaviour  of  the 
ri'O  true,  even  if  i‘plnux  and  td  can  he  simulated 
i,  ntist  r.  turves  hi,  MJ  and  MS  indicate  that  the 
real  viscous  solutions  simulate  the  experimental 
ia salts  reasonably  well,  when  one  considers  t In¬ 
flow  complex  1 1  ies  .it  the  spark  discharge  and  exp  1  inl¬ 
ine  wire.  turves  Ml,  MJ  and  Mb  almost  simulate  the 
slat  a  for  series  l,  l  l  and  IV.  respectively. 

Ihi  general  features  of  the  results  can  be 
i i r i .  ed  as  follows: 

!he  dec  a  v  index  of  p':;;av  evaluated  from  the 
M'i’U't  1  and  II  spark  data  is  about  1  .45,  while  for 
*:u  " :  ..ul.ited  spark  case  Ml,  n  -  1.40  at  r  -  Jl.om 
and  u  -  l.Ul  at  r  I'd:::.  lhe  deviation  of  the 
n  a  1  in-  fro;1:  t!ie  linear  acoust  ic  theory  m  -  1)  is 
i  ...u  1  .»t  ed  reasonably  well  and  can  he  attributed  to 
tin  dissipative  effects  ot‘  viscosity  and  vibrational 
:  !  i  b  :  i  ;t.  for  uxy  gen  . 

die  half  duration  t  j  ir.crea.se>  with  r  .fig. 

I .  a  .  it-s  rate  of  increase  is  about  1  sec  :v. 
uid  simulates  the  experimental  data.  I  he  increase 
of  tj  may  a  1  so  be  attributed  to  dissipative  effects 
of  viscosity  and  vibrational  noneijtii  1  ibrium  for 
e\>  gen . 

5  lhe  spark  data  iserics  1  and  111  and  their 
simulation  for  rise  time  (cases  Ml  and  MJ )  show 
that  the  shorter  relaxation  time  t.5.5  .sec,  case  MJ) 
gives  tin  longer  rise  time. and  the  longer  relaxation 
tic  "Cv.  ,  case  Ml  I  gives  the  shorter  rise 

time.  Ibis  tendency  is  due  to  the  nonstat lonary 
eldest.  I  lie  long  relaxation  tint  gives  tin  slow 
rate  ot  change  of  the  shock  thickness,  as  discussed 
m  Wfvtions  5 . I  and  4.5.5,  m*  that  »he  rise  time 
remains  short  even  for  weak  waves,  lurther  discus¬ 
sion  about  the  effect  of  relaxation  time  on  t  ,•  will 
he  ijin'ii  m  the  succeeding  sect  ion. 

I1  lhe  exploding  wire  data  (series  IV)  and  their 
simulation  'case  M5 )  shew,  by  comparison  with  the 
spa  r k  data,  that  the  stronger  explosion  and  longer 
viur.it  ion  give  the  longer  rise  time  for  the  same 
overpressure  ll-ig.  l.Jl'.  ihis  is  again  due  t  o 
the  non--.!  ,i  t  i  Pii.i  ry  effects,  lhe  strong  explosion 
gives  a  slower  rate  of  change  of  tin*  maximum  over- 
pi1  s-siire  for  the  same  overpressure  (see  fig.  4.18) 
so  that  the  rise  time  has  enough  time  to  increase, 
furthermore,  a  longer  duration  provides  a  margin 
tor  increasing  the  rise  time.  Ihis  effect  will  he 
dismissed  m  more  detail  in  Section  4.4.5. 

1  igures  l  . JJ(a  i  - ( c }  show  the  pressure,  tempera 
Tuia-  and  v  ib  rat  ioiia  1  -temperature  profiles  at 
several  locuti  ns  for  real  viscous  cases,  Ml,  MJ 
and  MS,  respe  ‘ively.  ITie  solid  lines  indicate  the 
pressure  and  temperature  profiles,  which  are  the 
same  for  weak  waves  in  normalized  forms  of  (  p) 

p  i ... ,  and  '1  )  /  t  l  ),mix  •  lhe  broken  lines  indicate 
rfie  v  1 1)  ra  t  i  oim  1  temperature  in  a  normalized  form  of 


figure  l.JJiai  shows  the  simulation  for  the 
spark  iseries  1 i  experiment.  lhe  wave  profiles  a  re 
shown  at  r  -  M.J"m,  J.llm,  ".dim,  4.”Sm,  lJ.(>4m, 
15.<>m,  lR.oin  and  Jl.tirn,  and  the  maximum  overpres¬ 
sures  are  1  <><>.’  Pa,  *  >5  Pa,  50.0  Pa,  J5.0  Pa,  l-."' 

Pa,  15.J  Pa,  1 1  * .  1  Pa  and  S . V  Pa,  respectively.  As 
seen,  the  peak  pressure  and  T emperat ure  become 
gradually  blunted  due  to  the  energy  transfer  from 
the  translation.il  and  rotational  modes  to  the 
vibrational  mode,  while  t lie  back  (expansion) 
pressure  and  temperature  profile  becomes  gradually 
rounded  due  to  tlie  reverse  energy  transfer  from 
the  vibrational  mode  to  t  lie  t  runs  1  at  i  >na  1  and  rota¬ 
tional  modes.  Ibis  arises  owing  to  slow- tv  lax ing 
behaviour  of  the  vibrational  energy  and  leads  to 
an  elongation  of  the  half  duration.  lhe  shock 
thickness  or  the  rise  time  is  mainly  controlled  by 
the  dissipative  effect  of  viscosity,  though  it  is 
only  partly  affected  by  the  energy  t ran<fer  from 
the  translational  and  rotational  modes  to  the 
vibrational  mode  for  waves  at  r  ivom  and  Jl.om. 

In  a  sense  that  the  shock  front  is  mainly  controlled 
by  viscous  dissipation,  these  wave4*  may  be  called 
partly  dispersed  waves. 

figure  J.JJib  exhibits  the  simulation  for  tin 
spark  (series  11)  experiment.  lhe  wave  profiles 
are  shown  at  r  -  O.jOm,  l.aOn,  J.aSm,  t>.J"i:. ,  S.5~"., 

1  u  .  t » "*  m ,  IJ.'.im,  1  .  1  in  and  14.0m,  and  the  maximum 
ov  erpressures  are  lonJ  Pa,  151  Pa,  .ii’J1  Pa,  J5.0  Pa, 
1  ~ ~  Pa,  15.J  Pa,  10.4  Pa,  S..»  Pa  and  0.55  Pa, 
respectively,  approximately  in  accordance  with  the 
maximum  overpressures  in  fig.  4  .  J  J i a  » .  Hie  differ¬ 
ence  in  profiles  bet  ween  cases  Ml  and  MJ  can  readily 
be  seen.  In  case  MJ,  the  process  of  peak -blunt  ing 
occurs  between  r  -  !  .  am  and  4.5Sr:  or  (  pinj;ix  -  a  a .  1 
Pa  and  151  Pa,  while  in  case  Ml  it  occurs  between 
r  ” .  54  m  and  J1  .bin  or  >  pl|litlN  -  .  5  Pa  and  55.1 

Pa.  lh  i  s  can  he  attributed  to  the  differences  in 
vibrational  re laxa t  i on  time  for  oxygen  and  initial 
temperature:  ,  !.■.<•.  sec,  lj  _“5  k  in  case  Ml, 

and  a. ('..sec,  \\  JS4  y  m  case  MJ.  In  case 

of  the  shorter  rel.ixat  ion  time,  the  peak  blunting 
occurs  in  the  earlier  stage  when  the  shock  thickness 
is  still  relatively  thin,  lurthermore,  m  case  of 
the  higher  initial  temperature,  more  energy  is 
required  to  excite  the  vibrational  mode  so  that 
the  effect  of  vibrational  excitation  appears  for 
waves  at  higher  maximum  overpressure.  In  the 
range  r  -  S  .  5~m~ 1 4  .  On:  or  (  pln>.,x  o.]“  Pa  1  "  .  J  Pa, 
the  front  structures  are  mainly  controlled  by 
vibrational  excitation  and  the  wave  profiles  nearly 
follow  the  vibrational  temperature  profiles  owing 
to  tlie  energy  transfer  t  o  the  vibrational  node.  In 
this  sense,  the  waves  may  be  called  fu 1 1 y • d l spersed 
waves  in  this  range.  However,  it  should  be  noted 
that  the  viscous  dissipation  also  plays  an  important 
role  in  increasing  the  shock  thickness  or  rise  time, 
by  contrast  with  ''toady,  fully  dispersed,  plane 
waves,  as  seen  in  figs.  1 . jo  and  l.Jl  in  which  the 
real  -viscous  solutions  M*Jl  are  compared  with  the 
perfect  viscous  and  real  mviscui  solutions  (BJ,  fJ) 

figure  J.JJlc)  shows  the  simulation  for  the 
exploding  wire  (senes  i\  1  experiment.  The  wave 
profiles  are  shown  at  r  f 1 .  J  s : . .  4.15m,  ('.5m,  15.0m, 
IS. 5m,  Jl.lm  and  J4.5m,  and  the  maximum  overpressure 
are  lb'O  pa ,  155  Pa,  XJ.l  pa,  55.4  Pa,  JS.n  Pa,  1 " . " 
Pa  and  11.1  l'a  ,  respectively  .  lhe  waves  below  55 
Pa  show  characteristic  features  of  fully-dispersed 
waves,  though  the  shock  thicknesses  are  different 
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from  the  ones  in  lig.  4.22(b).  I  ho  peak  -b  hint  i  ng 
occurs  between  r  =■  1.13m  and  h .  3  in .  I  inures  4.22ui- 
c  show  that  wide  variations  of  wave  profiles  are 
possible  depending  on  combinations  of  relaxation 
time,  initial  temperature,  ha  1 f-durat ion ,  and 
strength  of  explosion. 

In  i igs.  4.23  and  4.24,  the  calculated  pressure 
profiles  are  compared  with  those  observed  at  several 
locations  for  series  l  and  11,  respectively,  figure 
1.23  shows  a  comparison  between  case  Ml  and  series  I 
profiles  at  r  -  4.S5m,  13. tun  and  21. tun,  while  1  ig. 

1.21  shows  a  comparison  between  case  M2  and  series 
II  profiles  at  r  -  1 1  . "m  and  19m.  Ihe  solid  and 
broken  lines  indicate  the  numerical  and  experimental 
pressure  profiles,  respectively.  Ihe  shock  transi¬ 
tion  and  overall  profiles  are  simulated  reasonably 
well,  if  we  con>Kkr  the  difficult)  of  adjusting  j 
and  r,i  to  get  the  required  values  for  i  plniax  and  tvj 
.it  the  des  i  red  positions.  ■.  v'onq\»  re  with  figs.  3.10 
and  3. Id  for  the  Iighthill  and  modified  l.ighthill 
pressure  profile*. 

In  I igs.  1.23  1.2",  the  full  N  wave  profiles  of 
pressure,  temperature  and  vibrational  temperature 
aii  plotted  at  the  longest  distances  of  observation 
m  series  I,  11  and  1\  for  cases  MIA,  M2A  and  M3  A . 

In  order  To  save  computation  time,  the  full  N-wave 
solution*  were  obtained  with  larger  mesh  si;:es, 
which  were  two  or  three  times  as  large  as  those  for 
the  half  \  wave  solutions  for  Ml,  M2  and  M3  shown 
above.  Ihese  figures  show  that  the  transition 
profiles  and  rise  times  of  the  rear  shocks  are 
different  from  the  front  shock  due  to  the  difference 
in  vibrational  nonequ i l ibrium. 

In  lig.  4.2S,  the  calculated  full  N-wave 
profile  of  pressure  is  compared  with  the  observed 
one  at  2d. 3m  for  series  IV,  in  which  the  full  V- 
wave  profiles  were  obtained  (Kef.  b) .  Although  the 
calculated  ha l f-durat  urn  tj  is  20,  longer  than  the 
observed  one,  both  profiles  are  similar.  The 
precise  simulation  for  full  N- waves  would  require 
an  adjustment  of  the  initial  temperature  ratio  T4 j 
in  addition  to  the  finer  adjustment  of  l\j  \  and  . 

This  may  be  done  in  a  future  study. 

4.1.4  iffects  of  Vibrational  Relaxation  lime 

Ihe  purpose  of  this  section  is  to  show  the 
effect  of  vibrational  relaxation  time  more  clearly 
by  comparing  cases  M2  and  Ml.  The  initial  pressure 
ratio  P4 j  and  the  sphere  radius  rp  are  the  same  for 
both  cases,  but  the  initial  temperature  and  humidity 
are  different  and  gives  rise  to  a  p  of  5.54  and 
13. b  ..sec,  respectively.  (The  initial  temperature 
and  humidity  of  case  Ml  correspond  to  case  Ml.) 

figures  1. 2d  4.31  show  (*pln,ax.  *d  and  tr  as 
functions  of  r.  fhe  discontinuous  change  of  tr  in 
1 ig.  4.31  is  also  attributed  to  the  randomness 
which  appears  in  the  RFM  solutions,  figure  4.32 
shows  a  comparison  of  both  pressure  profiles  at  the 
same  distance  r  1dm. 

The  attenuation  behaviour  o!*  1  o ',mr, Y  (Fig. 

1.29)  is  slightly  affected  by  the  vibrational  relax¬ 
ation  times  for  these  cases.  The  decay  indices  at 
r  =  20m  are  n  =  1.31  for  M2  and  n  =  l .4u  for  Ml. 

As  easily  seen  in  Figs.  1.29  1.31,  (  p)m;ix  ,?u^  *d 
are  not  affected  appreciably  by  the  difference  in 
■  y,  but  tr  is  very  much  affected  by  it.  The  rise 
tune  for  case  !>4  with  a  longer  relaxation  time  is 
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shorter  at  .1  fixed  distance  than  the  rise  time  for 
cast-  M2  with  .1  shorter  relaxation  time.  This 
tendeiicv  can  be  explained  b>  the  nonstationary 
eft  ect 

lot  steady,  dispersed  plane  waves,  as  shown 
in  Section  3.4,  tile  longer  relaxation  time  gives  a 
thicker  trail-- it  ion  or  a  longer  rise  time,  since  tr 
is  proportional  to  p  | see  I q .  «3 .42  1]  However, 
the  shock  thickening  tune,  which  was  defined  by 
the  time  of  approach  of  an  impulsive  step  wave  to 
the  final  steady  state,  is  also  proportional  to 
in  the  modified  I ighthill  solution  for  a  non station 
ary  fully  dispersed  wave,  furthermore,  it  was 
shown  in  Section  4.3.3  that  this  shock  thickening 
tune  was  in  close  agreement  with  the  RFM  solution. 

That  is,  the  longer  the  relaxation  time,  the  slower1 
is  the  rate  of  change  of  shock  thickness,  for  case 
M4  with  the  longer  relaxation  time,  the  wave  still 
remains  a  part  1 y  d i spersed  wave  whose  shock  transi¬ 
tion  is  mainly  controlled  by  viscous  action,  while 
for  case  M2  with  a  shorter  relaxation  time  the  wave 
becomes  a  ful 1 y -d i spersed  wave  whose  transition  is 
mainly  controlled  by  vibrational  nonequ i 1 ibr ium . 

Hus  is  the  reason  why  the  longer  relaxation  time 
gives  us  the  shorter  rise  time  for  the  weak  spheri¬ 
cal  waves  in  contrast  with  steady  plane  waves. 

1.4.3  I :  f f cc t s  of  N - U a vc  l H 1  ra t  ion 

In  t-i-  section,  the  effects  of  the  duration 
of  the  N-wave  on  the  decay  rate  of  I'plntux*  the  rate 
of  increase  of  tj  and  the  rise  time  tr  are  inves¬ 
tigated  by  changing  the  radius  i*u  of  the  pressurized 
air  sphere.  in  cases  M3  and  Mo,  real -viscous  solu¬ 
tions  are  obtained  for  the  same  conditions  as  case 
M2  except  for  the  sphere  radius.  Ihe  radii  for 
cases  M3  and  Mb  are  ten  and  fifty  times,  respectively, 
as  large  as  the  radius  for  M2.  Consequent ly ,  the 
hu ) f- Jurat  ions  of  the  generated  \ -waves  for  M3  and 
Mb  are  about  ten  and  fifty  times  as  long  as  the 
ha  1 f -Jurat  ion  for  M2,  furthermore,  the  distances 
travelled  by  the  wave  fronts  in  cases  M3  and  Mb  are 
about  ten  and  fifty  times  as  long  as  the  distance 
in  case  M2  t o  reach  nearly  the  same  maximum  over¬ 
pressure  . 

In  cases  M3  and  Mb,  the  vibrational  relaxation 
time  for  oxygen  is  much  shorter  than  the  half-dura¬ 
tions  of  the  N -waves,  and  'n'tj  ID"-  for  M5  and 
Jxl O'0  for  Mb  i  ()  -  3 . b  sec),  by  contrast  with  M2, 
where  q  tL|  o.s  at  r  -  19m.  According  to  the 
classification  described  :n  the  last  part  of  Section 
3,  the  former  cases  correspond  to  quasi -equi 1 ibrium 
waves,  while  the  latter  corresponds  to  a  moderate 
nonequ i 1 ibrium  wave.  In  this  section,  the  bulk- 
viscositv  concept  is  introduced  to  evaluate  the 
vibrational  relaxation  for  oxygen  instead  of  solving 
the  relaxation  equation  for  oxygen.  This  assumpt  .  n 
is  reasonable  for  these  cases  due  to  the  fact  that 
the  relaxation  time  or  length  for  oxygen  is  much 
shorter  than  the  character) st ic -flow  time  or  length, 
such  as  the  halt-duration  or  N-wave  length.  In 
practice,  a  typical  time  step  t  for  MS  and  Mb 
becomes  longer  than  q  (  t  -  11.2  ..sec  for  M5  and 
3b. 2  .sec  for  Mb l .  In  the  basic  equations  (Fq. 

(4.1  )|.  the  coefficient  of  viscosity  +  ,.r  was 
replaced  by  .  ♦  ♦  (..vi<),  and  the  method  of  solving 

the  equation*  for  perfect  viscous  flows  was  used. 

That  is,  only  spherical  and  viscous  corrections  were 
carried  out  in  the  operator-splitting  technique. 

Petal  Is  of  the  bulk  viscosity  analysis  are  shown 
1 n  Append  1 x  f  . 

Figure  1.33  shows  the  attenuation  of  (,‘.p)m?x  vs 
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the  normal  ized  dis'unve  r*  t  or  oases  Am  I > J ,  P3  and 
Pm  Case  An  is  a  perfect  -  inv  iso  id  solution  tor  the 
same  1*4  j  as  eases  PJ,  Pi>  and  Po .  In  the  initial 
stage  ot'  |  |'l,!UX  Ion,  the  decay  indices  are  almost 
the  same  for  all  eases  at  higher  !  p),.,.ix  hot  vary  as 
the  waves  weaken.  At  r*  -  J,000,  n  -  1  for  ease  Am 
1.1J  for  ease  PJ,  1.13  for  ease  P5  and  1  .On  for 
ease  Pm  The  decay  index  n  decreases  at  a  fixed 
distance  as  the  ha  1 f -dura t ion  tj  increases.  That 
is,  the  effect  of  vibrational  nonequ 1 1 ibr nun  on  the 
decay  rate  of  i  pim;lx  is  weakened  as  the  wave  ap¬ 
proaches  equilibrium.  Since  n  -  l  far  a  weak  frozen 
wave  (case  AM,  the  maximum  value  of  n  would  exist 
for  a  moderate  nonequ i 1 ibrium  wave.  Among  the  cases 
shown  in  this  section,  the  maximum  value  of  n  was 
obtained  for  case  PJ. 

figure  -1.34  shows  the  normalized  half  duration 
tj  as  a  function  of  r* ,  where  t j  is  defined  by  tj  ; 
ajtj  r,).  Ihe  ha  1  f  -durat  ions  rapidly  increase  in  the 
initial  stage  for  all  cases  due  to  the  nonlinear 
effect,  but  are  quite  different  at  the  later  stage 
of  weakened  waves  depending  on  the  degree  of  v ibra 
tii'iial  nonequ  1 1  ibnutn.  lor  weak  haves  below  l»Ui  Pa, 
t tie  rate  of  increase  is  zero  for  a  frozen  wave 
case  Ao  ;  ,  but  positive  for  nonequ  1 1  1  b  r  mm  waves 
leases  PJ,  P:>,  Po »  .  Ihe  raximun.  rate  of  increase 
of  tj  vs  r*  is  obtained  for  the  .vioJerute  noncquili 
brium  wave  ;case  PJ 1  |as  well  as  for  the  decay  rate 
1  I’  111’..  1.3.’ |  . 

figure  J .  .'o  shows  the  rise  time  tr  as  .1  tunc 
tion  of  t  p'nt.ix-  The  broken  line  indicates  the 
modified  Taylor  solution  for  steady,  roa 1 -v iscous, 
plane  waves,  in  which  the  coefficient  of  viscositv 
♦  .  r  ♦  l  •  x  1^  is  used  instead  of  .  in  the  lav  lot- 
solution.  Ihe  chain  lines  indicate  the  modified 
lighthill  solutions  tor  real  viscous  \-waves  with 
tj  =  ~P  ..sec,  ;>"P  ..sec  and  J,lHU>  ..sec,  respect  1  ve  1  \  , 
which  correspond  to  the  half  durations  at  r*  -  J,otUi 
in  cases  PJ,  PS  and  Po.  It  can  be  seen  m  !  ig.  1.3a 
that  the  rise  time  increases  with  increasing  tj  at  a 
fixed  I’p'r.uix  and  approaches  the  modified  lavl«>r  or 
lighthill  value.  In  case  Po,  the  rise  time  tr  over¬ 
shoots  the  modified  lay  lor  value  for  tin-  higher 
1  '  P  *  mil  x  •  ^  1  s  overshoot  would  correspond  to  the 

overshoot  of  tr  above  the  la>  lor  value  for  an  in 
pulsive  step  wave  described  m  Section  1.3.J,  and 
can  be  improved  by  using  a  finer  mesh  size.  As 
described  for  plane  waves  in  Section  1.3,  the 
present  method  of  calculation  gives  good  results  for 
the  transient  behaviour  of  tr,  but  lias  the  defect 
that  there  appear  some  overshoots  of  tr  above  the 
Taylor  value  for  quasi -steady  waves.  Some  improve 
men t  will  be  required  in  the  calculations  in  a 
future  study.  The  nonstat iona ry  effect  clearlv 
appears  for  1  pimax  TP  fa  in  case  Pa  and  for 
l  p)m,ix  IP  Pa  for  c.ise  PM  so  that  the  tr  values 
tend  to  freeze  and  have  lower  values  than  the  modi 
fled  lay  lor  and  lighthill  values. 

figure  1.3o  shews  the  rise  time  tr  as  a  fane 
turn  of  1  pim;lx  for  comparison  between  eases  Po 
and  PoA.  In  ease  Po\,  the  relaxation  equation  for 
oxygen  was  solved  without  using  the  hulk  viscositv 
concent  for  the  same  parameters  as  case  PM  Ihe* 
black  and  white  circles  correspond  to  cases  Po  and 
PoA,  respectively.  Hu*  broken  line  indicates  the 
modified  laylor  solution.  Ihe  rise  time  values  of 
I >b A  result  in  a  higher  overshoot  of  the  modified 
Taylor  values  than  those  of  PM  Ibis  is  the  reason 
why  the  hulk  viscosity  concept  was  used  in  the 
analvsis  of  \  waves  of  lorn:  durations. 


figure  1.3'  slnuss  the-  normalized  o .  ej  pressure 
pro! lies  t ’p/pi 1  plotted  against  tlu-  normalized 
time  it  -  t  tqi  tor  the  same  maximum  overpressure 
tor  cases  PJ,  !'.>  and  PM  As  the  wave  approaches 
equilibrium,  the  peek  'rise  timei  becomes  sharp  and 
the  hack  bectus.es  straight. 

l-l.o  !  Meets  o!  '■  1  t  ]v  tii  Vibrational  Relaxation 

If  ■  Meets  0}  .ibratjiMj.il  rela'  iiiiM.  for 
nitrogen  01  Mu  it  t  •  tuj -t  1  .mi  of  p.,..<JV,  the  half 
duration  t  j(  ai  :  tK-  rise  time  t  r  an  1  nv  est  1  gat  ed 
b>  introducing  the  vibrational  relaxation  equation 
tor  nitrogen  into  the  governing  equations  in  J  in 
t!u  previous  section.  Ihe  vibtuTiiual  relaxation 
tor  o.wgen  is  taken  into  account  through  the  bulk- 
’.  isoosity  concept  .  Ihe  real  v  1 :  eoijs  ca  1  cu  1  a  t  1  on 
was  carried  out  h>  introducing  the  real  gas  correc 
t  ion  for  nitrogen  in  the  operator  splitting  Muhniqm 
The  details  can  be  seen  m  Appendix  i 

1  he  lower  1  ti  1 1  1  1 1  pres.uie  r.iti>-  of  Tj;  l.-'s 
wis  .  hi'Scn  tor  v  a  -« ■  •  P  uni  Pv  .  •pna  ,  t  u.o 
j'e.  to!  tll.it  tin  ot!'.-,  t  ot  1  1 ' '  I  .» t  1  *  -Jia  i  It  1  .1  Xa  t  I  oi; 
for  nitrog.en  would  appeal  for  tin  lowii  z.a  v  ;  :.:uz. 
overpressure.  •.  .is,.  !'"  is  ,  •-«-.»  1  v  i  ■k.>.h  •.  In:  i. 
tor  the  same  condition-  a*  I1'  ex*r|t  i  or  Ijt 

ill  po  ,  on  1  \  tlu  v  ihjat  mr.al  rvlacg  i.'ii  :  V-  r  ox;,  gen 
m  Taken  into  account  b»  .‘i-mg  tlu  r  »= j  ]  >-  v  iscoitv 
concept  ,.  3  .  1  ic  .  a  -a  Ps  i !  leal  .  i  s<  ous 

so lut  ion  f«T  the  s.i-i-  cond  1 1  i .  n -  . i >  c.ise  P"  iiuludin 
th.e  vibr.ition.il  re  !  i  xa  t  i  on  foj  nitrogen  p  1  ...1 
msec  •  .  In  1 1  l  a  1  I  >  ,  it  wa  Jc  •  i  i  «:•  1  v  t <  >  c  » i  m  .  : ;  t  t  he 
calculation  for  tlu  el  feet  ot  nitrogen  h>  adding  the 
vibrational  e.xv.  1 1  a  *  i  er  foi  nitr.gcii  to  va>«-  Po. 
,ii>w«ver.  it  w .  i  *»  ;  etiiid  tfi.it  :  ook  a  lone  c  •jata*i  i 
t  iii.e  for  the  wave  to  i  in  t!u  -;i.i  ^  :  :;,u;z  o.  erji'cs  sur-.  , 
enouef.  to  show  tlu  efttu*  >f  tutispgeti.  lnstt.i.l,  tile 
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j  igure  l  11  shows  .i  somp.ii  i  on  between  the 
pji'».uri*  protilt-.  at  r*  l.P.”'  t’vr  cases  P“  and. 

Ps,  .aid  I  ig.  1  1  '  -d-.ows  piot' die-  ot  tenjuratuiu 

and  vibiational  t  enpc  r.i  1  un  t\:i  mtnvo',  at  the 
same  1  ovation  !>t  !,v  ihese  tigures  indicate  that 

the  wav<*  of  P  ‘  at  r*  l,P3u  corresptMids  Ic  a  higlilv 

nonequ  i  1  i  b  r  l  nr.  wave  ter  vibrational  excitation  of 
nit  logon,  and  thus  tlu  -heck  transition  l s  mainly 
controlled  tu  the  vibiational  ivlaxatu'-i  v>t  v'xvgen. 
Ilu-  in  t  roger.  nonntu  i  I  i  hr  ium  aits  onI\  to  lower  t  hv' 
max  limn;  overpri'ssm  e .  It  r.a;  be  said  that  the 
dissipatin'  etfect  of  vibrational  nonequ  i  1  l  hr  i  uni 
for  nitrogen  m  case  PS  pl.ivs  tole  like  ovvgen  m 


1 


the  increase  of  shock  thickness  cannot  follow  the 
change  in  maximum  overpressure  tnonstat ionary 
effect),  furthermore,  a  longer  relaxation  time 
results  in  a  shorter  rise  time  in  contrast  to  a 
steady  wave,  since  the  shock -th  ickening  time  is 
nearly  proportional  to  the  relaxation  time  for 
dispersed  waves  (nonstationary  effect). 


i  ~  i  As  tile  duration  increases,  the  rise  time 
approaches  the  modified  Taylor  value  for  steady 
plane  waxes  or  the  modified  1 ighthill  value  for 
quas 1 -stat ionary  N-waves,  which  is  obtained  by 
introducing  the  bulk  viscosity  concept  into  the 
viscous-flow  analysis.  lor  waves  with  longer  dura¬ 
tions,  the  nonstationary  effect  on  rise  time 
appears  only  for  lower  maximum  overpressures. 


case  Ml,  while  the  vibrational  relaxation  of  oxygen 
acts  as  a  bulk  viscosity  in  case  08.  figure  4.4  2 
shows  that  the  wave  is  a  partly-dispersed  wave  for 
nitrogen,  fully-dispersed  waves  for  nitrogen  could 
be  obtained  for  waves  with  longer  duration  and  lower 
mu  x  i  mum  o ve  rp  res  su  re  . 


5.  coxa  US  ION'S 

Ihe  foregoing  results  can  he  summarized  as 
f o 1 1 ow  s  : 

( I )  It  was  shown  that  the  transient  shock  structures 
of  weak  plane  and  spherical  waves  in  air  can  be 
analysed  by  solving  the  unsteady,  compressible 
Nav  ier-Stokes  equations  with  a  vibrat ional -relaxa • 
t ion  equation  for  oxygen  or  nitrogen,  using  the 
random-choice  method  (RCM)  with  an  operator-split¬ 
ting  technique. 

id)  The  perfect -viscous  and  re a  1  - inv isc id  solutions 
for  impulsive  step-waves  show  that  the  smearing 
processes  due  to  dissipative  effects  of  viscosity 
and  vibrational  nonequilibrium  for  shock  fronts  are 
m  reasonable  agreement  with  analysis.  However, 
there  is  some  randomness  in  the  shock  thickness  or 
the  rise  time  value  and  there  are  some  overshoots 
above  the  steady-state  value. 

i  Ihe  initial  X-wave  formation  process  was  estab¬ 
lished  for  a  perfect- inviscid  wave  for  exploding 
pressurized  air  spheres.  It  was  found  that  the 
attenuation  of  the  maximum  (peak)  overpressure  and 
the  ha  1 f -dura t ion  of  an  X-wave  in  the  far  field 
can  be  controlled  by  a  proper  choice  of  sphere 
radius  and  initial  diaphragm-pressure  ratio. 

ill  The  perfect- inviscid,  perfect -vi scous ,  rea  1  - 
inviscid  and  real -viscous  far-field  solutions  for 
weak  spherical  waves  in  air  were  compared.  It  was 
found  that  the  dissipative  effects  of  viscosity  and 
vibrational  nonequilibrium  of  oxygen  on  the  decay 
of  the  maximum  overpressure,  half-duration  and 
X-wave  rise  time  become  distinguishable  for  values 
of  i  p ) ma x  1M  l’a  . 

(3)  Ihe  numerical  simulations  were  carried  out  for 
weak  spherical  X-waves  generated  in  atmospheric  air 
from  sparks  and  exploding  wires.  The  numerical 
results  show  good  agreement  with  the  experimental 
data  with  regard  to  the  decay  rate  of  (.‘p)max»  the 
increasing  rate  of  tj,  the  rise  time  tr  and  the 
wave  profiles.  The  results  indicate  that  the 
observed  shock  structures  of  weak  spherical  X-waves 
are  controlled  by  the  coupled  dissipative  effects 
of  viscosity  and  vibrational  nonequilibrium  of 
oxygen . 

l(>)  The  calculated  and  observed  rise  times  (or 
shock  thicknesses)  for  weak  spherical  X-waves  are 
mostly  much  smaller  than  those  predicted  analytic¬ 
ally  for  steady  plane  waves.  It  is  found  that 
this  phenomenon  is  attributed  to  the  and  the 

*  «  "■'»  coupled  with  vibrational 

relaxation  of  oxygen.  The  shorter  half-durations 
give  shorter  rise  times  for  the  same  maximum  over¬ 
pressures  due  to  flow  expansion  behind  the  front 
shock  of  an  N  wave  (X-wave  effect).  A  more  rapid 
decrease  of  the  maximum  overpressure  also  results 
m  a  shorter  rise  time  for  the  same  maximum  over¬ 
pressure,  since  the  shock-thickening  time  becomes 
i nc rea singly  long  as  a  wa ve  is  wea k en cd ,  so  t ha t 


(8)  Ihe  deca\  index  n,  which  denotes  the  local  decay 
rate  of  i  plmax.  defined  by  f'pjmax  r~n,  is  equal 
to  unity  for  a  classical,  linear  acoustic  wave,  but 
increases  due  to  the  dissipative  effects  of  viscosity 
and  vibrational  noneqm 1 ibrium  for  moderate  non- 
equ i 1 ibrium,  weak  spherical  X-waves.  It  approaches 
unity  tor  quas i -equ i 1 ibrium  waves  of  long  duration. 


finally,  as  noted  in  Section  !,  the  original 
motivation  tor  the  present  study  was  to  answer  the 
question  whether  X-waves  generated  by  sparks  or 
exploding  wires  can  simulate  SS I -generated  N-waves 
including  real  gas  effects  on  X-wave  rise  times.  Ihe 
answer  is  'no',  since  unlike  SSI -generated  N-waves 
produced  hv  sparks  and  exploding  wires,  the  rise 
times  are  determined  hv  \-wu\e  duration  and  the  non- 
stationary  distance  travelled  from  the  source. 
Nevertheless,  the  present  study  is  important  since 
it  has  succeeded  in  providing  appropriate  explana¬ 
tions  for  the  rise  times  of  spark  and  exploding-wi re 
generated  N-waves  by  using  the  concepts  of  the 

*•  and  the  •  •  .*  ’  *  **  with  the  aid  of 

very  good  Rt'M  simulations  of  tlu*  actual  experiments. 
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V ibrat  iona  1  Relaxat  ion  T i me s  t  or  Oxygen  (  *  )  and  \  1 1  rogen  (  ^ ) 
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I  1 -room  temperature,  RH- relative  humidity. 
All  absolute  humidity,  l.W - expl od ing  wire 


Parameters  t'ur  Computation  ot  Spherical  have? 


(at  Perfect  -  In  vise  id  Mows 
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(a)  Spark 


1  SPARK  AND  iiXPLODINC',-WIRl:.-Gi;.Ni:RATnD  N-WAVLS. 
(a'  STRIP'S  I  -  SPARK 

S  =  b.O  KV,  r  =  21. 6nw  (Ap)max  =  8.2  Pa, 
tj  =  72  us,  tr  =  11.9  us . 

(b)  STRIPS  IV  -  EXPLODING  WIRIi 

S  =  b . 0  KV,  r  =  29.3m;  (Ap)max  =  20.2  Pa 
td  =  122  us,  tr  =  15.2  us. 


» 
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FIG.  2.3  HALF-DURATION  DATA  trf  vs  r. 


SPARK:  SERIES  I  -  O  AND  II  -  A  ; 

EXPLODING  WIRES:  SERIES  III  -  O  AND  IV  -  D  . 


I0'5  2  5  10“'  2  5  I0‘3 

(Ap^/p, 


i-IG.  3.2  NORMA!.  I  Zi;i)  TAYLOR  T1IICKN1.SS  ( *.x )  «  /  (  7;i  j  )  OR  NORMAL I  ILL 
TAYLOR  R I  SI:  T I  Ml:  'r()/(  i?a  i  - )  P'  C  r.H  AGAINST  SHOCK 
STRI-NCiTl  I  (\p),/p 


<AP>max  <Po) 


COMPARISON  Bi:  IVVI-r.N  l.XPI.RIMI  MAI.  AND  IllliiRI  ill  \! 
RISC  TIMIiS  PI  OTTHD  AC.AINST  MAXIMUM  OVI  KPK!  ^HRI 


SPARK:  SCRIPS  I  -  O  AND  II  -  A  ; 
CXPLODINO  KIRC:  SCRIPS  III  O  AND  I\ 
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I  I  lit  ITIII  1,1,  N-KAVF.S  (ASYMPTOTIC  FORMS  OF'  PUFSFS  WITH  ZFRO 
MASS  FLOW)  FOR  RFYNOI.DS  NUMRTRS  Re  =  30,  10,  3,  1,  0.3 


I.  Ii;m  III  I.[  N  —  IV  A  VIS  IN  A  I’i.OT  Of  Till  NORMA  II  II:P  OVliRPRIJSSURT 
I  p  )  /  I  p  )  max  vs  1111  I'ISI  ANCi  I’ARAMI  1 1  R  I  I  OR  DURATION 
I’ARAMI  HRS’ Zj  -  ■,  d() .  4  ,  JO.  9,  3..H7,  1.93. 


Taylor 


l0-5  2  5  |0-4  2  5  |0- 

(^P^max^Pi 


FIG.  3.7  NORMALIZED  SHOCK  TliICKNLSS  AX/(  3/aj)  OR  NORMALIZED  RISE  TIME 
tr/(6/ai2)  VS  NORMALIZED  MAXIMUM  OVERPRESSURE  (Ap)max/pj  FOR 
NORMALIZED  DURATION  tj/(6/ai2)  =  5x10$,  2x10$,  105,  5xl04 
FOR  LIGHTHILL  N-WAVES. 


AP  (Pa)  “  AP  (Pa) 


SERIES  I  -  SPARK  DATA 

.IGHTHILL:  (Ap)max  =  8 . 52  Pa ,  tr  =  1 2  us 

SERIES  I:  r  =  21.6m,  (Ap)max  =  8-52  Pa > 


tj  =  72  us. 

tr  =  11.4  us,  td  -  72  us. 


(b)  SERIES  II  -  SPARK  DATA 

I. IGHTHILL:  (Ap)max  =  5.83  Pa,  tr  =  16.7  us,  td  =  76.8  us. 

SERIES  II  r  =  19.0m,  (Ap)max  =  5.83  Pa,  tr  =  15.5  us,  td  =  76.8  us. 


FIG.  3.10  COMPARISON  BETWEEN  EXPERIMENTAL  -  AND  THEORETICAL 

(LIGHTHILL)  PRESSURE  PROFILES  OF  N-WAVES. 


AP  (Pa) 


pressure 


temperature 


(a)  Fully  Dispersed  Wave 


I(i .  3.14  PRHSSURC  p,  AND  TFMPFRATURF.  (TRANSLATION  +  ROTATION  T, 
AND  VIBRATION  1'v)  PROFI  I.FS  OF  SHOCK  TRANSITIONS  WITH 
VIBRATIONAL  I.XCITATION  FOR  FULLY  AND  PARTLY  DISPFRSHD 
WAVFS . 
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01-  INITIAL  TI  MI’LRATURi:  Tj  I  OR  AIR. 
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1  I'lj  I  P  I  c  r ,  i  l()K  Bill  I  Y  AND  I’ARTIA -l>  I  SI'IIRSID  KAY!  S. 
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COMPARISON  H!  i'W i  I  X  I  XI'I  RIMI  NTAI.  |  SCR  I  l;.S  IV)  ,  IiW  ■  AND  TllliORP.TICAI. 
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(c)  Series-I[I:  Exploding- wire  data 

(c)  smri i:s  ni:  i:\ri.om.Nt.  mm  i •  m  \ 

Monirn n  i.ioimn  1 1 ;  i ' p)nuix  :  1_  |,;1-  t,.  u.s  ..s.  td  =  hs.o  ..s 
Si:RIIS  III:  r  -  27  .cm,  I  '  p  l„u,x  =  1"  I'a,  t  ,•  in.;.  .  s,  td  -  1 1 3 .  o  ..s 


t(^sec) 


(d)  Senes-JZ:  Exploding -wire  data 


(d)  SCRIPS  IV:  r.XPl.ODI  NC-IVI RP  DATA 

MOnilir.n  I.IC.IIIIIILI,:  ('P)itiux  =  15.3  Pa,  tr  =  37.4  us,  tj  =  105.3  us. 
SPKIHS  IV:  r  =  29.3m,  (.'.p)max  =  17  Pa,  tr  =  18.7  us,  tj  =  105.3  us. 


10.  3.19  -  CONTINUPI)  -  COMPARISON  BPTWPPN  CXPCR IMCNTAI.  AND  THCORCTICAP 

(MOD  mill)  I.IGHTI11 1. 1,)  PRPSSURP  PROFIPPS  OP  N-WAVPS. 
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(a)  Quasi-Equilibrium  Wave 


(b)  Moderately-Nonequilibrium  Wave 


I  (e)  Quasi-Frozen  Wave 


0  CLASSIFICATION  OF  N-WAVLS  ACCORDING  TO  Till-;  DFCRLF.  OF 
NONF.QU I  LIBRIUM. 
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R(i'.  4.4  SHOCK-TUBF  PROBI.FM  USING  RANDOM-CIIOICF  MliTHOf)  FOR  A 

PFRFF.CT-VI  SCOWS  FLOW  (P41  =  2.0,  T4  j  =  I  .0,  Ax*  =  1/40, 
X()  =  0.001  cm). 
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TRANS  II. Yl  BIIIAV10IIK  01  SHOCK  lllIiKMSs  I  OR  SIIOlK-IHBI  I’ROBI.I.M 

ns i m;  rcm  ior  a  itriici  visions  inm  irn  -’.o,  l.j]  =  >■". 

1  i  -  K ,  p|  '  1 01  .  A  Kl’a  I  . 
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(c)  COMPARISON  BETWEEN  ROM  AND  MF'M  (MacCORMACK ' S  FINITE-DI FFLRI NCI 
METHOD) 

+  X  =  1.25.\l()-5  cm  (MacCORMACK) 

C  'X  =  1.25x10-5  cm,  a  ,\X  =  0.75xl0-5  cm  (RANDOM-CHOICE) 


F10 .  4.6  -  CONTINUED  -  TRANSIENT  BEHAVIOUR  OF  SHOCK  THICKNESS  FOR  S1IOCK-TUBI 

PROBLEM  US  I  NO  RCM  FOR  A  PERFECT-VISCOUS  FI  .OK 
d’41  =  -•«,  T41  =  1.0,  T,  =  275  K,  p!  =  101.3  KPa). 
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I  Id.  4.7  TRANSIENT  SHOCK  PRESSURE  PROFILES  FOR  SHOCK-TUBE  PROBLEM 

MS  INC  RCM  FOR  A  PERFECT-VISCOUS  FLOW.  (Pm  =  2.0,  T4 \  =  1.0, 
Tj  =  2-3  K,  pi  =  101.3  KPa,  x  =  1.23x10-5  cm). 

- 1.1  dllTI  1 1 1. 1. ,  •  i  =  0.99,  a  i  =  45.0,  o  r  =  58.3. 


FIG. 


4.9  IRAN'S  I  F.N  i'  BF.MAY  I  OUR  01'  SHOCK  TIIICKNFSS  FOR  S110CK-TUBL  1’ROBl  FM 
USING  RCM  FOR  A  UFA). -I  NY  ISC  II)  FLOW  j  P4  j  =  1.0018,  I4  [  =  1.0, 
Pi  =  101.3  Kl’a ,  Tj  =  303.15  K,  RH  =  90", ,  Me  =  1.0004,  Cp)->’  = 
91.1  Pa | .  “ 

A  'X  =  0.075  cm,  O  X  =  0.0 US  cm. 
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FIG.  4.10 


STF.ADY  SHOCK- PRFSSURF  PRO!  I  I  F 
USING  RCM  FOR  A  RF.AI.-  INY1SCI0 
1 4 1  =  1-0,  pi  =  101.3  Kl’a ,  T] 
Me  =  1 -0004,  (  ' p ) j  =  91.1  Pa, 

0  =  1-04  sec,  xq  =  0.5  cm,  Ax 


FOR  SHOGK-TUYP;  PR0B1.FM 
FI. OK  |P4i  =  1.0018, 

=  505.15  K,  Rll  =  90", 
(■FI J  =  0.0777  K, 

=  0.0175  cm] . 


=  77.0 


-  ANALYTICAL,  •  .  =  J5.0,  o  i 
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III..  1.11  COM  I  Mil  l> 

MAR-tini'  snnnioN  oi  ixi’iosiov  01  a  i’ki ssuk i zii>  air  smi'Ri 
iisim;  ram  h  )M  ciioici  Minion  ion  a  muiri  iwisi  in  i-mn  asi 

\l,  I- ,  -'.II.  I  l.u,  r*  i.  im. 

1 1>  ■  1 1  IRMA  MOV  m  \  HAW 


(  0\J  I  XU}  D  -  )  }  H  CIS  Of'  JMTI.AI.  PRISSIIRI:  AM)  TIMPPRATURl  RATIOS 
>4i  AM)  ON  N-W-WI  PRO!  III.  (  ‘  r*  =  1/10). 


CONTI NU!-;i>  -  CM  liCTS  OF  INITIAL  PRHSSURI  AND  THMPFRATIJRi:  RATIOS 
f>4 1  AM)  I41  ON  N-WAVM  PROFU.fi  (V.r*  =  1/10). 

,d)  P41  =  9.0,  I4!  =  9.0  (CASH  A 4 ) 


Sound-wave  Path 


4-13  PAT"  0F  SI,0^F  FRONT  FOR  A  PERFl-CT-INVISCID 


FLOW  (CASE  A5) 
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COMPARISON,  HI.  Ihl.l  \  PPRM  <T- I \V  I  SCI  i>  (AS  -  a),  I'iRWCT -VISCOUS  IH1 
RHAI.-1NVISCID  (Cl  -  +),  AND  RHAl, -VISCOUS  (HI  -•),  I-AR-!  1  f- 1.1)  RCM 
SOI. II 1  1 ONS :  1  IAI.H-DURAT  I  ON  t<j  VS  HI  STANCH.  r. 


COMPARISON  BHTWHH.N  R1SI:  I  1  MIAS  tr  VS  IU  STANCH.  r  FOR  PHRHHCT- I NV1 SC  1 D 
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RPAI. -VISCOUS  (HI  -•),  PAR  I  II  ID  RCM  SOLUTIONS:  RISK  T1MP  tr  VS 
I)  I  STANCH  r. 


(AP)/(AP)mox  IAPM  AP)max  (AP)/(AP)mox  (AP)/(AP)max 


(a)  Perfect -inviscid  (Case  A5) 


( b)  Perfect  -  viscous  (Case  B 1 ) 


(d)  Real-viscous  (Case  D1)  5mm 


HO.  4.17  COMPARISON’  or  f’l  K’fl  CT  I.WLSCID,  1’LR) JIC  J -VI SCOUS,  RFAL- INVISCID 
AND  RI.AI. -VISCOUS  PRLSSURT  PROI  I  I  I  S  AT  SLVLRAl.  LOCATIONS  FOR 
RI1AI, -VISCOUS,  FAR-I 1ILD  ROM  SOLUTIONS. 
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1  .'< . 

Mi'IIKUM:  i  p )  ni.i  x  i’-* .  tr  <>.('5  ..s.  tj  = 

i  M’!  k  I  Ml  \  I  \l. :  'p)max  la. ‘A  I'a  ,  tr  '  .  1 ..s,  tj  =■ 
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\l  r  1  .  I'm 

Mi'll  k  1  CAI  :  <  P )  max  1,;'  >  t  ■  <>.'*)  ..s,  tj  -  SJ." 

I  M'l  k  I  Ml  M  M  :  Cp)ni!ix  S.K.  Pa,  tr  II."  ,.s,  tj  =  .. 


( 'i i'll’ \k i si i\  Hi  M'\kk  phi  ssi ik i  pkOMiis  oi  coMpimp  \-iv-\vi 

Willi  I  \  PI  li  I  'll  \  !  \l  P  A  I  A  :  MIR  1 1  s  |  VS  CASI  PI. 
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i..::,  i  ni  i  \ -wAvi  I’poi  mis  ni  puissum,  n  mitr-viuki  ami 

VIHRAIIiIVM  IIMITHATimi  Al  r  J 1  .  (,i:i  I  OR  CASI  IMA. 


AT,  (ATV)  K  AP  (Pa) 
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COMPARISON'  Of  p  VS  t  IOR  CASPS  DJ  AND  D!  (OR  f  I  f  f  J  f  Kf  \  | 

VI  BRATIOXAI  Rl TAXATION  PIMPS  (D_':  S  ",1  ;  n.j  •  1  , 

..see  IOR  ).  ’ 
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r  =  1  P  'max  =  0.35  P;i ,  tr  =  1  (, .  “H  :.soc  tj  = 
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COMPAR  I  SON  01-  (n)m.lx/pi  VS  r*  IOR  PI  Rif  CT  - 1  XV  1  SC  1 P  CASH  Ac  O 
AM)  RIAL  VISCOUS  CASPS  I)  J  a  ,  DSD  AND  D<>  ■  (  P ,  ,  =  ].s, 
r„  ^  1.15  cm  IOR  D J). 
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4.34  COMPARISON  OH  tj*  VS  r*  H'OR  CASHS  A(>  I  P I  RI- I:C  J  -  I XV I  SC  1 1  > )  , 

D:,  1)3  AM)  Do  ( RI  AL -VISCOUS)  COR  PI  FFF.RHNT  HAI.1-  DURATIONS 
(P41  =  1.8,  r0  =  1.15  cm  TOR  f>2,  iU :  11.5  cm,  1)0:  5'. 5  cm 
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0.  4.35  COMPARISON  OH  tr  VS  (.\p)m  HOR  PIFFHRHNT  HALT  DURATIONS 
FOR  RHAL-VISCOUS  CASHS  D Z,  D5  AND  D6  (P4I  =  1.8,  r..  =  1. 
cm  FOR  1)2.  1)5:  11.5  cm,  06:  5  7.5  cm).  0 


AiTl.M'  I  X  A 


1YAIUAI10N  Ol-  p  ,  d  ;\M »  d  ,  IN  i HI  IH.H1HIII  N  V\ A\  i  Mil'lfiuN 

- - -  m  •  i  x  - -  d  -  -  •  —  —  —  •  -  -  ------  - . 


Ihe  centre  el  the  shock  front  of  .in  \-uau*  can 


e;i  1  cul.it ed  from  the  iterative  equation- 
tier  1  cod  f  n*ui  l  q  .  1  '■.13  ,  as. 


i  J  .  n  •  e  x  p  \  Re  1  - 1 


■  n  exp v  He 


which  I**  denied  from  lq.  i3.15i,  where  .  is  the 
at  T  =  n.;,  in  the  shock  front. 


Ihe  peak  point  of  an  N-wave  can  he  obtained 

from  i  dl'  d;„  .  ;  i'.  Ill  is  gives  the  relation 

He = const 


J  ■  n  -  exp  «  He 


exp » He  i  1  -  i‘  “  - 1  !  ex p  i  V- d  1 


where  k  is  the  u umber  of  iteration 
values  of  iteration  used  were  '  c/ 


where  'm  is  the  at  I’  -  l’max  .  l.quat  ion  :A.Ji  is 
rewritten  as 


,  ex pi  He)  -l 
=  _  •  n  *  — - - 


to  evaluate  n;.  Ihe  method  of  successive  iteration 
is  used  tv'  solve  !  q .  i.A.oi  for  tn- .  Hie  centre 
value  was  used  for  the  initial  value  of  item 
tton.  Ihe  value  of  is  obtained  by  substituting 

•  '  m  into  1  q .  i  s  .  1 5  i  . 


Ihe  shock  thickne* 
’ j  are  obtained  from 


and  t he  ha  1 f -dura t i on 


o’  d  '  1 


where  j  and  <,  are  the  at  i  -  0.1  Pm.lv  and  o.O 
Pjiuix .  respectively.  Ihe  values  of  j  and  \>  are 


ihe  value*  of  d  and  dq  are  uhtui 
s  t  1 1  u  t  i  n  c  r,lutx,  and  q  so  obtained 
fa.il;  and  ,3 . 1 3  •  .  Die  computer  piorr 
l n p  i'jna \  *  -  and  dq  is  given  m  !  i  stin 

plot  of  dq  vs  He  is  shown  in  lip.  A.l. 

In  order  to  compare  the  expv r  i"en 
with  the  lighthlll  N  w.i  v  e  sol  ut  i  on  s  ,  : 
to  establish  the  procedure  tor  ohtaini 
spend init  lighthlll  rise  time  fiom  the 
values  of  the  maximir:  .  a  rr;  l'rssun  p 
observed  ha  1  f -dural  i on  t  j.  \t  t  :  r - 1  . 
mined  t’rom  lq.  3  .  1  j  .  .  then  thi  qrn 
of  He  is  read  out  from  the  dq-Kc  in 
|<y  repeating  the  process  of  lq".  A.l 
approximate  value  of  dj  is  found.  ilu 
process  is  repeated  when  dq  so  obtairu 
coincide  with  the  exact  dq .  finally, 
spending  value  of  He  is  determined  aru 
obtained.  Ihe  actual  shock  thickness 
t’rom  fq.  i3.!4i.  rile  computer  pr.vnr 
procedure  i s  given  in  I  isting  A.d. 
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C2n'  =  “’50RT(lTAv1  ) 

0C52C 

p  •iax=gZ’’,/(  i  .o*:e  *•' 

(  f  TAV1  /?.  ''  )  / p 4 9  ) 

00330 

E 4=0.9BF“AX 

00340 

El  =  C  .  1  aP“AX 

00350 

0  29= G 1  * 

00360 

■  1 

6Z9A=rSiRT ( 2.0aDL~ 

•;,(  ERRBI  GZ9/P9-1  .  C  )  )  ) 

00370 

IE  (  "  A  R  5  (  G  2  9  A  -  G  2  9  )  . 

LT.l. 0P-06I  GO  TO  30 

0  0  3  c  0 

C-  2  9  =  5  Z  9  A 

0  0390 

G  j  T  t’  31 

0  0400 

3  0 

C21=?.0«G2C-S29A 

0  0  4  10 

41 

G21A  =  n,Si.RT(?,0aDL" 

G (ERRaiG 21/ R 1-1.0)  1  ) 

0  0  4  2  0 

IP (DA«S(6Z1A.GZ1 ) • 

LT.1.0D-0«>  GO  TO  40 

00430 

0  2  1  =  '•  Z  1 4 

0  0  4-0 

GO  T"  41 

0  0450 

4  0 

'  G/  =  GZ I  A -GZ 9  A 

00460 

C  2 '  =  G  Z  1  A 

0047  1 

[  Z  =  G  ”  D  G  Z  a  P  M  A  x 

0  0  4  -  c 

ZO  =  GBGZr  aP"Ax 

0  0  4  4  7 

..-1  I  TP  (  6»50  )  R9,EZ» 

2" 

G  0  50  C 

7  0 

c:\tim.e 

C051C 

10 

f  j  a  m  a  t  c  1 1  n  i 

00520 

1? 

r jkmaT( 7F 1 n.2 ) 

C  0  5  3  0 
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EDR“AT(//1h  ,  10X»’ 

aaa  UGHTkILL  ANAlITICAL  SOLUTION 

00540  aSA^CTfR  AND  THICKNESS  PARAMETER  •  ,  /  1  9  ,  1  X » 

0055 C  a  '"P  N-..AVE  E R " GIVEN  REYNOLDS  NUMBER  aaa’, 

C  0 5 1 0  a  ///in  »22x, ’R1 ,11X1 ’PZ 1 .1 3X, • IC / > 

00570  5  C  EM-'-ATtlH  *10X»3D15.7) 

005^0  STEP 

0  0  5  9  0  EM 

f NC  CF  2atA 
tNC  S 

SAVE  0  TO  "AT  A  SET  (  '  TN3D0  00  .  L-<  l  .  POST  •  ) 

READY 
TSLOG  end 

a  a  a  a  a  a  a  a  a  a  «  a  a !!  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  it  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a  a 
a  USE0  =  TN30000  a 

»  p«r^„,,uRE  =  L  °60N2  a 

a  TSIDG  ENDED  T I “E  =  1 1 : 53 !  12  DAT E =62- 1 2-0 8  a 

aaaaaaaBaaaaaBaaitaKaaaaaaaaaaaaaaaaaaaaaBaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


AIM’ I  NI)  IX  B 


IM  RIVA'l  ION  OI;  ANAI  YTluA!  HI.l.A’I  IONS  IN  SliCUON  S.  i 


.1  i  'on  v  at  i  on  o  f  lq .  ( '  ri  t  i  c al  Ov  e r 

pressure 

donsider  a  steady  normal  shock  wave  in  a  ga< 
with  vibrational  excitation  ami  assume  the  equili¬ 
brium  states  tor  both  sides  of  the  shod'  front. 

hen,  the  equations  of  continuity,  momentum  and 
energy  are  given  by 


Pi 


i  H  .  I  i 


t  b  . :  s 


<  R  .  a  • 


wiieia-  file  subscripts  1  and  J  denote  the  states 
ahead  of  and  behind  the  shock  front,  respectively; 
.  ,  density;  u,  velocity;  p,  pressure;  1,  tempera¬ 
ture;  ,  vibrational  energy,  1  ^ ,  specific  heat  at 
constant  nre-'sure  for  t  ran- lat  i  onal  and  rotational 
energy,  a  sumed  constant.  ihc  equat ion  of  state 
and  the  vibrational  energy  arc  assumed  to  be  ex 
pressed  by 


■  i- .  \  » 


,tt  t  v >  !' 
•'.ltd.  . 


fie  g.i  *  constant;  c,,  vibrational  spec  i  I 
;  molecule,  assumed  constant  across 


B.l  .,nd  •  i. .  J 


p  .  . 


1U1  1 


l  l 


'  l 


“l 


on  ’’.ich  nuroho r  .  >!<  f :  nod  by 


Substituting  l.q  .  tB.6)  into  l.q.  ,  B  .  "  i  obtanv- 

.  r 


$ •  v 


1  u 


:  V  1 


i’  . 

R 


1 


'/ 


1  1 

”IIY 


,  - 1 1 


V  •  si 

Substituting  l.q.  i  B  .  d  j  into  i  q.  .  .  t  >  ,  ,  (hf ; 


‘•P*>  P-i-P, 


1 


/[•♦I-  ■ 


Putting  1,  the  expression  for  the  crii; 

pressure  i'Pvr.i  *  * 


1  i 

B.J  derivation  ot’  I  .•'.SI.:  I > l  f f;;.'  .  •.  :  t  > 

i qua t ion  S.J(>  for  •  can  bt  fwMttv 


-  *  J 


■  b  . ;  i 


«l  it;  •  ;  the  m!  io  of  specific  heat  in  vieritiOii 

a  ;  squ  i  1  i  i.r  i  m;i;  VA  ,  -j'octl'ic  heal  if  constant  \oluru 
tor  translational  and  rotational  eiierg>. 


r  :ng  the  re lat i 


on  turn  t rom 


.  R  .  i  .' 


uid  •  B  .  1 , 


n.  1  it  ;  on  .  B  . 
•  •  w  r  i  r  t  in  .i 


Hid  i  B  .  (>  ;  ,  I  qs  .  :  |:  .  s 


Per  i  vat  •  on  ■'! 


I  ..  :um»  :  e*  k 


ihc  parameter  k  appearing  tn  lq.  .  .  J  l  ■  cat; 
row  r  i  M  on  a  • 
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MICROCOPY  RLSOLUTION  USI  CHARI 

N.Vii'NAl  HllKI  All  Of  s’ ANl’AKOS  1°*'  *  A 


(6. 13) 


ae  al  '  2  al  a,  ■  2yM 


(Ap) 2  1  (AP)- 


introducing  the  relation  (3.29),  then 


ae  al  2  al  Qj  2yMf  Pi  2Y  Pi 
Substituting  Eq.  (B.13)  into  Eq.  (B . 12),  obtain 
1  .  yl  ^P)2  Y  ,  (AP)2  Pi  _  (Ap 


(*P)2  Pj  (Ap)2 


k  2>  Pl 


( v - 1 )  Cj  Pl  *-AP^cr,j  ^Ap^cr,j 


Y+1  f2(Y"1)2cj 

(Ap)2-| 

2y  L  Y+1 

pl  -1 

(Ap)cr,j  r  . 

fAp)2  ' 

pl  1- 

(Ap) 

K'cr,j  ■ 

B.4  Derivation  of  Eq.  (3.34)  from  Eq .  (3.24): 
Polyakova  et  al  (Ref.  21) 

Equation  (3.24)  can  be  rewritten  as 


(3.32)  The  terms  in  Eq.  (3.37)  can  now  be  rewritten  as 

f  2  2  G2l  G2  y+1  tAp)cr  if  (Ap 

1+Ym/-(y+1)m/  I  -±  5  1  ♦  Txzrr 


P  .  fi  -  *  j  tn  (l  + 

kt.  {  k  )  ( 


f  -  hi  h 

\i  >  *•  >  K  V 

0  0 


Using  the  notations  for  a  normal  shock  in  B.l 


Vu  =  v+vQ 


v  Vu 
1  -  —  =  2  — -  =  2 


P1  _  2  (Ap) 

^2  Pi  (AP)2 


Using  the  expression  (3.32)  for  k 


y  _  Y+1  ^AP^2  Y  y  Y+1  al  y  ^Ap^2_  Y+1 

Pl"  (Y-l)2Cj  V  2Y  T^T  pl 


where  Z  is  defined  as 


aj  y  ( Ap ) 2 

(«v) .  "p^ 


Substituting  Eqs.  (3.33),  (B.16)  and  (3. 35)  into  Eq . 
(B. 14)  obtains 

IT  “-V-F-tS— ]'"[■-%] 


-  -  IT  ~  1  2n 

L  (Ap:icr,j  L  (AP)2 


zo  ’  $T  [n-  *  !  H  '  const  ( 

B.S  Derivation  of  Eq.  (3.34)  from  Eq.  (3.37): 
Johannesen  and  Hodgson  (Ref.  12) 

From  Eq.  (B.10) 

2  [Y*l+2(Y-l)c.] (Ap),/p  +  2[Y*(Y-1)c.] 

Mf  * - 2y  [  1  +  (Y~i~)~c  j  1  - L  < 


-(Y+1  )Mf 2  —  fl  -  ~ 


(Ap)cr.jf1 

(Ap)  - 

pi  L 

+  (Ap)  .j 

r  Cr,JJ 

Y  +  l  (Ap)2 

(B. 20) 

f  R  1 

Y  Pl 

»  . LI  ) 

Mf  [(y*1).2(y-1)c.]  L  u2  K  (Y*1)X  (AP)2 


2yVj  Pl 


Thus,  Eq.  (3.37)  can  be  rewritten  as 


(3.15) 

(Y+l)x  (Ap)2 

Y+1  ^Ap^  2 

(3.33) 

2Ya i T  j  Pj 

Y  Pj 

Y+1  ^Ap^cr,j  1 

1  .  (Ap)2  1  *n 

r  i  (Ap)z 

III  i 

2y 

Pj 

r+l  (Ap)cr,) 

.  W'.r.J  J 

r  lflp)2  l 

_ymJ  pi 
r  i  (Ap>2 

(B.16) 

2Y  Pl 

1  TKp)  1 

1  cr.j  J 

Um2  pl 

(3.35) 

1 

2y_  _x_  ____2__ 

Y+1  at  (Ap) 

1  J  cr.j 


z  2Y  2(Ap)2 
0  '  Y+1  (Ap)cr 


2-  l  +  ,n{JL^}l  = 

j  l  Wm2  pi  y j 


const 

(B.24) 


then  obtain 


(Z-Z  )  =  1  +  ]  In  fl  -  -frEL 

2y  0  L  (Ap)cr,j  J  L  (Ap)2  J 
-  (Ap)z  1  en  fUEL' 

^cr.iJ  L<Ap>2J 


I 


Substituting  F.q.  (B.19)  into  Fq.  (B.25),  obtain 


B.6  Derivation  of  Eq.  (3.41):  Frozen-Shock 
Overpressure 

Generally,  the  frozen  overpressure  can  be  ex¬ 
pressed  as 


(Ap)f  (Ap) 2 


(Ap) 

P 


cr,j 

1 


(Ap)f 


2y 

Y  +  l 


(Mf2-1) 


(B.25) 


(AP)  £ 

t wr2 


(Ap) 


(Ap). 


££«2 


(3.41 
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APPENDIX  C 


PROGRAM  LISTING  FOR  RANDOM-CHOICE  METHOD 


The  program  for  solving  Eq.  (4.1)  using  the 
RCM  with  an  operator-splitting  technique  is  given 
in  this  section.  The  normalized  variables  used  for 
computation  are: 


E'  =  E/fp^Tj),  v'  =  v/.'RTj,  p*  =  p/pj, 


where  Lq  is  a  reference  length,  taken  as  Lq  =  5r0, 
in  most  calculations. 

The  time  step  is  determined  from  the  maximum 
value  for  the  local  stability  criterion  (CFL 
condition)  at  each  time  step: 


P'  =P/Pj,  T'  =  T/Tj,  oj  =  o  /(RTj), 


r1  =  r/LQ,  t’  =  ajt/C/y  LQ) 


At1  =  max [Ar '/{ |  V  |  +  v'yp'/c  ' )] 


The  program  listing  given  below  was  used  for 
the  computations  of  real-viscous  spherical  waves. 


C-l 


I 


n  n  non 


RANDOM  CHOICE  METHOD  VARIATION  H2. 

#  SPHERICAL  WAVE  « 

K  REAL  >  VISCOUS  # 

IMPLICIT  REALS6  (A-H*P-2) 

DIMENSION  KT1(10)»PT1(10),TT1(10),TV1(10) 

REAL  XARRAV  (416) » Y ARRAY ( 416 ) 

COMMON/ l K/KR1>KL> ISTP»KR»NP1, ITT,N 

COMMON//DT  >P.L>UL>PL>R>U>P>E>RR>UR>PR>Xl>Y  >GAM>SQL  >S0S  >S0R 
COMMON/OUT/TIME>DX>RHO<416) ,PRE<416) >UX(416) >ENGC416) >XR(416) 
l.PRFAC 

COMMON/RAD/ETA >REO>PRAN>TA( 416) >U2<416) 

COMMON/RELAX/SO (4 16) >UZ(416)  >El>THl>TAU 
COMMON/TSU/ISK, 1SS,  ILM 
INTEGER  TSTP 

DATA  READING 


READ ( 5 

>81) 

NPRINS 

READ  C  5 

>81) 

ISTART 

READ ( 5 

>81) 

NQQT 

READ ( 5 

>61) 

IQ 

READ ( 5 

>61) 

NSTOP 

READ! 5 

>81) 

JCT 

READ(5 

>81) 

JD 

READ (5 

>61 ) 

N 

READ ( 5 

>81) 

NHALF 

READ ( 5 

>81) 

NQQ 

READ ( 5 

>81) 

I XYP 

READ ( 5 

>81) 

I  NCR 

READ ( 5 

>61) 

ISK 

READ ( 5 

>61) 

ISS 

READ ( 5 

>82) 

TMAX 

READ ( 5 

>62) 

TMIN 

READ (5 

>62) 

PMAX 

READ ( 5 

>62) 

PM  I  N 

REAC ( 5 

>82 ) 

XP1 

READ ( 5 

>62) 

XP2 

READ ( 5 

>e2) 

XFAC 

READ! 5 

>S2) 

RMAX 

READ ( 5 

>82) 

PRFAC 

REAOI5 

>82) 

ESS 

READ ( 5 

>62) 

ETA 

REAO(5 

>02) 

WL 

READ! 5 

>82) 

PL 

READ ( 5 

>82) 

RL 

REA0I5 

>32) 

TO 

READ ( 5 

>32) 

RH 

REA0C5 

>82) 

COEP 

REAC ( 5 

>82) 

COET 

81 

FORMAT ( I  10) 

62 

FORMAT (F 15. 

7) 

C  COEFFICIENT  OF  XYPLOT 
YP1=-PMIN 

YP2=<PMAX*PMIN)/12.0 
YP3=-TM I N 

YP4=(TMAX^TMIN)/12.0 

JCTM= JCT -JO 

LMT  =  1 

NP1=N*1 

NP2=N* 1 5 

NPM=N-1 

NPX=N-5 


C-2 


I 


APPENDIX  D 


PROGRAM  OF  MacCORMACK'S  FINITE-DIFPERENCE  METHOD 


In  Section  4.3.2,  the  RCM  solutions  are  com¬ 
pared  with  MacCormack' s  solution  for  a  perfect- 
viscous  plane  wave.  In  this  section,  the  scheme 
and  the  program  of  the  MacCormack  method  are  given 
for  the  perfect-viscous  plane  wave. 

The  basic  equation  (4.1)  can  be  written  for 
perfect-viscous  plane  waves  as 


SU 

St 


SF 

3x 


S2C 

TP 


=  0 


1 

pv 

0  ■ 

u  * 

PV 

E 

,  F  = 

pv%p 
_ (E+p) v 

,  C  = 

2yv 

_  vr+uv2  _ 

E  =  0  |e  +  i  v2 


|  RT, 


p  =  pRT 


and  the  corrector  step 

Un+1  =  i  (Un+1  +  u")  -  *  (Fn+1  -  Fn+!) 
l  2  '  i  i'  2Ax  '  l  l-l ' 

-  _At  (en+l  _  2en,l  +  p+U 
2  1  l+l  i  i-l; 


The  corresponding  finite  difference  scheme  of 
the  MacCormack  method  are  the  predictor  step: 

nn+1  -  fin  _  dt  ic"  _  .  dt  r„n  ^ n  „n  ^ 

1  'Vs  lFi+l  Fi}  *  ^7  ^ci+i  '  2Ci  *  ci-i) 


(Ax) 

The  normalized  variables  used  for  computation 


E'  =  E/CpjYRTj),  v'  =  v/ / (yRTj) ,  p'  «  p/Pj, 


p  =  p/Pj,  T1  =  T/Tj ,  x’  =  x/Lq,  t’  =  ajt/LQ 


where  Lo  is  the  reference  length,  and  put  as  Lq  *  5x0. 

The  time  step  is  determined  from  90%  of  the 
maximum  value  for  the  local  stability  criterion  (CFL 
condition)  at  each  time  step: 


t'  =  max[0.9  Ax' / (/f7  +  |v'|)] 


D-l 


if  USERID  =  TN3D000  » 

if  PROCEDURE  =  L0G0N2  « 

#  TSLOG  STARTED  TIME=10J57:0fi  DATE=82-12-08  * 

»»if)fitififlf»ififififitiiiftfitifififif!tififitifififitiiififififit#«ifififif{tititifififififif»ifitif!fitit#|«#ififififif#«»a##«it*if# 


READY 

E  MAC  F  7  <  E  I  ) 

E 

LIST 

00010  C  PROGRAM  LIST  OF  MACCORMACK  METHOD 

00020  C 

00030  C  if  SHOCK  TUBE  if 

00040  C  »  PERFECT*  VISCOUS  » 

00050  C  if  MACCORMACK  if 

00060  IMPLICIT  REALMS  (A-H,P-2> 

000  70  DIMENSION  U1 ( 40 1 . 3 ) , U2( 40 1 » 3 ) , V2< 40 1 ) * X ( 401 )  »PA< 41 5 ) 

00060  1  »KT1 ( 10  >  »PT1 ( 10  ) 

OOORC  2»VC(401 )  »TC(401  ) 

00100  REAL  XARFAY(401)»YARRAY(401) 

00110  C  «  DATA  READING  if 

0012C  KJ=3P9 

0013C  ^  A  X  =  2  2  0 

00140  L  N  0  a  4 

00150  10=161 

00160  1 DC=2 

00170  \UGsO 


00150 
CO  190 
00200 
00210 
00220 
00230 
00240 
00250 
00260 
00270 
0  0  250 
00290 
00300 
00310 
00320 
0C330 
00340 
00350 
00360 
00370 
00360 
00390 
00400 
00410 
00420 
00430 
00440 
00450 
00460 
00470 
00480 
00490 
00500 
00510 
00520 
00530 
00540 
00550 


NPR  I  NT a  1 
ESS=0.03 
XFAC=2.5 
PF  AC=  1  •  0  0 
CFAC=O.ROOO 
P4=2.0000 
T4= 1 , OCOO 
wL=0.C05 

C  it  CONSTANTS  if 
KJlaKJ+1 
KJ2=KJ+2 
K J JsK  J* 1 3 
K50= I 1  9 
KL0=ID-20 
RPX=KJl-5 
NC0r4,0«NN0 
ESlaFSSKO.Ol 
CX=1.0/DFLOAT(KJ1) 

GF=1.4000 
G?  =  GF  « ( GF” 1 • 0 ) 

Gl=l . 0/G2 
VISC=1. 50-05 

RE0=( 0.10 130*04 >itWL/340.0/V!  SC 
GFPR=0.7000«(GF-1.0) 

C  #  PRINTING  OF  CONSTANTS  it 
WRITE<6»111) 

111  FORMAT(//lH0*10X,*»«»#w#»»»««»»HHHt»#»««»' */lH  »10X» 

it  •»  SHOCK  TUBE  #V1H  » 1  OX  *  • »  perfect  *  VISCOUS  »’ 

if  *  / 1 H  *  1 0  X  *  •  »  MACCORMACK  if  •  , 

if  /1H  »10X»  •  ifitififitififififif  if  ifififififififififif  it  •  ) 

WRITE (6 .112)  KJ1.NMAX.NN0* 10* I NC » NQQ . NPR I  NT, ESS.XFAC , 
if  PFAC.CFAC.P4.T4.WL. VISC 

112  F0RMAT(//1H  ,’KJ1  =  '.I3,'»  UMAX  =  *  » I  3 . 1 ,  NNO=', 

if  I  2  »  *  »  I  Da  '  .  I  3  »  '  »  I  NCs  •  .  I  2 »  '  .  NGQs  '  . 

if  13. '»  NPR  I  NT=  •  »  1 2  »/lH  »  'ESSs  •  »F7*5» ' »  XF  AC=  ’  » 

if  F5.2* ' »  PFAC=  '  »F5. 3.  •  »  CFAC='»F7.5. 

#  •,  P4  s,.F10»5.'.  T4a • .F10.5.  •  »  WLa**F7.5» 

#  /1H  , 'VISCOSITY*' ,015.7) 


D-2 
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00560  C  «  mesh  » 

00570  X( 1 )=-0.5«nx 

005B0  E'O  190  J  =  1*KJ1 

00590  190  X( J.l  )=X< JHDX 

00600  C  if  INITIAL  CONDITIONS  if 

00610  DO  206  J  =  1  »  I  U 

C  0620  U1<J»2)=P4/T4 

CObiO  U1(J»?)=0.0 

00640  V2(J)=0.0 

00650  VC( J)=C.0 

C  0  6  6  0  TC ( J ) =T4 

00670  01( J»3)=G1»P4 

00660  206  PA( J)=P4 

00690  lDl=ln+l 

00700  DU  207  J=mi,KJ2 

00710  L1(J»1)=1.0 

00720  U1(J»2)=0.0 

00730  V2( J)=0.0 

00740  VC(J)=0.0 

00750  TC ( J  )  =  1  -  C 

00760  Ul( J*3)=G1 

00770  207  PA ( J  )  =  1  •  0 

007;  o  DC  2070  J  =  KJ2*KJJ 

00790  2070  P  A  (  J  )  =  1 »  0 

00600  DO  20P  1=1*3 

00610  208  ul( IP* I )=0.5«(U1< m*l, I  )4U1 ( ID-1,  1 >  ) 

0  0  620  PA< in>*0.5#<PA( ID+1)*PA<  m-1 >  > 

0  0630  TC( ID)=C.5«(TC( I D ♦ 1 )+TC( ID-1  )  > 

00640  Y=G.C 

00650  NNi.sNKC 

00660  NCAsf.CO 

00670  C  if  PLOT  1 «  ---FOR  SLO-v  PLOTTED 

00660  CALL  DEVICE! *XYPLQT  '.OjO.O.O) 

00"  90  CALL  P  a  I  NP  (  0  *0*200  if  80*260  if  60) 

00900  CALL  VSINI (0.0,0.0*20.0*26.0) 

00910  DJ  250  I = 1  * K J 

00920  XARRAY(  I  )  =XF  ACif  (  FL DAT  (  l -1 ) /FLO *T < K J 1 ) *0 . 5«DX } 

C0930  250  YAPRAYI I )=PA( 1*1 )-l,0 

00940  YARRAY< 1  )=-0.7 

00950  CALL  PLOT<4.0»4.0»-3) 

00960  CALL  SCALE(XARRAY,12.5»KJ,1) 

00970  CALL  SCALE ( Y ARRAY , 15 . 5  ,KJ *  1 ) 

00980  CALL  AX  I  5 ( 0 . 0 » 0 . 0 »6HX- AX  I  5  * -6 » 1 2. 5 » 0 . 0 » 

00990  if  XARRAY(KJ1  )  ,XARRAY(KJ2)  ) 

01000  CALL  AXIS(0.0,0.0*12HDVERP9ESSURE, 12, 15. 5*90.0* 

01010  if  YARRAY  (KJl  )  »YARR4Y(KJ?)  ) 

01020  YARRAY (  1  )  =  1  •  0 

01030  CALL  LINE(XARRAY»YARRAY,KJ»1,C»0) 

01040  CALL  SYMPCLC  1.2*17.0»0.3»37H"ACCC!RMACK  k‘FTHDD  (  NR  =  1  60  *  TX  =  0 . 9»CFL  )  » 

01050  K  0.0*37) 

01060  CALL  SYMBOL ( 1.2* 16.0 » 0 . 3, 33HPERFECT * V I SCPUS < P41  =  2 . 00 »T4 1  =  1.0 ) , 

01070  if  0.0,33) 

01080  C  if  NACCOPMACK  if 

01090  00  209  J=1*KJ2 

01100  00  209  1=1,3 

OHIO  209  U2  (  J  *  I  >=U1(  J*  I  ) 

01120  KS1=<S0 

01130  KL1=KL0 

01140  00  360  N=1,NMAX 

01150  CFL1=1.0 

OlloO  DO  300  J=KL1*KS1 

01170  CFL2=1.C/(PSQRT(V2( J) >*DSQRT<0Afl5<PA<J>/Ul( J»l>> ) ) 

01180  IF(CFL1.LT.CFL2)  Gn  to  300 

01190  CFL1=CFL2 

01200  JCFL=J 

01210  300  C  JUT  I \UE 
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01220  UTX=CFL1*CEAC 

01230  DT=OTX*DX 

01240  D«g  =  2.0#PT 

01250  03  =  0 • 5X0TX 

01260  C4=0TX/DX/RE0 

01270  P5=D4/GFP« 

01280  06=0.5*04 

012*30  0  7  =  0.5*05 

01300  Y=Y*OT 

01310  00  302  J=KL1»KS1 

01320  DOD=OTX 

01330  U2 ( J  *  1 )=UKJ*1)-DTX«U1(J*1»2)*DPD*UKJ»2> 

01340  U2<  J *2)  =U1< J>2>-DTX«U1< J*1»1)»V2< J+l) *OPnjtUl  (  J,  1 )  »V2<  J > -PTX* 

01350  1  (PA( J*1)-PA( J) )/GF 

01360  *  +  C4»(VC(  J*1  )-2.0«VC(  JKVCI  J-l)  ) 

01370  U2<  J»3)=U1  (  J,3)-DTX*UK  J+1»?)#(UK  J  +  l  »3)*PA(  J*1)/GF  )/l  1<  J*l,l) 

01380  1  ♦ODO»Ul( J»2)»(U1( J»3)+PA( J)/GF)/U1( J.l ) 

0  1 3  <3  0  *  ♦D5*(  TC(  J*1  )-2.0iiTC(  J)+Tf(  J-l )  )  +04*  (  V?  (  J+ 1 ) -2. 0*  V  2  (  J  )  *V2  (  J-  1  )  ) 

0140C  302  CONTINUE 

01410  00  303  J=KL1.KS1 

01420  V3  =  U2(  J»2>/U2U»1  > 

01430  IF(DABS(VB).LT.(0.in-0n>)  GO  TO  3010 

01440  V2(J>=VP**2 

01450  GO  TP  3011 

01460  3010  V2 ( J  )  =  0  •  0 

01470  3011  PA( J)=G2*(U2( J»3)-0.5*V2( J)*U2( J»1 ) ) 

01480  VC< J)=2.0«VB 

014*30  TC(  J)=PA(  J)/U2(  J*1  ) 

01500  303  Continue 

01510  00  304  1=1*3 

01520  304  U 2 ( 1 » 1  )  =U2 ( 2 » I  ) 

01530  1)2  (  1  *  2  )  =  -U2  (1*2) 

01540  V  2  < 1 )  =  V  2  (  2  ) 

01550  PA  < 1 ) =PA ( 2  ) 

01560  VC ( 1 ) = VC ( 2 ) 

01570  TC(1)=TC<2) 

01580  DO  306  JsKLltKSl 

0 1 5*30  000  =  03 

01600  UKJ*1 )  =0.5*  (Ul<  J  *  1 )  ♦  U  2  <  J  *  l )  )-OOP«U2<  J»2)*03#U2<  J- 1  *  2  ) 

01610  Ul( J*2)=0.5»(01( J»2)*U2( J»2) )-00D«U2< J»1)*V?( J)*D3»U?( J-l*l) 

01620  1  »V2(  J-1)-03*(PA(  J)-PA<  J-D1/GF 

01630  »  +P6»(VC( J+l )-2.0«VC( J)*VC(  J-l )) 

01640  Ul( J>3)=0.5»(U1< J,3)*U2( J*3) >-DDD*U2( J»?)«(U2( J*3)*PA( J)/GF)/U2( J 

01650  1  *l)*n3«U2(J-l»2)»(U2(J-l»3)*PA(J-l)/GF)/U2(J-l»l) 

01660  »  +n7*(TC( J*1  )-2.0»TC( J)*TC( J-l ) ) ♦06* ( V2 ( J* 1 ) -2 . 0*V2 ( J ) +V2 ( J- 1 ) ) 

01670  306  CONTINUE 

01680  DO  307  J=KL1,KS1 

01690  V8  =  Ul( J»2)/U1< J»1  ) 

01700  IF(DABS(VB). LT . <  0 . 10-08 )  )  GO  TO  347 

01710  V2(J)=V8«»2 

01720  GO  TP  348 

01730  347  V2(J)=0.0 

01740  348  PA(J)=G2*(Ul{J»3)-0.5*v2(J)*Ul(J»l)) 

01750  VC(J)=?.0*VB 

01780  TC(J)=PA(J)/U1(J»1) 

01770  307  CONTINUE 

01780  DO  314  1=1*3 

01790  314  Ul( 1* I )=Ul<?»  I  ) 

01800  UK1»2)=-U1(1»2> 

01810  V2( 1 ) =V( 2 ) 

01820  P A ( 1 ) =PA ( 2  > 

01830  VC(1>=VC(2) 

01840  TC ( 1 ) =TC ( 2 ) 

01850  00  3145  I =KL 1  * KSl 

01860  J=KS1-I*K11 

01870  PAA=PA( J)-1.0000 
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018-0  IF(PAA.C-T.FSl)  GU  TO  3150 

Cl 890  U1(J»1>=1.0000 

01900  L'l(  J»2)=0. 000000 

01910  Ul( J*3)=G1 

0 1 9  ?  0  PA(J)  =  1. 000000 

01930  VC(J)=0. 00000 

01940  TC(J)=1. 00000 

01980  3148  M2( J)=0. 000000 

0  1960  C  B  STEP  CONTROL  R 

01970  M-'O  IF(Kll.fcQ.2>  GO  T11  318 

019,80  KLl=KLl-l 

01990  318  IF  (KSl  ,C.F  .<J1  )  GO  T H  31c 

02000  KS1  =  *M*1 

02010  316  IF(^n{»,f,Nv).':[,?)  GO  TO  *60 

02U20  1F(\. LT.Ni.fi)  C-r'  TO  8ar' 

02030  .-.RITE  (  6  » 2  0  0 1  )  r. ,  Y  ,ny  x  >  JCFL 

02040  2  00  1  F  jR*fT<  1hO»5X»?h\=,  I  ?  »  5X  ,?py  =  »F  1  *  .  5  »5X  »4Hr' TX  =  ,E  1  3 . 5 , 5  '  ,  5h  JCC  L  =  *  18  ) 

0  20  8  0  K  l=(* Sl-KLl  )/10 

02060  K^lsiYM  +  1 

02070  K,yi2  =  k  L1*KV 

0  20  3  0  Du  ?800  [=KLl»<f'-2*  fPRINT 

02090  f  J  2400  J=l»10 

02100 

021  10  2400  FT1 ( J)=PrAC«(PA(<TK J)  )-1.0 > 

02120  r.»<  ITF(6»?801)  ( KT1 ( J  )  »PT]  (J)»J  =  1»10) 

021  30  2800  COXTI’.uF. 

0  2140  28  01  FDRMST(l4»F0»4»l4*F“,4*l4»l:''.<+»I4»F4#4*i4»F-,4jI4*F>,#6»l4»F'  •  4  >  I  4  * 

02180  1  F6,4»I4»Fn.4»!4,Fu.4) 

02160  C  »  P L T  2  K 

02170  P5A=PA(UPX)-l.o 

02180  IF(PCA.GT.E5S)  GP  TO  999 

02190  DO  260  |:1.yj 

02200  YABPAYt  1  )sPA(  l  +  U-1.0 

02210  2 60  CJNTI'.UF 

02220  CALL  LINE ( X  A  9  4  A  Y »  Y  A  9  R  A  Y  »  K  J  » 1 *0*0  ) 

022  *0  ft-"1  IF(M'P(N,NCA).Nt.O)  39  TO  360 

02240  NNIXs' NN  +  NNO 

02280  f.CAs*-CA*RCB'JCO 

C2260  I  3  C  =  I  \  C  ♦  1 

02270  360  OAT  INUE 

02280  999  CALL  V$TCR*(0»0> 

02290  CALL  GPSL T M 

02300  C  «  E  if 

02310  *  RITE (6*611) 

02320  61  1  F  u  R  V  A  T  (  /  /  /  4  0  X  *  1 5  H  if  R  B  R  R  END  RRRRif/// ) 

02330  STOP 

02340  EGO 

END  OF  DATA 

INPUT 

E 

END  S 

SAVLC  TO  DATA  set  ( *T\3roOO. VAC. FORT'  ) 

READY 
TSLuG  END 

RBRRRRRRBRRRRRRR  BRR  if  If  RRR  if  if  Rif  RRR  if  RRRRRR  if  if  if  RRRR  if  RRRRRR  RBRRRRRRR  RRRRRR RRRR  R  ififitifif 
R  USFR I D  =  T’  30000  r 

B  PROCEDURE  =  Lf  GON?  B 

B  TSLOG  ENDED  T I  ME  =  1 0 : 88 ! 30  DATE  =82-12-08  B 

»  B  if  R  if  >;  if  if  B  if  if  if  if  if  B»  if  if  if  if  if  if  if  if  if  B  R  if  if  R  if  if  if  it  R  R  R  R  R  R  R  R  if  R  R  R  R  R  it  if  if  R  R  R  if  R  R  R  R  R  R  R  it  if  if  R  if  B  it  R  B  if  if  if  B  B  B 
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APPENDIX  E 


COMPARISON  BETWEEN  NEAR-FIELD  SOLUTIONS 
OF  THE  EXPLOSION  OF  A  PRESSURIZED  AIR  SPHERE 
USING  FAX,  MacCORMACK  AND  RANDOM-CHOICE  METHODS  ( RCM) 
FOR  A  PERFECT- IN V I  SC  I D  F LOW 


In  the  initial  stage  ot'  the  present  study, 
several  numerical  methods  were  tried  to  solve  the 
problem  of  the  explosion  of  a  pressurized  air 
sphere.  Some  of  the  results  are  presented  here 
to  show  the  superiority  of  the  RCM  over  other 
methods  for  analysing  shock-transitions  of  spher¬ 
ical  N-waves. 

The  near-field  solutions  using  Fax,  MacCormack 
and  RCM  for  the  same  case  as  A1  (P41  =  2.0,  Tj]  =  1.0) 
are  shown  in  Figs.  E.l,  1 .2  and  E.3,  respectively. 

In  Figs.  E.l  and  E.2  iLax  and  MacCormack  methods), 
the  time  steps  were  selected  to  be  80'i  of  the  CF1. 
condition  to  avoid  undesirable  oscillations  of 
numerical  values.  As  seen  in  Figs.  E.l(a)  and 


E.2(a),  the  Lax  and  MacCormack  solutions  give 
smoothed  shock-transitions  due  to  the  effect  of 
artificial  viscosity  in  a  rough  mesh  size  of  ;.r*  - 
1/80.  By  using  the  finer  mesh  sizes  |.'.r*  =  1/320, 
Figs.  E.l(b)  and  E.2(b)j,  this  smoothing  is  improved, 
and  the  Fax  method  gives  a  better  result.  However, 
the  smoothing  at  the  front  shock  still  remains. 

The  RCM  solutions  [Figs.  F.2(a)  and  (b)[  show 
discontinuous  shock  fronts  irrespective  of  mesh 
sizes  (.'.r*  -  1/40,  1/80),  though  some  randomnesses 
appear  in  the  expansion  part  of  a  pressure  profile. 

In  our  analysis  of  shock  transition,  it  is  necessary 
to  clarify  the  effects  of  viscosity  and  vibrational 
nonequilibrium  on  shock  thickness  without  the  effect 
of  artificial  viscosity.  Consequently,  we  adopted 
the  RCM. 


FIG.  F. .  1(h)  NF.AR-FIF.I.D  SOLUTION  OF  IMPLOSION  OF  A  PRF.SSUR  I ZHD 
AIR  SPHFRF.  USING  LAX  MET!' GO  "OR  A  PERFECT- 1 NV I  SC  I!) 
FLOW  (CAST:  Al).  MESH  SIZE  A-  =  1/320. 


SI  AR  rm.lt  SOLUTION  01  IXITOSION  01  \  I'UI  SSIIKI  :.l  I 

AIR  sn H  R!  IISINO  M.uT'ORMAOk  Ml  IIIOU  I  OR  \  IT  R!  I  01 
IWISOIU  I  l.(H\  lOASI  All.  Ml. SI!  Sill  '  r'  I  SJU. 


NI  AR-1  I !  l.I>  SOI.UTION  OF  F.XPI.OSION  OF  A  PKF.SSIIK  I  -FI) 

air  si’iifrf:  using  random-cmoicf:  mftiiop  for  a 

PFRI'FCT- 1 NV I  SC  I D  FI.OIV  (CASH  All.  MFSII  SCI!  'r*  = 


Ill)  Nl  All  -  111  II*  S()II  1 1  ION  01  I  \l’l  OS  1 1<\  01-  A  I’KI  SSIIR  I  7.1  l> 
Mil  SI’IIIRI  IISIM;  KANi'OM  OIIOIOI  Ml  1001'  1011  A 

pi  hi  re  i  - 1  nv  i  so  i  o  ii. ok  ions'  mi  si  i  si::i  r‘ 


APRLNUl X  I 


BULK  VISCOSITY  ANALYSIS  FOR  VIBRATIONAL  RF  LAX AT  I  ON  FOR  OXYC.F.N 


In  Sections  4.4.S  and  4.4.6,  the  bulk  viscosity 
concept  is  introduced  to  evaluate  the  vibrational 
relaxation  for  oxygen  instead  of  solving  the  relax¬ 
ation  equation  for  oxygen.  The  basic  equations  arc 
shown  in  some  detail  as  follows: 


dJ 

•t 


(  ii  ♦  i  ,1 1 

.‘r  1  ,  2  r  ,<r  I 
■  >r 


instead  of  Fqs.  (4.1)  and  (4.2),  where  ,.e  is  an 
effective  viscosity  including  the  bulk  viscosity 
(Uv)o  for  oxygen,  defined  by 


The  bulk  viscosity  U<v)n  is  evaluated  from  Fq . 
(3. 25): 
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C  .  0.209  exp  -  ( F . 5 ) 
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p  =  ,RT,  F  =  .i  le  *  -  v* 
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The  operator-splitting  technique  was  applied  to 
Fq.  (F.l)  as  well  as  Fq.  (4.1).  The  effect  of  vib¬ 
rational  relaxation  for  oxygen  was  taken  into 
account  in  the  step  of  viscous  correction  [Step  S; 
Fq.  (4.111)  of  the  operator  splitting  through  Fqs. 
(F.j)-(F.S).  More  precisely,  in  the  first  step, 
the  RCM  solution  should  be  obtained  by  solving  the 
Riemann  problem  for  oxygen  in  vibrational  equili¬ 
brium,  since  the  whole  flow  field  may  be  considered 
for  oxygen  as  in  quits i -equi  1  ibrium.  However,  in 
the  present  report,  the  effects  of  oxygen  vibra¬ 
tional  excitation  is  taken  into  account  only  through 
the  bulk  viscosity,  since  its  contribution  to  the 
internal  energy  specific  heats  of  the  air  molecules 
may  be  considered  as  very  small  as  long  as  it  is 
nearly  in  equilibrium  at  room  temperature.  Thus, 
the  RUM  solutions  were  obtained  by  using  the  invis- 
cid-frocen  program,  excluding  the  term  ( ’o^e  *n  F.q. 

( F .  2 )  . 
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